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ABSTRACT 


Results are reported of a program encompassing an analytical, design, 
and experimental effort to establish the fundamental feasibility of the 
fluorine lithium hydrogen tri propellant combination. This effort com- 
prised: (l) a thermodynamic-kinetic analysis of the Li, F^H system, in- 

cluding theoretical performance calculations, (i») studies of liquid-metal 
atomization md combustion, ( 3 ) design and buildup of a liquid-lithium 
(LLi) facility, (fi) design of thrust chamber, including a fluorine-rich, 
hydrogen-fluorine gas generator and two lithium injectors, and ( 5 ) experi- 
mental firings of I F 2 /LLi GIF, The test matrix covered a range of chamber 
pressure (350 to 520 psia), F, Li mixture ratio (2.4 to 5.2), and propor- 
tion of hydrogen (15 to hk percent of total flowrate), and included 
variations in point of hydrogen addition and in chamber length. Corrected 
characteristic velocity efficiencies were 95 *, 100 percent, showing that 

use of this tripropellant combination is feasible and that highly efficient 
combustion can be achieved. 
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SUMMARY 


This report presents the results of an analytical and experimental invest! 
gation of the combustion ch .acteristics of the Li/F,,/]^ tripropellant 
combination. The primary overall program objective was to determine 
whether the high theoretical performance of the Li/F 0 /lf,combination is 
experimentally attainable and to investigate the necessary design require- 
ments. Characteristic velocity (c*) efficiency of 95 to 100 percent 
was obtained in experimental tripropellant firings, with injector and 
combustion chamber designs based on the results of analytical and experi- 
mental studies of the atomization and combustion of liquid metals and of 
means of carrying out these processes efficiently. Nominal design point, 
parameters were as follows: 

1. Chamber pressure, 500 psia 

2 . mixture ratio, 2.7k (stoichiometric) 

3. Hydrogen, 15 to ^lO percent of total propellant flowrate 

'*• Nominal altitude thrust (30 percent H Q , € = 60), 2900 pounds 

The program was divided into two major tasks: 

I. Analysis and design 
II. Combustion performance tests 

The initial effort of Task I was an investigation of the chemical charac- 
teristics of the Li/Fy/Hy system including an examination of the kinetics 
of the chemical reactions involved and detailed theoretical performance cal 
culations.. The kinetic studies showed that gas phase reactions of fluorine 
with hydrogen, fluorine with lithium, and lithium with hydrogen fluoride, 
even at the lowest temperature likely to exist in the combustion chamber, a 
more rapid by several orders of magnitude than the physical processes of 
propellant mixing and lithium vaporization. Hence combustion efficiency 
will be significantly affected by physical processes rather than by 
chemical kinetics, as is usually the case in rocket engine combustion. 


Theoretical performance calculations were made for LF 0 /LLi*in the mixture 
ratio range 2.20 to 3 . 30 , at chamber pressures of 200 to 1000 psia and 
with gaseous hydrogen addition ranging from 0 to 40 percent of total pro- 
pellant flowrate, as well as for U^/SLi^ under similar conditions. 
Calculations were also made to determine the effects on performance param- 
eters of incomplete condensation of LiF in the combustion products. 

Results of the computations are presented in a series of curves which show 
the variations in important performance parameters (i , c*, T , Cp) with 
proportion of hydrogen, chamber pressure, mixture ratio, and type of 
expansion. For LF^/hLi , the maximum vacuum specific impulse at 500 psia 
chamber pressure is 458 lbf-sec/lbm at € = 60. For LF 0 /LLi/GH, } , the max- 
imum vacuum specific impulse, which occurs at or near stoichiometric 
^/Li mixture ratio and with approximately 30 — percent hydrogen, is 
542 lbf-sec/lbm at 500-psia chamber pressure and € = 60 . The maximum 
vacuum specific impulse for LF^SLi/u^, under the same conditions, is 
520 lbf-sec/lbm. Chamber pressure variation (200 to 1000 psia) has little 
effect on vacuum specific impulse. 

The second Task I analysis and design effort involved establishment of 
requirements for high combustion efficiency in the LF^/LLi/GF^ system. 

It consisted of the following two parts: (l) on analytical study of 

lithium combustion based on a one-dimensional, vaporization rate-limited, 
combustion model, and (2) experimental investigations of liquid-stream 
penetration into gas jets and of the gas-augmenied atomization of a 
liquid metal . 


Because of the low burning rate of liquid lithium in gaseous fluorine 
(about one-tenth that of more common propellants), provision of very 
small lithium droplets is an essential requirement for high combustion 
efficiency. An extensive series of calculations was carried out to 


*The abbreviations LLi and SLi 


refer to liquid and solid lithium, 


respectively. 
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establish the effect on combustion efficiency of variation in lithium 
aroolet size, as a function of characteristic chamber length, chamber 
pressure, and point of hydrogen addition. It was found that a lithium 
droplet volume mean diameter of 20 microns would provide c* efficiency 
m excess of 0.96 with an F^/Li combustion zone L* of 40 inches; 50-micron 
droplets require an L* of 180 inches for the same result. Curves showing 
the calculated effects on combustion efficiency of varying the axial point 
of hydrogen addition were also derived. It was found that for maximum c* 
efficiency the hydrogen should be injected downstream of the point of 
fluorine-lithium injection. However, this effect was found to be very 
small for small diameter lithium droplets. 

Consideration of the atomization characteristics of liquid lithium showed 
that mean drop sizes on the order of 20 to 50 microns could be practically 
achieved only by gas-augmented liquid atomization. To aid in establishing 
specific design crixeria for a gas-liquid system, a photographic study 
was made of the penetration of gas jets by liquid streams in a triplet 
element, using nitrogen and water or liquid Cerrosafe (a proprietary, 
low-melting-point, lead-bismuth alloy which simulates lithium). 'These 
experiments indicated that the degree of penetration is a function of 
liquid stream diameter, liquid/gas velocity ratio, liquid/gas density ratio, 
and the angle between the liquid stream and the injector face. A satisfactory 
empirical correlation was established. 

Measurements were made of the particle sizes obtained by atomization of 
liquid Cerrosafe by gaseous nitrogen. Observed volume mean droplet size 
was adequately expressed by the Wol fe-Andersen equation, modified to 
account for the reduction in gas momentum due to acceleration of the 
atomized droplets. It was also confirmed that a conditional parameter 
based on the Weber and Reynolds numbers must be met to obtain the most 
efficient atomization. 
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A third effort conducted during Task I consisted cf the experimental eval- 
uation of two systems which are basic to the tripropellant thrust chamber 
design concept employed. The first was a high pressure gas generator in 
which liquid fluorine is reacted with gaseous hydrogen at high mixture ratio 
to produce fluorine gas at an elevated temperature. Its design consisted 
of separate fluorine and hydrogen injectors and a two-segment chamber, with 
each segment being 5 inches long. Firings of this assembly were carried 
out to verify that combustion efficiency in the gas generator was close to 
theoretical and that product gas at the desired temperature and pressure 
was consistently obtainable. 

The second system evaluation was directed to the design, fabrication, and 
operation of a facility to store, heat, and deliver liquid lithium, fol- 
lowed by liquid lithium flow studies. The system was designed, assembled, 
and successfully operated. Test flows of liquid lithium through orifices 
of various sizes and orientations indicated that minimum orifice size for 
reliable free flow of lithium was 0.043-inch and that no face adherence 
was obtained with orifice-injector face angles of 45 degrees or larger. 

The final Task I effort consisted of the design of the F 0 /Li/H 0 thrust 
chamber assembly on the basis of the criteria developed. The concept 
consisted of the following: 

!• f 2 / H 2 ® as generator operating at a nominal mixture ratio of 

about 500:1 to provide hot (960 R) gas for liquid lithium atomi- 
zation 

2. Lithium injector containing showerhead F 2 orifices which produced 

gas jets into which high-velocity liquid lithium streams were 
injected 


k 


Primary F 0 /Li combustion chamber 
Hydrogen injector 

F^/Li combustion product-hydrogen mixing ehambei 
Convergent-divergent nozzle 


Tho entire assembly was uncooled and most components were graphite or 
graphite-lined. 


The Task II (combustion performance tests) effort consisted of a series 
of experimental firings to establish the combustion characteristics of 
the F 2 /Li/ll 2 tripropellant combination. The results of these tests 
defined performance as a function of percent hydrogen added, point of 
hydrogen addition, and chamber length. The configurations employed ere 
sketched below: 





The experimental results clearly demonstrated that high combustion effici- 
encies can be achieved with the F,/u/ll. t tripropellant combination. Char- 
acteristic exhaust velocity efficiencies of essentially 100 percent were 
obtained with all three chamber configurations. In the last two at the 
lowest hydrogen flowrates (l 5 to L>0 percent of total propellant flowrate), 
c* efficienc y decreased to approximately 95 percent because of nonunion/ 
distribution due to incomplete penetration of the F Q / L i combustion gas 
stream by the hydrogen. 
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INTRODUCTION 


An analytical, design, and experimental program has been conducted to 
determine the combustion characteristics of the l.i/l' () /H ( , propellant com- 
bination. The combustion of fluorine with a light metal such as lithium 
releases an extremely large amount of heat per pound; if this thermal 
energy is used to beat a very light working fluid such ns hydrogen, a 
significant increase in specific impulse over that provided by bipropel- 
lant combinations is theoretically avai lable ( TTO-second theoretical 
specific impulse values are predicted). The objective of this program 
was to determine experimentally the combustion efficiency actually attain 
able with these propellants to ascertain the extent to which the desired 
reactions occur. A 13-month program, divided into the following two 
tasks, was completed: Task 1 was an analytical and design effort to 

define injector-thrust chamber design criteria for high performance; 
lask II consisted of experimental combustion performance testing. 

During the past sevcm years a significant an.ou.it of basic and applied 
combustion research has been carried out in the field of light-metal fuel 
additives and extended to the application of heavy metals as well. Metal 
additives have become standard components of various sol id-propcl l.Jit 
formulations, some of whim now prof i tartly use aluminum in quantities up 
to 110 percent by weight. In the liquid propellant field, effort is being 
oriented toward the use in the Titan II engine of aluminum in a gelled 
suspension with hydrazine in place of the 30-30 mixture of hydrazine 
(N L ,II ri ) and unsvmme t r ica 1 dimethyl hydrazine (UD.'UI). In this application, 
there -s a 3-percent increase in performance and. more important, a sig- 
nificant increase in propellant density. In addition, the t r: propel 1«„, 
system lie, 0^ II,, has recently been studied under Air Force sponsorship. 
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In general, the previous multicomponent luel studies were conducted to 
evaluate overall engine performance. No serious attempt was made to 
define injector design criteria for high combustion efficiency independent 
of potential performance deficiencies due to nozzle flow kinetics. Con- 
sequently, if overall performance were low, the causes of the overall 
inefficiency could only be postulated. 

The basic premises of the present program were as follows: (l) for maxi- 

mum performance with the Li/F^ system, it is necessary that the lithium 
and fluorine be completely reacted prior to addition of hydrogen; ( 2 ) injec- 
tion of hydrogen into a stream of LiF gas and subsequent attainment of 
thermal equilibrium do not present major problems; and therefore, ( 3 ) 
the most important consideration is the design of injection methods which 
vill result in complete Li-F^ reaction. The major portion of this program 
was consequently devoted to the investigation of means whereby high Li-F () 
combustion efficiency could be achieved. 

A model for metal combustion (Ref. 1 ) was used which agrees conceptually 
with available hydrocarbon droplet experimental results. This model suggests 
that liquid metal combustion is basically a vaporization rate-limited 
process, complicated by the formation of a metal oxide (or fluoride) on 
the metal surface. rhis oxide layer tends to reduce liquid metal vapori- 
zation although it is considered porous. In general, the problem is 
aggravated if the initial drop or particle size is large, if the metal 
liquid vapor pressure is low, and if the oxide or fluoride melting tem- 
perature is significantly higher than that of the metal itself. Two 
distinct classes of metals, volatile and nonvolatile, were considered in 
Ref. 1 . The metals A1 and Be, for example, are in the latter category. 
According to Ref. 1 , "the oxides of A 1 and Be are protective materials 
which adhere well to the metal surface. It is expected that an oxide 
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layer on the surface of these metals largely inhibits oxidation over a 
\*ide temperature range. Uhen the surface temperature cannot rise above 
the melting point of the oxide, the solid shell thickens and further oxi- 
dation is inhibited." Ihtis. combustion of nonvolatile metals which have 
oxide melting temperatures higher than the metal boiling point would be 
very slow. Conversely, volatile metals such as lithium are characterized 
as follows (Ref. 1 ) : "All of these metals have a vapor pressure of about 

1 mm Ilg or more at 1000 K....In view of the criterion for vapor-phase 
burning, it is expected that these metals can all burn in diffusion 
f lames .... Their melting points are low compared to realistic propellant 
flame temperatures. Accordingly, these metals should burn as liquid drop- 
lets. It is further pointed out that if large drops (80 to 100 microns) 
exist, some diiliculty may arise due to oxide formation on the surface 
during the time when the comparatively large mass is being heated. 

these studies indicated that efficient reaction of liquid lithium with 
fluorine in a rocket combustion chamber is feasible, and that its success- 
ful realization depends primarily upon the use of properly designed 
injection systems. To develop such systems, the same basic approach to 
the successful definition of injector design criteria was followed as 
that which has been thoroughly proved in previous performance programs 
(c.g.. Ref. J and "5 ) . Only the details of application differed. In 
♦his approach the fundamental injector spray parameters, propellant dis- 
ti ibution and atomization, are considered to be performance controlling. 
Ihus. lor the 1 luor ine 1 i t Ilium combination, combustion efficiency wi'l 1 be 
a function of local fluorine and lithium mass and mixture distributions, 
ns well as liquid lithiiun drop size. Control of the latter is considered 
the most critical spray requirement in the fluorine-lithium reaction. 

Data presented in this report indicate the necessity of achieving mean 
lithiiun droplet sizes no larger than about U0 microns for efficient 
combiis t i on . 
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Consideration of the atomization characteristics of liquid lithium showed 
that simple impingement of liquid streams of practical size and velocity 
would not produce drop sizes small enough for realistic engine designs 
because of the extremely high lithium surface tension. Therefore, injec- 
tor design concepts which achieve liquid atomization by use of high- 
velocity gas were considered. In principle, this is identical to the con- 
cept that was successfully employed in the NASA-sponsored F 2 /H 2 nozzle 
program, in which efficiencies approaching theoretical were achieved 
over a mixture ratio range of 9 to 15 and a chamber pressure range of 
50 to 200 psia (Ref. 2). In that case, performance was considered to be 
liquid fluorine vaporization-rate limited, and liquid fluorine atomization 
was accomplished by use of high-velocity hydrogen gas. 

Lithium atomization by high-velocity hydrogen was not considered desirable 
because this would result in unfavorable propellant distribution, i.e., 
lithium drops in a hydrogen atmosphere. Instead, the concepts studied 
provide for generation of high-velocity F 2 and HF gases for lithium atom- 
ization. During Task I (Analysis and Design), cold-flow experiments were 
conducted with simulating fluids to define qualitatively the gas forces 
necessary to achieve 10- to 50-micron drop sizes. Hot-gas generator experi- 
ments were also carried out during this task to investigate the capability 
of delivering the hot-gas force necessary for atomization. These studies 
provided sound bases upon which injector design could proceed. 

The Task II effort (Combustion Performance Tests) included determination 
of hot-firing c* efficiency as a function of percent hydrogen added (15 to 
40 percent of total flow), hydrogen injection station, and chamber volume. 
The results obtained clear.’y demonstrated the feasibility of obtaining 
high combustion efficiencies with this tripropellant system. 
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The Task I (Analysis and Design) and Task II (Combustion Performance Tests) 
efforts are reported under separate headings within the text. The Task I 
discussion is subdivided as follows: 

** Thermodynamic and kinetic analyses of the Li/p ^ system 

2. Analytical evaluation of Li/F^ combustion 

3. Liquid-metal atomization 

4. Thrust chamber concept and design of components 

5. Experimental evaluation of fluorine-rich gas generator 

6. Liquid lithium facility and flow characterization 

The Task II discussion is divided as follows: 

1. Experimental facilities and procedures 

2. Experimental evaluation of Li/l’/H combustion 



TASK I 


THERMODYNAMIC AND KINETIC ANALYSES 
OF THE Li-F f) -H () SYSTEM 


REACTION KINETICS 

A study was made of the kinetics of the various competing reactions involved 
in the Li-F^-H 2 system to determine the sensitivity of combustion efficiency 
to the reaction rates. For purposes of this study, it was assumed that all 
three components are initially present and interacting. Under these condi- 
tions, it is important to ascertain to what extent lithium and hydrogen will 
compete for the fluorine and whether lithium would react rapidly with HF. 
if this were the first product formed. 


Hydrogen/Fiuorine Reaction 


Because hydrogen and fluorine ignite hypergolical ly under confined condi- 
tions (Ref. 'i), it might be expected that the reaction of gaseous hydrogen 
and fluorine would be quite rapid. However, Levy and Copeland (Ref. 5 ) 
have shown the reaction of purified gaseous fluorine and gaseous hydrogen 
to be slow at +230 F. with a half life of 32 seconds. They found the 
reaction rate under these conditions to be independent of the hydrogen 
concentration and proportional to the fluorine concentration (i.e. , first 
order in fluorine). Brokaw (Ref. 6) has proposed the following energy 
chain mechanism to explain the results of Levy and Copeland: 


F + H Q — * HF + H 


H + F 0 
HF* + F f> 


allF* + (l-a) IIF + F 


( 1 ) 

( 2 ) 

(3) 


HF + 2F 
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HF* + M 


(4) 

(5) 




— ► HF + M 
2F + M — F + Ni 

An energy chain mechanism has also been proposed by Kapralova, et al. (Ref. 

38 ), whose work on fluorine/hydrogen ignition limits seems to support this 
nostu late . 

The kinetic results of Levy and Copeland do not suggest autocutalysis by 
the product HF, although Grosse and Kirshenbaum (Ref. 7 ) found it neces- 
sary to eliminate HF (and transition metal surfaces) to prevent spontaneous 
ignition of mixtures of hydrogen and fluorine at low temperatures. 

To extrapolate the kinetic results to temperatures of interest in the 
present study, it is necessary to know the activation energy. Levy and 
Copeland (lief. 8) present data which give a half life at 162 F of 
li minutes. This half life and their earlier value of 32 seconds at 
230 F (Ref. 3 ) give an activation energy of 17 .ii Kca 1/mole. Brokaw’s 
mechanism predicts an activation energy, , equal to 2 Ej+E -E -E ; (Ref. 6 ). 
Available estimates of E } are 1.3 Kcal/mole (Ref. 9 ) and 5^7 Lai mole 
(Ref. 10). Because the other three reactions appear to have very small 
activation energies, a small overall activation energy on the order of a 
few kilocalories is predicted from Brokaw's mechanism. However, since 
the rate constant for Reaction 3 represents an integrated average over 
oil the vibrational energy levels of HF*, the activation energy of 
Reaction 3 may actually be large enough to permit agreement between 
Brokaw’s mechanism and Levy’s activation energy. 

The lowest temperature encountered in the Li-F„-H (> system would be that 
occurring if the hydrogen and fluorine were to react to equilibrium, with 
the liquid lithium remaining unreaiteu and, together with the excess 
hydrogen, acting as a heat sink. This would result in a theoretical 
temperature of about 1400 K, with the largest proportion of hydrogen 
being considered (40 percent). Using an activation energy of 17.4 Kcal 
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to extrapolate the results of Levy and Copela-d, the predicted half life 
at 1400 K is 2 microseconds. It thus appears that the hydrogen/fluorine 
reaction will he quite rapid even under the most adverse temperature con- 
ditions likely to be encountered in the combustion chamber. This would 
not have been the case had the activation energy been found to be much 
smaller. For example, a value of 3 Kcal would have led to a predicted 
half life of 0.3 seconds. 

It is of interest that lithium vapor would be expected to inhibit the 
" J «rogen, fluorine reaction, which is apparently a chain reaction in which 
fluorine atoms are one of the chain carriers. If an appreciable amount 
of lithium should vaporize before the hydrogen/fluorine reaction is com- 
plete, the fluorine atoms would react more readily with lithium atoms 
than with hydrogen molecules, thereby inhibiting the chain reaction. It 
is expected, however, that the reaction of hydrogen with fluorine would 
reach completion before the lithium is vaporized. 


Lithium/HF Reaction 


The above discussion suggests that the gaseous lithium atoms will find 
themselves in an atmosphere of ItF and excess hydrogen os they evaporate 
from the surface of the parent liquid lithium droplets. It is of impor- 
tance, therefore, that they react readily with HF if the desired reaction 
is to go to completion within the combustion chamber. An attempt was 
made to estimate the minimum rat‘e of the reaction: 

Li(g) + HF ► UF(g) + H 4H = -1.0 Kcal/mole (6) 
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Exchange reactions of this type normally have activation energies on the 
order of 5 to 10 Keal. Using the empirical rule of Semenov (Ref. li), an 
activation energy of 11.2 Kcal/mole is obtained. At U00 K, an error of 


5 Kcal in the activation energy causes an error of only a factor of six in 
the estimated rate. Estimation of the p-e-exponential factor is a more 
likely source of error, hut most exchange reactions have pre-exponential 
factors between IQ 8 ! 1 ' 1 and lO 10 ? 1 2 (mole- sec)" 1 
value of should be approximate! v 


Therefore , the minimum 


8 1 

= 10 T exp(-ll ,000 ItT) (mole-sec) 


-1 


At I'iOO K this estimated rate expression will give a value of 

7.0 x 10' (mole-sec)" 1 for k . 

0 

At chamber conditions of MOO K. 500 psi, and 20 mole percent Iff, the Iff 
concentration would be 0.058 mole liter. In this cose, the maximum pre- 
dicted half life of a lithium atom would be 0.2 microseconds with the 
value of k^ estimated above. Even though the Iff concentration will be 
much smaller as the reaction nears completion and a number of half lives 
are required to approach complete combustion of the lithium, this result 
indicates that the reaction of lithium with Iff will be virtually complete 
m less than 10 microseconds even under the most adverse conditions which 
might be obtained (note that as the Iff concentration decreases, the 
temperature and reaction rate increase). 


It is believed that the lowest likely rate constant was used in calculating 
the rate of Reaction 6. If this reaction should have a very abnormally 
low rate for some reason, the reaction half life will be correspondingly 
longer, but this is very unlikely. 
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Fluorine, Lithium Yauor Reaction 

The kinetics of the reaction of lithium vapor with 
studied . 


fluorine has not been 


Li v g) + ► Lii + I* 4 II = -09 Keal/mole ( 7 ) 

Empirical formulas for estimating activation energies (Ref. 12 ) predict 
virtually no activation energy for Reaction 7. since a relatively weak 
bond is being broken and the reaction is highly exothermic. Assuming a 
stenc factor at the low end of the "normal" range of steric factors for 
reactions of this type, and following the same reasoning used above for 

the lithium/ 111 ’ reaction, the minimum rate constant for Reaction 7 should 
be approximately 


k-, - 10 T (mole-sec ) ^ 

This assumes that about one collision in 100 is effective in leading to 
chemical reaction. At H00 K and a fluorine concentration of 0.03 mole'' 
liter, this predicts a lithium atom half life of only 0.006 microseconds. 
It thus appears that if the fluorine does not react rapidly with hydrogen 
for some reason (such as inhibition by Li atoms), the reaction of lithium 
atoms and fluorine will he very rapid even at low temperatures. The rate 
of this reaction will be relatively insensitive to temperature. 


Condensed Lithium Reactions 

A direct reaction between fluorine or Iff and liquid lithium may occur at 
the surface of the lithium droplets. This would only be likely hefore 
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the droplet lias ignited or been heated to a temperature at which it vol- 
atilizes rapidly. Otherwise, the reaction of lithium vapor will markedly 


deplete the concent rax i un of 


oxidizer at th< 


surface of the d r on 1 pfs. 


Reaction at the liquid surface would contribute to self-heatinc: of the 
surface . 


No experimental data are available on the rate of reaction of liquid 
lithium with HF or fluorine. It is anticipated that HP 1 and particularly 
fluorine would react moderately rapidly at the liquid surface. However, 
the LiF layer that would form melts at 1143 K and would not be expected 
to reduce the rate of lithium reaction above this temperature. 


Conclusions 

For the concurrent reactions of gaseous fluorine, gaseous hydrogen, and 
liquid lithium, on the basis of available and estimated rate constants, it 
is concluded that the hydrogen/fluorine reaction would be very rapid even 
at the lowest temperature likely to exist in the combustion chamber, with 
predicted half life of about 2 microseconds at 2520 R (1400 K). As the 
lithium vaporizes, it will react extremely rapidly with the HF (predicted 
half life of about 0.2 microseconds at 2320 R) to form LiF. If the lithium 
vaporizes before HF formation is complete, the lithium may inhibit further 
HF formation; under these conditions lithium vapor will react nearly 
instantaneously with fluorine to form LiF (predicted half life of about 
0.006 microseconds at 2520 R). In any -event , the equilibrium condition 
following lithium vaporization is that the fluorine is nearly completely 
reacted with the lithium rather than with the hydrogen. Further, within 
the certainty of the estimated rate constants, the desired reactions will 
proceed much more rapidly than the physical processes of mixing and lithium 
vaporization, so that the latter will be the limiting factors in the 
attainment of complete combustion; reaction kinetics will not significantly 
affect combustion efficiency. 
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Tho experimental procedure adopted in the present investigation was to 
permit the fluorine and lithium to react prior to addition of hydrogen, 
at combustion zone temperatures of about 0800 R. However, although 
addition of all three propellants at the same axial position would result 
in substantially lower theoretical combustion temperatures (loOO to >'t00 R, 
depending upon the proportion of hydrogen), they are still high enough so 
that combustion efficiency will not be affected by reaction kinetics. 

THEORETICAL PERFORMANCE ANALYSES 

Theoretical performance calculations for the Li F„ and Li F, ( H t> systems 
were made with the Rocketdyne N- Element computer program (Ref. ll), in 
which the propellants are reacted isentha lpica 1 ly at a specified pressure 
to form atomic and molecular species in thermal and chemical equilibrium, 
defined by minimizing the (iihbs free energy of the system. The most 
recent JANAF propellant property tabulations are used for input data. 

The combust ion process is irreversible, and therefore non- i sent ropic , hut 
all other nonequi 1 ibrium effects (chemical kinetics, transport processes, 
wall effects, etc.) are excluded. After combustion, the gaseous products 
expand by ideal onc-d imens i nnn 1 . isentropic flow through a nozzle to a 
specified exit pressure. Continuous chemical, thermal, kinetic, and 
phase equilibria are maintained among all species and phases. Temperature, 
pressure, density, velocity, and composition are uniform across any section 
normal to the velocity gradient. The nozzle is assumed to be fed from an 
infinite, stagnant , chamber and to discharge in parallel flow, so that 
calculated specific impulse is maximized. 
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In addition to the shifting equilibrium calculations, the model permits 
assumption of frozen composition; that is, chemical composition may be 
frozen in the chamber or at any point in the expansion, downstream of 

vh.ch the combustion product composition and phases remain invariant at 
the freezing point values. 

Performance calculations were made to cover the LF 0 LLi mixture ratio 
range of 2.0 to at chamber pressures of 200 to 1000 psia, with the 
addition of various amounts of hydrogen ranging from 0 to >i0 percent of 
the total tripropellant flowrate. Limited calculations were also made 
With solid lithium injected at 36 R with liquid hydrogen. Auxiliary- 
calculations included some in which no condensation of LiF was permitted, 
and others in which only extremely small amounts of hydrogen were present, 
as would be the case in experimental Li > o firings (whore F /h i s 
employed in the gas generator at very high mixture ratios)." 

Results of the theoretical performance computations are presented in a 
series of curves which show the variations in important performance 
parameters at specified conditions of chamber pressure, mixture ratio, 
and type of expansion. The data are then compared to establish the effects 
on the various performance parameters of chamber pressure, mixture ratio, 
proportion of hydrogen present, type of expansion, propellant state at 
injection, and incomplete condensation of LiF. 


LF, , LLi Performance Calculations 

Results of the first series of calculations (Table 1 ), for liquid fluorine 
at 153 P- and liquid lithium at 060 R, are illustrated in Fig. 1 through >, . 
The range of chamber pressure is 200 to 500 psia, and the range of F r /Li 
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mixture ratio is L’.O to 5.4. Maximum theoretical c* occurs at a mixture 
ratio of approximately 2.2 at 500 psia and 2.0 at 200 psia. On the other 
hand, chamber temperature, vacuum specific impulse, and thrust coefficient 
reach maxima at. or very near, the stoichiometric mixture ratio (2.74). 

As shown in Fig. 5. maximum vacuum specific impulse at 500 psia chamber 
pressure is 450 lbf-sec/lbm at an expansion ratio of 40 and 488 lbf-sec/lbm 
at an expansion ratio of 500. 


Previous calculations of theoretical performance parameters for the liquid 
1 it hi urn/ liquid fluorine propellant combination were published by NACA in 
1051 (Ref. 14). Inthalpies of formation used for those calculations were 
-5.050 Real mole for LF n (153 R) , and 2.305 Kcal/mole for LLi ( 1000 R) . 

Both frozen a:id shifting equilibrium performances were computed, for 
500 psia chamber pressure and expansion to 1 atmosphere. Two sets of 
computations were made: one in which the combustion products (Li, F, LiF, 

Li > 1 , e ) included ionized substances, and the other in which only 
nonionized combustion gases were considered (Li, F, LiF). It was found 
that the effects of ionization on the performance parameters are negligible. 


Propellant inlet conditions for the performance calculations of the present 
study (Ll^ at 153 R and LLi at 960 R) were comparable to those of Ref. 14. 
The combustion products considered were Li, F, LiF, Li 0 , and Li Q F 0 . 

Although the last two species were no! included in the computations of 
Ref. 14 , they are present in only very small amounts (mole fractions less 
than 0.<>l), so that the inlet conditions and combustion products in both 
sets of calculations are comparable. It is therefore possible to compare 
the results obtained (Fig 5 and (>). Shifting equilibrium c* and I 

8 

reported in Ref, 14 are somewhat higher than those calculated in the 


present study, but values of frozen composition c* and C p in both compu- 
tations are very close. Since the present calculations are based on newer 
thermochemical data and were obtained with the most recent modification of 




21 



i 


/ 


1 




An amount of hydrogen equivalent to F 0 /lI () mixture ratios of 100 to 100 vas 
expected to he used for reaction with fluorine in the gas generator upstream 
of the lithium injector in the experimental U/F 0 firings of this investiga- 
tion. Therefore, theoretical performance calculations were made which 
included these proportions of hydrogen to provide a base for comparison 
o* experiments 1 results. It vo> found that inclusion of such small amounts 
of hydrogen resulted in very minor changes in the performance parameters. 

The changes, for Li mixture ratio of 2.7'» and chamber pressure of 
210 D«:a, are shown in Fig. 7 and amount to about 0.1 percent for c* 

and I . 
s 



Performance Calculations 

An extensive series of performance calculations was carried out for the 
R)/LLi(960 R)/0H^(l37 R) tripropellant combination. The various 
parameters are illustrated in a series of curves which are listed in 
Table 2 and shown in Fig. 8 through Vi. Important results are 
summarized as follows: 


1. At 500 psia chamber pressure and with 11 to '»0 percent added 
hydrogen (Fig. 8 ), vacuum specific impulse is maximum at, or 
very near, the stoichiometric F f) Li mixture ratio (2.7'j); varia- 
tion of specific impulse with mixture ratio above the stoichio- 
metric point is more pronounced at the higher hydrogen levels. 
Optimum percentage hydrogen, which varies with mixture ratio and 
expansion ratio, is 10 to 11 percent at e = >»0 (Fig. 8 ), de- 
creasing to 2h to 28 percent at € = 100 (Fig. 10 ). Maximum 
vacuum specific impulse at 100 psia chamber pressure is 111 lbf- 
sec lbra (l’ 2 I.i MR = 2.7'*, 11 percent) at * = '.(> (Fig. 8 ), 
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")/ 1 lbf-sec lbm (i - ,, I.i Mil = 2.7't. II,,: 10 percent) at c = 100 
(lift. «)). and >78 llf-sec lbm ( F, , I.i Mil - 2.7't. H (> : 28 percent) 
nt € = 100 (Fiir. 10). Chamber pressure variation in the rnnfte of 
'tOO to 1000 psia has a small effect on vacuum specific impulse at 
€ = 10 (Fi<t. 11 ) and slifthtly lnrfter effect at € ~ 100 (Fijt. 1>); 
in either case, the chanites at the higher percentages of liydroften 
are nearly negligible. Differences between shifting equilibrium 
and frozen exposition values of vacuum specific impulse decrease 
with increasing percentage of hydrogen (Fig. 22). At mixture 
ratio 2.7't and €= 10. the approximate difference is 11 lbf-sec ' 
lbm at 21 to 10 percent hydrogen. 17 lbf-sec lbm at 11-percent 
hydrogen, and 10 lbf-sec lbm at 'tO-perccnt hydrogen. 

2. Characteristic exhaust velocity ( c *) peaks at stoichiometric 
mixture ratio with the higher percentages of hydrogen (21 to 
10 percent) hut not at H-perccnt hydrogen (Fig. 21). At 100 psia 
chamber pressure, mixture rai io 2.7't. and 10-percent hydrogen, 
maximum c* is <>101 ft sec with shifting equilibrium (Fig. 2l). 
and 0011 ft sec with frozen composition (Fig. 2l); corresponding 
maxima nt 1000 psia chamber pressure are 0211 ft sec (Fig. 2(>) 
and 0118 ft sec (Fig. 27). 

1. Similar to o + . chamber temperature also peaks at the stoichio- 
metric mixture ratio with the higher percentages of added hvdro- 
gen (10 to ’*0 percent) but not with lesser amounts (Fig. 28). 

The substantial heat absorption capability of the hydrogen is 
shown by the chamber temperatures at 2.71 Mlt and 100 psia 
chamber pressure in the absence of hydrogen (0811 R. Fig. 2 ), 
with 1 1-porcent hydrogen ( !'.!<> li, Fig. 20), and with 'iO-percont 
hydrogen (ll2(> It, Fig. 20). 
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h. Vacuum thrust coefficients ((’p) at expansion ratios of 40 to 500 
(Fig. 31, • 33) are generally maximum at the stoichiometric 

mixture ratio with 25 to 53 percent added hydrogen, but not out- 
side this hydrogen addition range. As with specific impulse, the 
differences between shifting equilibrium and frozen composition 
values of vacuum thrust coefficient are more pronounced with 
15- to 30-percent hydrogen addition than at higher hydrogen 
percentage levels (Fig. 40). 

5. The mass fraction of condensed species in the chamber increases 
with increasing percent added hydrogen, with decreasing mixture 
ratio, and with increasing chamber pressure (Fig. 41 , 42 ). In 
the nozzle, condensation is complete at an expansion ratio that 
decreases with increasing percent hydrogen; for example, with 
- ) percent added hydrogen, at 500-psia chamber pressure and 
2.7 ; » F 2 /Li mixture ratio, maximum condensation (75 weight percent) 
is reached at € = 300. but with 40-percent hydrogen at the same 
chamber pressure and mixture ratio, maximum condensation (60 
weight percent) is reached at € = 3 (Fig. 43). 


A limited number of computations were made to estimate the effect on thrust 
chumber parameters of incomplete LiF condensation in the combustion chamber. 
Decreasing the amount of heat available for raising the hydrogen temperature 
would result in lowered chamber temperature, characteristic exhaust velocity 
and specific impulse. The effect is shown in Fig. 45 and 40. The abscissa 
represents the actual weight percent of condensed I.iF as a fraction of the 
theoretical condensed weight percent. These effects are largest with 
40-percent hydrogen addition (with 1 j- to 25-percent hydrogen, there is no 
LiF condensation in the chamber at 500 psia chamber pressure and 2.74 F /Li 
mixture ratio) Tn the extreme case of no LiF condensation, with 40-percent 
hydrogen, chamber temperature decrease would be approximately 8.5 percent, 
c* decrease, approximately 4.4 percent, and I (vac, € - 40 ) decrease, 
about 6,2 percent. 
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IZ. ^SIi, jJ| 4 > _P_erf ormance Calculations 

A set of theoretical performance calculations was carried out to establish 
the effects of injecting the lithium as a solid at liquid hydrogen tempera- 
ture instead of as a heated liquid. The tripropellant combination employed 
was LF\,(l3"> R) SI.i()b R) R), at a single chamber pressure (300 psia) 

Results are illustrated in a series of curves which are listed in Table 3 
and shown in I ig. 'i? through 73. They are briefly summarized as follows: 


1 Maximum theoretical vacuum specific impulse occurs at or near the 
stoichiometric Li mixture ratio for all conditions of percent- 
age hydrogen addition and expansion area ratio, as with liquid 
lithium (fig. '• 1 . '<8. and '»')). Optimum percentage hydrogen de- 
creases from about 30 at € = '»() to about 23 at € = 300 (Fig. 30, 

31- and 32). Maximum vacuum specific impulse at 300 psia chamber 
pressure is 312 lbf-sec Ibm ( F., I.i MR = 2.7'. . 11,,: 30 percent) at 
€ = '.0 (Fig. 30). 330 lbf-sec Ibm (F„ I.i MR = 277^ - H.,: 23 percent) 
at € = 300 (l ig. il). and 338 lbf-sec Ibm (F () l.i MR = 2.7'i. 

• - * pet* toil ) at € >00 (tig, >2). which represent decreases 

of about percent from tin* liquid lithium gaseous hydrogen 
system. Differences between shifting equilibrium and frozen 
composition vacuum specific impulse are nearly identical for the 
Sl.i I Jl. , and I.I -i tUI () combinations (Fig, 22 and 30). 

2. Maximum characteristic velocity occurs at the stoichiometric F I,i 
mixture ratio with 23- to *iO- percent added hydrogen, but is above 
mixture ratio ’t with 1">- percent hydrogen (Fig. (>()). At 300 psia 
chamber pressure and I.i mixture ratio ,,f 2.7’t, maximum c* 

(873t> ft sec ) occur* with ~>3-pMeent added hydrogen r or shifting 
equilibrium (I ig. (>! ; but with bO-percent hvdrogon for frozen 


(8)82 ft/sec. Fig. 6]). Maximum c* in both cases was with 40- 
percent hydrogen for the LLi/GII 0 system (Fig. 23 and 24). The 
difference in c* between the IJLi/GH,, and SLi/LH,, systems is 
approximately 3 percent. 

3. Variation of chamber temperature with F n , Li mixture ratio and 
with percent hydrogen (Fig. 03 and 64). at 500-psia chamber 
pressure, is the same in both systems. Decrease in temperature 
of the injected lithium and hydrogen is reflected in chamber 
temperature decreases of about 4 percent with 15-percent added 
hydrogen, 3 percent with 30-percent added hydrogen, and 8 percent 
with 40- percent addee hydrogen. 

4. Vacuum thrust coefficients of the SLi/L» 0 system differ very 
slightly from those of the LLi/GII f) system. The abrupt drop in 
C p at about 25 - percent added hydrogen is evident in both cases, 
for F^/Li mixture ratios below stoichiometric (Fig. 65. 66, and 
67 ). Again, the difference between shifting equilibrium and 
frozen composition C p decreases significantly with increasing 
percentage of hydrogen (Fig. 73). 

5. Because of the lower chamber temperatures occurring with the 
SLi/ ITI q system, the percentage of condensed mass in the chamber 
(Fig* 74) is somewhat greater than in the LLi/GH () case under 
comparable conditions (Fig. 41 ). Similarly, maximum condensation 
occurs in the nozzle at a lower expansion ratio with SLi/LH 
(Fig. 75) than with LLi/GH,, (Fig. 43). 
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ANALYTICAL EVALUATION OF LITHIUM FLUORINE- HYDROGEN COMBUSTION 


Design of a LF 0 /LLi GH } rocket engine requires information regarding: 

(l) the mechanism by which lithium burns in a fluorine and/or fluorine- 
hydrogen combustion product environment, and ( 2 ) the design parameters 
which affect the reaction completion. Accordingly, present models of 
liquid metal combustion were reviewed, and it was determined that for 
small lithium droplet sizes, the rate- limit ing mechanism inhibiting com- 
plete combustion would be the liquid lithium vaporization rote. 

With th i 3 as a basis, a modified, one-dimensional, computer program based 
on a vaporization rote-limited combustion model was used to predict char- 
acteristic exhaust velocity efficiency as a function of engine design vari- 
ables (contraction area ratio and chamber length), operating chamber pres- 
sure, amount and location of hydrogen addition, and lithium droplet size. 
These computations suggested that c* efficiency is extremely sensitive 
to drop size and chamber residence time (L*), with a lesser but significant 
dependency upon hydrogen injection location and amount. The predicted 
effect of pressure was minor, r onsiderat ion of the lithium combustion 
mechanism suggested that the dependency of c* efficiency upon drop size 
may he greater than predicted because the combustion of larger drops may 
not he vaporization rate-limited. 

LTTHTUM- FLUORINE COMBUSTION MECHANISM 

Brzustowski and Glassman (kef. l) hnv^ observed that liquid metals which 
have, (l) vapor pressures of 1 mm Hg or more tt 1000 K, ( 2 ) boiling points 
rfucb lower than the corresponding oxide boiling point, (“>) low melting 
points compared to the flame temperature, and (/j) porous oxides tend to 
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burn with vapor-phase diffusion flames. These conditions are satisfied 
by lithium burning in fluorine because; (l) the vapor pressure of lithium 
is greater than 1 mm Hg at 1000 K, (2) the fluoride boiling point is >050 F 
as opposed to the 2450 F boiling point of liquid lithium, (5) the melting 
point of liquid lithium is low (36O F), and (4) the Pil ling-Bedvorth ratio* 
of 0„7 indicates a porous fluoride. 

On the basis of reasoning similar to that presented in Re^. 1 , the follow- 
ing theoretical mechanism is suggested for combustion of lithium droplets 
injected into a Li-F^ combustion chamber: 

1. Upon injection, the liquid lithium droplet is considerably below 
its boiling point, and fluorine diffuses to the surface of the 
droplet where it reacts to form a porous layer of LiF. 

2. The droplet continues to react as the fluorine diffuses through 
the porous LiF layer. The reaction heat quickly causes the 
droplet to heat up to the melting point of LiF, at which point 
the latter becomes a molten layer coating the lithium droplet. 

1. Only small amounts of LiF have been formed on the droplet to this 
point so that the molten LiF layer is thin, particularly on the 
upstream stagnation point of the droplet during acceleration in 
the combustion gas stream. Chemical reaction continues by diffu- 
sion of fluorine through the thin liquid layer. 

As the temperature of the droplet increases, the vapor pressures 
of both the lithium and the LiF increase. Because the LiF is 
exposed to the gas phase, it gradually evaporates. 

5. The liquid lithium becomes exposed to the gaseous fluorine 

atmosphere and has significant vapor pressure, so that it begins 
to burn in a diffusion flame around the droplet. The burning is 
now much more rapid than before. 

♦The Pill ing-Bedworth ratio is the ratio of the volume of metal oxide or 
fluoride formed to the volume of metal reacted with oxygen or fluorine 
Vhen this number is lens than unity, the oxide or fluoride is porous 
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6 . The droplet becomes further heated until its temperature approaches 
the lithium boiling point, and the burning rate of the droplet 
reaches a quasi-steady state. 

It is important to note that the theoretical mechanism described above is 
applicable only to small droplets of lithium, vhich have high ratios of 
surface area to mass, so that the molten LiF layer formed is thin and 
easily evaporated in step 4. On large droplets, with lower ratios of 
surface area to mass, the molten LiF layer formed in steps 1, 2, and 3 is 
relatively thick, and it physically inhibits the heat-producing reaction 
of lithium and fluorine; this results in the layer actually increasing in 
thickness rather than evaporating as in step 4. For further clarification, 
the large and small drop combustion mechanisms are schematically illustrated 
in Fig. 76. 

It is interesting to note that if all of the heat of the Li-Fg reaction 
were used to heat a 500 F lithium droplet to its boiling point, approxi- 
mately 2 percent of its mass would be reacted with fluorine. This would 
result in a layer of LiF on the droplet with an average thickness of approx- 
imately 1.5 percent of the droplet diameter (actually, somewhat more LiF 
will be formed, because of heat losses from the droplet during the heating 
period). For initial droplet sizes of 1, 5, 10, 30, 60 , and 100 microns, 
a 1.5-percent layer of LiF has a thickness of 0.015, O.O 75 , 0.15, 0.45, 

0.9 and 1.5 microns, respectively. If a 1-micron-thick layer of LiF is 
arbitrarily selected as the maximum thickness which will allow a droplet 
to burn by the small droplet mechanism, then droplets smaller than approxi- 
mately 70 microns will burn by the mechanism illustrated in Fig. 76 -A, and 
larger droplets will burn by the mechanism of Fig. 76 -B. 
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Although maximum droplet size for reaction as a small drop is not known, 
the determining factor will he the amount of LiF formation on the droplet 
surface. Greater quantities of LiF will be formed if: 

1* A solld rather than liquid is injected. A larger time interval 
will elapse prior to steady-state combustion (time required to 
reach the metal boiling point). 

2. The atmosphere contains only dilute concentrations of fluorine. 
Under this condition, the rate-limited chemical reactions dis- 
cussed previously (specifically, steps 1, 2, and 3) will be slow, 
and greater quantities of LiF formation are necessary to achieve 
droplet heating because of increased heat losses from the droplet 
to the surroundings during the extended heating period. 

3. The surrounding atmosphere is cold. Again, additional LiF must 
be formed to compensate for heat losses. 

Thus, if the lithium is injected as a solid particle into a cold, dilute 
atmosphere of fluorine, it will have a greater tendency to burn by the 
large droplet mechanism. On the contrary, if the lithium is injected as 
a molten droplet into a hot, concentrated fluorine atmosphere, it will 
have a greater tendency to burn by the small droplet mechanism. There- 
fore, which of the two mechanisms prevails depends not only upon the size 
of the particle, but also upon the injection temperature of the lithium, 
the concentration of fluorine, and the prevailing temperature in the 
chamber. 
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MODEL OF THE LITHIUM DROPLET COMBUSTION MECHANISM 

Brzustuvski and Classman (lief, l-j) have proposed an analytical model for 
the vaporization rate-limited combustion of liquid metal droplets, which 
resembles „.e models developed for propellants such os hydrocarbons 
(Ref. 16 through 19), except that effects of thermal radiation to the 
droplets and the possibility of condensed products in the combustion zone 
around the droplet were considered. If the effects of thermal radiation 
to the droplet are not included in the combustion analysis, the model of 
Ref. I 5 becomes very similar to those of Ref. lb through 19. A schematic 
description of the single droplet burning model is given in Fig. 77 . 

The equations (Ref. 15) used in the present study to model single metal- 
droplet burning under stagnant gas conditions are: 


with 



In 


(l 1 B / 
^o 


( 2 ) 


Data used to calculate droplet burning rates by Eq. 1 and 2 are given i 
Table A, which lists all the parameters required except for the temper 
ture of the reaction zone surrounding the droplet, Tg. In some instance 
this temperature ha* been taken as equal to the metallic oxide boiling 
,'oint for metals burning in air (Ref. 15). However, for lithium burnin 
stoichiometries 1 lv in gaseous fiuorine, it has been calculated that no 
condensed products are produced and therefore T Q should not be equal 
the condensation temperature of LiF. Furthermore, the possibility of L 
condensation is decreased in the present case because the gaseous fluor 
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Vi 11 be preheated during gassif ication in the gas generator and by the 
early stages of combustion with lithium. For these reasons, the value of 
T B vas calculated separately by the following expression (Ref. 1.°.): 

1.76T +64,200 I (l + Y /2.74p*^ - l] 

T_ = ^ J , 

1.76 ( 1 + Y /2. 74) 

The theoretical derivation upon which Eq. t is based neglects the effects 
of combustion product dissociation so that it tends to give values of T 
*hicli are somewhat high, but tins, in turn, compensates for the earlier 
assumption of no radiant hc-at transfer from the reaction zone to the liquid 


If the droplet has readied its equilibrium temperature and the surrounding 
conditions do not change, the stagnation burning rate, Id, remains constant 
throughout the combustion of the droplet and is given byV.e following 
expression (itef. 15): 

- ' A L i, _ 8 \W 

dt " s " P T TT (4) 

E p 


With the data given in Table '1 and IV;. !, 2, 5, and 
constant for lithiim droplets under stagnation 
calculated a « : 


» the burning rate 
conditions was 



O.OOOVi in sec 




For cases in which hydrogen is added to the lithium-fluorine system, it 
was assumed that the total amount of hydrogen is added at a given axial 
location, and that it attains instantaneous equilibrium with the other 
species present. Addition of hydrogen results in the formation of HF, 
which also diffuses to the lithium droplet and reacts with lithium vapor 
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The higher diffusivity of the HF molecule (compared to F () ) tends to be 
offset by the formation of increased amounts of combustion products which 


musi diffuse in a countercurrent direction away from the vapor phase 
reaction zone surrounding the lithium droplets. The primary effect of 
hydrogen addition, therefore, is the dilution and the decreased tempera- 
ture of the reactant atmosphere. The addition of hydrogen affects the 
droplet combustion rate primarily as shown in Eq. 3. A 20 to 30 percent 
decrease in for lithium droplet combustion was generally found to 
occur for the amounts of hydrogen addition (l5 to >,0 percent) considered 
when all three propellants are added at the same point. 


DESCRIPTION OF THE COMPUTER PROGRAM 

The single uroplet burning rate constant, , is directly applicable to 
'“ c 0T iginal version of the Rocketdyne combustion model (Ref. 20), Only 
two minor changes to the model were required. One was the incorporation 
of Eq 3, and the other concerned alteration of the computer program to 
consider a three-propellant combination rather +hs„ ine original two- 
propellant combination. 

The combustion model considers the liquid propellants to be distributed 
into several droplet size groups, with all droplets in each group having 
the some diameter. The propellants were distributed into three size groups 
at the injection point by selection of the volume mean diameter, P.^, and 
use of a Nukiyama-Tanasawa distribution. 


The following major processes are considered in the model: (l) droplet 

vaporization under forced convection, (2) droplet drag and acceleration, 
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and (l) compressible flow dynamics of combustion gases. These three 
processes are interrelated and result in "bootstrap" combustion in the 
rocket engine combustion chamber. Droplet vaporization increases the gas 
veiucaty, ana mcmsed gas velocity increases convective heat transfer 
to, and vaporization of , the droplets. As a balancing factor, increased 
gas velocity causes greater droplet acceleration and hence shorter 
residence time in the combustion chamber. 


lor droplet vaporization under forced convection, the program computes an 
increased droplet burning rale constant by applying a convective heat 
transfer factor to the stagnation burning rate constant; 
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Droplet drag and acceleration are accounted for by use of a standard drag 
equation, together with drag coefficients specifically applicable to 
liquid droplets (Ref. 20); 
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The compressible flow dynamics of the combustion gases are handled by 
classical methods (Ref. 21 ), except that the interchange of momentum 
between the compressible combustion gases and the liquid droplets is 
considered (Ref. 20 ). 

The combustion model used considers one-dimensional, steady-state, rocket 
engine combustion. It takes no account of the ignition processes, nor 
does it consider the combustion which occurs in the poorly mixed (two- 
dimensional) regions near the injector face. This version of the pro- 
gram requires the following boundary conditions: 


1. Conditions at the start of the well— mixed, one— dimensional zone. 

a. Propeiiant flowrate 

b. Droplet size distributions and velocities 

c. Combustion gas velocity 

12. Sonic gas velocity at the throat 


The computer program starts calculation with the first boundary condition, 
it then applies the appropriate equations, in finite difference form, to 
determine the combustion processes which occur in the chamber by calculations 
at short increments (about 0.05 inch) until the throat is reached. 

The final results of the computations define the amount of lithium vapor- 
ized as a function of pertinent variables such as drop-size, chamber 
length, etc. Assuming that the amount vaporized equals the amount burned, 
characteristic exhaust velocity efficiency is express *! as a function of 
the amount ot burned and unburned propellants: 
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For the Li/F,/^ propellant combination, the fluorine and hydrogen are 
considered to be injected as gaseous propellants and are therefore 
assumed to be 100 percent burned in all cases. The only unburned pro- 
pellant which enters into the calculation of T}* by £ q . 8 is any lithium 
which does not become vaporized in passing through the combustion chamber. 

For all of the combustion model calcula*,ions , it was assumed that the 
well-mixed zone begins 2.5 inches downstream of the injector face, at 
which point 15 percent of the lithium is reacted and the velocity of all 
droplets is 50 ft/sec (Fig, 7«). The combustion gas velocity assumed at 
the 2. 5- inch location is consistent with a Fg/Li mixture ratio of 2.74, 
fluorine injected in the gaseous state, and 15 percent of the lithium 
reacted at that point. For cases in which hydrogen was assumed to be 
injected upstream of the 2.5-inch location, its presence was also included 
in the assumed gas velocity at 2.5 inches, to satisfy the mass continuity 
equation. The hydrogen injection velocity was assumed to have no effect 
on the combustion gas momentum because the specific direction in which 
the hydrogen can be injected is not fixed. At the poinc of injection, the 
hydrogen was assumed to be instantaneously mixed and completely reacted 
with the other gaseous species present. 

Propellant injection rates were estimated so that the portion which reacted 
would result in sonic gas velocity with the desired throat area. If the 
propellants burned more completely or less completely than anticipated, the 
computer program automatically adjusted the contraction ratio to the proper 
value. If the calculated contraction ratio was not sufficiently close to 
the desired value, the propellant injection rates were adjusted, and the 
entire set of calculations was repeated* 

The somewhat arbitrary assumptions of 15 percent lithium reaction and 
uniform conditions at the 2.5-inch location do not profoundly affect .‘he 
final results of the combustion model calculation. If, in fact, the 
lithium is 15 percent reacted at 1 . 25 -inches downstream of the injector 
rather than 2.5 inches, then the same overall result may be obtained by 
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merely constructing a combustion chamber 1.25 inches shorter. On the 
other hand, if the true distance is 5 inches rather than 2.5 inches, then 
a chamber 2.5 inches longer than anticipated would be needed. 

The relationship between c* efficiency and percentage of lithium burned 
is illustrated in Fig. 79, which shows that incomplete lithium vaporization 
has relatively little effect on Tj^ when hydrogen is present, compared to 
the case for zero percent hydrogen. Losses in c* are small when hydrogen is 
present because it can react with the unused fluorine. If none of the 
lithium were to vaporize and react, c* efficiency would be entirely deter- 
mined by the reaction between fluorine and hydrogen. In the Li/F,,/H 2 
combustor as presently conceived, fluorine and hydrogen are injected as 
gases; therefore, it is likely that the F^/h^ reaction will reach equilib- 
rium. Figure 80 shows the *7 C * which would be obtained if the lithium 
were partially reacted and the F 2 /H 2 reaction were complete. For hydrogen 
addition in excess of about 7 percent, the shaded region in Fig. 80, in 
which lithium reaction is incomplete, represents an range of 
approximately 20 percent. 

The characteristics of the Li/F^/ll 0 propellant combination shown in Fig. 79 
and 80 are important in interpreting the combustion model results, reported 
herein as c* efficiencies. 

COMPUTER PROGRAM RESULTS AND DISCUSSION 

Table _> lists the results of the combustion model calculations in terms 

of the following parameters: injector-to-throat chamber length, liquid 

lithium volume mean droplet diameter, chamber pressure, distance of point 

of hydrogen addition from the injector face, amount of hydrogen added 

(expressed as weight percent of the total propellant flovrate), chamber 

contraction ratio, chamber characteristic length, percentage of lithium 

reacted before reaching the throat, and c* efficiency. For all the cases 

listed in Table 5, the F„/Li injection mixture ratio was 2 . 74 . 
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The calculated efficiencies given in Table 3 for lithium volume mean 
droplet sizes of 100, l^O, and 200 microns may not be realistic, since 
these large droplets very likely do not burn by the vaporization rate- 
limited process which is characteristic of smaller sizes. Nevertheless, 
the vaporization rate- limited process was assumed to occur even for these 
large droplets, and the calculations are included in Table 3 to show the 
detrimental effects of large droplet size on combustion efficiency. The 
large sizes will probably give even lower combustion efficiencies than 
indicated, because they may actually burn by the slower mechanism. 

Figure 81 displays most of the data from Table 5 for the case of zero- 
percent hydrogen on the basis of characteristic chamber length (l*), which 
is a convenient parameter because it effectively normalizes the various 
chamber lengths and chamber contraction ratios to a single plot. The 
parametric curves represent lithium volume mean droplet sizes (with 
Nukiyama-Tanr.sawa distribution) from 20 to 200 microns. It is apparent 
that lithium volume mean droplet sizes must be less than 50 microns to 
achieve effirient combustion with a practical chamber L*. 

Figure 82 shows the effect of injecting hydrogen at the Li/F 9 injector 
face. The addition of small amounts of hydrogen actually results in 
incrersed which is a result of the F tt /\l l} reaction; the increase 

occurs even though less lithium is reacted. Hydrogen addition decreases 
the hi/F^ flame temperature and thereby retards the lithium vaporization 
process . 

Figure 83 shovi: the effect of adding hydrogen 5- inches downstream of the 
Li/F 0 injector. Comparison of Fig. M2 and 83 reveals that as much as 
percent gain in r? c# is obtained by moving the hydrogen injection point 
downstream to the 5- inch location. This gain occurs because the absence 
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ot hydrogen near the injector face results in an increased flame tempera- 
ture, and thereby accelerates the iithiuu vaporization. Downstream 
hydrogen injection also gives slower gas velocities near the Li/F 0 
injector lace, which results in a longer residence time for the lithium 
droplets and hence their more complete vaporization. 

The effect of increasing the total chamber length to 15 inches, with 
hydrogen injection at the 10-inch location, is shown in Fig 84. The 5 - 
inch increase in upstream zone length, where Li-Fg combustion 1.1 the ab- 
sence of hydrogen occurs, results in l? c# gains of approximately 2 per- 
cent. The effect of an auditional 5-inch increase in upstream zone 
length is shown in Fig 85 , in which overall chamber length is increased 
to 20 inches and hydrogen is injected at 15 inches. Comparison with 

Fig 84 shows gain in ^ of only 1 percent or less due to the 5 -i„ch 
l ncreaae. 

The effect of increasing zone length downstream of the hydrogen injec- 
tion point by 10 inches (by increasing the overall chamber length to 
30 inches while retaining the hydrogen injection point at 15 inches) is 
shown in fig 86 . which indicates that even with lithium droplet sizes 

as large as 50 microns, it is possible to obtain 99 percent combustion 
efficiency. 

Figure 87 shows the effect of adding hydrogen at the Li/F 2 injector face 
in a 20-inch-long combustion chamber. Commparison with Fig. 85 shows 
that hydrogen injection at this point results in up to 4 percent decrease 
in compared to injection at the 15-inch location. Figure 8 S shows 
that combustion efficiency is not a strong function of chamber pressure. 
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LIQUID-METAL ATOMIZATION 


The static burning rate of liquid lithium in a gaseous fluorine atmosphere 
is only about one-tenth that of conventional storable and cryogenic pro- 
pellants. To compensate for this inherently lower burning rate, the mean 
diameter of the 1 iquid— 1 i thium droplets must be approximately one— third 
that of more common propellant* for equivalent combustion in chambers of 
comparable size. 

Acomization of liquid lithium is relatively difficult because of its high 
surface tension (approximately 393 dynes/cm) compared to that of other 
liquid propellants (approximately 25 dynes/ cm). Mean drop size in liquid 
sprays is a square root function of the liquid surface tension, so that 
under normalized injection conditions the mean drop size of liquid-lithium 
droplets would be about four times that of more conventional propellants. 
Rocketdyne experience (Ref. 22) with conventional self-impinging doublet 
injectors for molten sodium indicated that a mean drop size of approximately 
180 microns was produced when extremely small (0. 0135-inch diameter) ori- 
fices were employed at high injection velocities (about 200 ft/sec). 

Because of its physical properties, liquid lithium injected under similar 
conditions would produce mean drop sizes about twice this value. It was 
therefore evident that a much more efficient atomization technique would 
be required to produce the small lithium droplet sizes required for high 
combustion efficiency. 


GAS-LIQUID ATOMIZATION 

Experience with gas-liquid atomization techniques indicated that the use 
®f high-velocity gas for aerodynamic utomization would probably provide 
the lithium drop size range required for efficient burning. With the 
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present tripropellant combination, either fluorine, hydrogen, or a combi- 
nation of the two could be considered for use as the atomizing gas. Use 
of hjdiogen alone would result in low flame temperatures and, therefore, 
slow vaporization and combustion of the lithium droplets. Addition of r. 
small amount of fluorine to the hydrogen to generate higher-temperature 
gases would result in an increase in lithium burning efficiency. However, 
the need to vary the amount of hydrogen in combination with the fluorine 
and lithium would be accompanied by significant changes in the resulting 
lithium dropsize. As a more desirable alternative, therefore, gaseous 
fluorine derived by precombustion of a small amount of hydrogen with the 
liquid fluorine was selected as the most feasible atomizing gas. A spe- 
cific advantage of this method is that the density and temperature of the 
gas can easily be varied by changing the pressure and mixture ratio in 
the gas generator. To maintain high lithium combustion temperature, the 

bulk of the hydrogen would be added downstream of the fluorine-lithium 
combustion zone. 


A 


A proven technique for gas-liquid atomization is the employment of an 
injector design in which a central showerhead gas jet is used to provide 
shear atomization of a pair of impinging liquid streams (Ref. 2. 23). 
Compared to the case of atomization by the mechanical impingement of two 
or more liquid streams, a basically different liquid breakup mechanism 
prevails when a liquid stream is exposed to high-velocity gases (Kef. 2',, 23), 
In the latter case, very small wavelength capillary waves are induced on 
the surface of the liquid body; these grow in amplitude and eventually 
sljear loose into small discrete droplets. 


The mean droplet size resulting from the breakup of a single large 
in a high-velocity gas flow into discrete, finely atomized, smaller 
lets is given by the following expression (Ref. 26): 
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Equation 9 is an expression for the mean droplet size which can be 
expected under somewhat idealized conditions, namely: (l) the original 

spherical liquid droplet, of diameter JD^, is essentially at rest before 
the atomization process occurs; (2) only a single droplet is exposed to a 
relatively large amount of gas so that there is essentially no gas velocity 
degradation during the entire atomization process due to the acceleration 
of the fine droplets; and (3) the ga3 flow is enclosed in a duct so that 
no gas momentum is lost due to interactions with the surrounding media. 

^30 symbo1 bracketed in Eq. 9 because of the ideal conditions 
associated with its derivation. 


To botimate the droplet size produced by injection of a liquid stream into 
a high-velocity gas jet, the characteristic diameter, in Eq. 9 was 

replaced by the liquid jet orifice diameter. 1^. Secondly, a correction 
was made for the loss of gas velocity due to the acceleration of t.hi fine 
particles which are produced. The assumptions required for this correc- 
tion are: (l) the droplets produced are sufficiently small so that they 

are instantaneously accelerated to the gas velocity, and (2) the atomiza- 
tion process is a steady-state process so that the liquid being atomized 
at any point is subjected to a gas velocity determined by its initial 
velocity and the amount of liquid which has been atomized and accelerated 
upstream of that point. With these assumptions, the gas velocity at any 
position can be approximated by a momentum balance: 
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( 10 ) 


Assuming that Eq. 9 and 10 apply to any given differential portion of the 
atomization process, the number of droplets, produced when an amount 

of liquid. dW^, is atomized, may be expressed by: 
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fc.quuti.on 11 can be integrated to obtain the total number of droplets pro 
duced. N, when the amount of liquid atomized, Wj , is equal to the total 
amount of liquid injected, W . The result is 
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The volume mean droplet diameter produced by the overall atomization 
process is given by: 
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If jet spreading (gas divergence) effects are unimportant. fclq.ljis an 

expression for the volume mean droplet size produced. It will be shown 

that it does, in fact, reasonably predict experimentally measured droplet 

sizes. If jet spreading effects are large, they may be estimated and 

applied by appropriate reduction of the relative gas velocity 0 

g 


It is convenient to express fclq . H in the following short form: 
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D 30 i * 1 k'q . 1 4 is the I’roplet size which would be produced if there were 

no gas velocity decrease due to momentum exchange between the gas and the 
atomized droplets. Hence, p is the increase in droplet size which results 
from the finite ratio of gas to liquid. Figure 89 is a plot of the <f> fac- 
tor as a function of the W^/W ratio. As W. /W approaches zero, <f> approaches 

unity; it does not become significantly greater than unity until the ratio 
* / • 

W^/W exceeds about 0 . 1 . 


Equation 14 does not completely express the breakup processes which occur, 
since the limiting condition at which the aerodynamic breakup of liquids 
ceases must also be considered. It has been shown (Ref. 127 ) that the aero 
dynamic mechanism of liquid breakup will prevail only under the following 
condition: 


^ > 1.0 
yifc " 


(15) 


That is to say, th*> aerodynamic breakup mechanism described by Eq . 9 
occurs and the results expressed by Eq . 14 are valid if the condition 
described by Eq.13 is first met. If the condition of Eq. 19 is marginal, 
aerodynamic breakup may start, but the gas velocity will dissipate suffici- 
ently during the process to stop it. Therefore, Eq . 14 may not be valid 
when Eq. 13 is marginally true. It should be noted that, the latter was 
obtained for application to spherical droplets, but has here been implied 
to cylindrical liquid streams. 

JET PENETRATION 

A basic tenet of the aerodynamic atomization expression (Eq. 14) is that 
the liquid stream is wholly exposed to the high-velocity gas. For a trip- 
let injector element, maximum initial exposure of the liquid surface area 




occurs when the streams are able to penetrate to the center of the gas jet. 
full penetration results in the most effective use of the available gas 
momentum. An additional desirable feature derived from full penetration 
by the liquid streams is the relative improvement in the mixing of gas 
and liquid. Without full liquid penetration, a gas-rich condition exists 
vi thin the core, resulting in degradation of combustion efficiency. Because 
of the criticality of having good penetration of the gas jet, a brief 

experimental program was conducted to define the parameters which affect 
this process. 


APPARATUS AND PROCtDURE 

A schematic drawing of a triplet injector element, with the liquid streams 
penetrating completely to the center of the gas jet, is shown in Fig. 90. 
Ihis condition i:. described by: 


ligure 91 is a schematic of the apparatus used to study jet penetration at 
ambient pressures. Photographs were taken of flow tests with injection 
conditions ranging from nonpenetration to complete penetration. The gas 
used was nitrogen, the flowrate of which was measured by a choked orifice. 
Water and Cerrosafe*were the liquids used. 

ITie tests covered a relatively large range ( 0.125 to 1. 5 0 inches) of gas 
jet diameter. Twc liquid stream diameters, 0.040 and 0.060 inch, were 
used. The angle cf the liquid streams with respect to the injector face 
was maintained constant to limit the number of test injectors. It is 
desirable to minimize this angle since smaller angles give better jet pene- 
t rat ion; 45 degrees was chisen because it is close to the smallest angle 
which can be used without causing backsplashing of liquid onto the injector 


*Cerroaafe propertiei are given in Table 7, 




face. rhe nitrogen density was 0.065 lb/ft 3 (ambient pressure and tempera- 
ture) for the majority of the tests; however, four test series were con- 
ducted at elevated pressures with nitrogen density as high as 1.21 lb/fi 3 . 
The higher densities corresponu roughly to the fluorine gas densities used 
in the Li-F 2 injector. Gas velocities were varied from 380 to 1100 ft/sec 
(Mach 1.0) with two liquid densities: 62.3 lb/ft 3 (water) and 590 lb/ft 3 

(Gerrosafe metal alley). The liquid Cerrosafe tests were carried out 
because these data would be of interest in connection with droplet size 
distribution tests. Water, on the other hand, has a density more compar- 
able to that of lithium and is more easily handled than Cerrosafe; it was 
therefore used for the majority of the tests. The critical liquid veloci- 
ties (i.e., those required for penetration of the liquid stream to the 
center of the gas jet) were obtained by suitably adjusting the liquidAP 
in each individual test. 


Results and Discussion 


Test data for the jet penetration experiments are given in Table 6, and 
typical photographs are shown in Fig. 92. Figure 9&A shows liquid-stream 
impingement in the absence of gas flow through the showerhead orifice of 
the triplet element; Fig.92-B shows a condition in which the liquid streams 
do not have adequate velocity for gas jet penetration, and Fig.92-C illu- 
strates the condition in which the liquid streams have sufficient velo- 
city to penetrate to the center of the gas jet. 

The data obtained from the cold-flow jet penetration studies indicate 
that the physical mechanism primarily responsible for the limitation of 
liquid-stream penetration into a flowing gas jet is the aerodynamically 
caused breakup of the penetrating stream into small fragments. Such dis- 
integration of the liquid streams is apparent from the test photographs 
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* ,,ch US IiK ' <)2 ~ B and 92_C - Breakup of the liquid stream results in drop- 
lets which are much smaller than the original diameter of the cylindrical 
stream. iheso small droplets are very rapidly accelerated by the gas jet, 
and their velocities therefore become almost identical to the gas velocity 
very soon after they are formed. The rapid acceleration of the small 
droplets effectively prevents further liquid penetration. 


Since breakup of the liquid stream is the predominant mechanism limiting 
jet penetration, correlation of the experimental data vas based on the 
flighl time" of the liquid (defined as the time during which it is exposed 
to the high-velocity gas flow), and the calculated time required for the 
jet to break up due to the aerodynamic forces. The flight time of the 
liquid is given by the following equation: 

X 

l f “ Vj cos Q (17) 


lor the breakup time of the cylindrical liquid stream, 
sion, derived and experimentally verified for applicat 

of relatively large spherical liquid droplets into sma 
was used: 


the following expres— 
ion to the breakup 
11 droplets (Ref. 26), 



( 18 ) 


The value of t f) given by Kq . 18 is the time between the first exposure of 
the liquid to the gas jet and the start of actual liquid breakup. IMs 
"dead time," or preparation time, is a result of the necessity to first 
generate disturbances (capillary waves) on the surface of the original 
liquid droplets. When the disturbances are critically large, liquid 
breakup begins. The lime required for cc-pletion of the liquid breakup 
is approximately equal n the initial preparation period. The total time 
from initial exposure of tie liquid to he gas to the completion cf the 
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breakup process is thereiore about twice the value of t^ given by Eq.18 
Since this equation was originally developed to be applicable to liquid 
droplets rather than to cylindrical liquid streams, some modifications 
to it might l,e required to account for the differences between the two. 
However, because breakup time was used only as an empirical correlating 
parameter, no consideration was given to the geometrical difference between 
spheres and cylinders. 


Equation 17 gives the exposed flight time of the liquid stream, which is 
equal to its total breakup time (or slightly more, due to continued pene- 
tration of the gas jet by the atomized liquid even after the actual cylin- 
drical liquid stream has been completely broken up). For any particular 
situation, then, it should be expected that t f will be slightly more than 2t 

Values of t f and t fe are given in Table 6 and are plotted in Fig. 93 where, 
in view of the approximate nature of this type of data, they are well cor- 
related by the equation 

t f = 2.5 t b (19) 

The coefficient 2.5 is consistent with the expectation that it should be 
slightly larger than 2. 


Combination of Eq. 17 through 19 gives: 
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Equation 20 essentially expresses the experimental results of the jet 
penetration studies. 
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In another series of experiments, schlieren photographs vere taken of two 
opposing 0.10-inch-diameter helium jets penetrating into a gaseous nitrogen 
jet exhausting from a 2.0-inch-diameter tube, to aid in development of 
design criteria for the hydrogen injector. Details are given in a subse- 
quent section of this report. The data from these experiments were used 
to evaluate a coefficient for an equation of the same form as Eq. 20: 


X 

JL_ 

D„ 


3.2 





• f 

cos $ 


( 21 ) 


where cos $ was 1.0. The numerical coefficient. 3.2, is close to the coef- 
ficient (2.5) found for liquid-stream penetration into a gaseous jet. 

This \»as not unexpected since the aerodynamic forces greatly overwhelm the 
surface tension and viscous properties of the liquids so that inertial 
resistance remains as the primary factor to resist breakup. 


EXPERIMENTAL STUDY OF LIQUID-METAL ATOMIZATION 

Following determination of the conditions required for efficient penetra- 
tion of a gas jet by a liqu.d stream, the more basic problem of attaining 
a high degree of liquid lithium atomization was investigated. The objec- 
tive of this study was to determine the validity of the atomization equa- 
tion of Wolfe and Andersen (Ref. 26) for specific practical applicability 
to the design of the gaseous fluorine/liquid lithium injector. 


The basic Wo If e-Andersen atomization equation has er been experimental ]v 
verified for liquids of very high surface tension such as liquid l.thium. 

In particular, experimental verification was required of the modified 
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Wolfe-Andersen equation derived herein for the case in which only a limited 
amount of gas is available for the liquid atomization. Experiments were 
therefore carried out to verify these atomization expressions, using a 

* suitable simulant for liquid lithium. 

f 

S' 

f; 

r Apparatus and Procedure 

l 

Cerrosafe* was used as lithium simulant for the droplet-size measurements. 
This had two advantages: (l) after atomization, the metal droplets which 

were formed were easily frozen and collected, and (2) Cerrosafe has high 
surface tension (400 dynes/cm, measured in this laboratory) approximately 
equal to that of lithium. No data are available for the viscosity of 
Cerrosafe. However, the viscosities of most liquid metals at their melt- 
ing points are in the approximate range of 1 to 3 centipoise. It was 
assumed for computational purposes that Cerrosafe has a viscosity of 
1 centipoise. Volume change data are included in Table 7 to show that no 
significant error is introduced by the assumption that solid particles are 
the same size as the liquid droplets from which they are formed. 

figure 94 is a schematic representation of the apparatus used for the 
liquid-metal atomization droplet size measurement studies. Cerrosafe was 
contained in a tank immersed in water maintained at approximately 200 F 
to ensure that the alloy was molten. The outlet of the Cerrosafe tank was 
connected through a valve to the triplet injector manifold supplying the 
liquid doublet orifices. Gaseous nitrogen was supplied to the element 
through an orificed line as indicated. 


♦Cerrosafe properties are given in Table 7 


The spray produced by the triplet element vas directed into a collection 
vessel, the bottom of which was covered with water. Cerrosafe droplets 
were exposed to ambient-temperature gas during their descent and were 
solidified when they entered the collection vessel, although some of the 
larger droplets probably did not freeze until actually contacting the water 
layer. After each test the Cerrosafe particles in the collection vessel 
were sampled for photomicrography. A mechanical shutter was provided 
between the triplet injector element and the collection vessel to divert 
the spray during the transient portions of the tests. Cerrosafe flowrates 
were calculated from the injection pressure drop (C D = 0.7), and nitrogen 
flowrates were calculated from the pressure drop across a choked orifice. 

Two triplet element injectors were used for the droplet size tests. Both 
had 0.060-inch-diameter liquid doublet orifices inclined 45 degrees to the 
injector face. The showerhead gas orifice diameter was 1.0 inch for one 
injector and 1.50 inch for the other. Preliminary experiments had shown 
0. 060-inch liquid orifices to be approximately the lowest usable size; 
intermittent plugging occurred with smaller diameters due to formation of 
oxides and other slag. The relatively large gas orifices were necessary 
to obtain sufficient total gas momentum without choking the flow of the 
ambient-pressure nitrogen at the orifice exit. Figure 95 is a photograph 
of a typical Cerrosafe atomization test. 

The solidifed Cerrosafe droplets sampled from the collection vessel were 
dried and a portion of the sample was photomicrographed at magnifica- 
tions ranging from 35 to 200 X, depending on the relative droplet sizes 
in the sample. A typical photomicrograph is shown in Fig. 96. It was 
found that collected particles smaller than 40 to 50 microns resembled 
chrome-plated spheres, whereas larger droplets were somewhat distorted 
from the spherical shape. 
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To measure droplet sizes from the photomicrographs, a transparent sheet was 
placed across the photograph in a random orientation. The transparent 
sheet had a very thin line scribed upon it, and the diameters of all drop- 
lets images vhich touched the scribed line were measured and recorded. 

The process was repeated until a minimum of 200 individual droplet images 
were measured. A weighting factor, equivalent to the inverse of the diam- 
eter, is applied to each droplet counted when mean diameters are calculated 
from such image measurements (Ref. 28 ). The defining equation for mean 
droplet size, in terms of the measured droplet images, is: 
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A computer program was written to sort the droplet measurements according 
to size and to calculate the mean droplet sizes for values of q and p from 
0 to 6 inclusive. The mean droplet size of greatest interest is the volume 
mean droplet size, D^ Q , obtained when q = 3 and p = 0: 
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Another factor calculated for the droplet sprays vas G(d), the total volume 
fraction of the spray contained in droplets with diameter less than D: 
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The volume distribution of the spray is obtained by graphical differenti- 
ation of G(D). 


The weighting factor, l/D^, has been retained separately in Eq.22 through 2-'i 
to emphasize its presence. 


It is apparent from l.q . 21 and 24 that the larger-diameter image measure- 
ments are the most important in the determination of and G(D). The 
measurement technique employed tends to increase the frequency of measure 
ments at the large end of the droplet size spectrum, and although the 
higher frequency of measurements at the large end is accounted for by the 
weighting factor, this in turn improves the accuracy of the summations of 
the larger droplets at the expense of the smaller. Since the large end 
of the size spectrum is the most significant, a definite advantage is 
gained by measuring only droplets which touch a randomly positioned line. 

rhe following expression (Ref. 29) vas used to estimate the accuracy of 
^qp' <i®fin®d by Eq. 22, as a function of measured particle sizes. 
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Equation -5 gives the estimated maximum error in D at the 93-percent 


confidence level which occurs because a finite number, M, of droplets are 
counted; it does not include errors from other sources (sampling technique, 


experimental errors, etc.). ihe measured D_ 0 values are estimated (Eq. 2~)} 


to have errors ranging from (> to lf> percent. 


Experimental Results and Discussion 

Kleven experiments were carried out to obtain drops ize measurements for 
comparison with the sizes predicted by the atomization equations described 
above. the salient parameters varied during this study were the Weber to 
Reynolds number ratio (Eq. 13) and the quantity of gas available per pound 
of .iquid to be atomized. Ihe effect of the latter parameter was evalu- 
ated in terms of the parameter ^(Eq. 14). In geneial. the data obtained 
showed that the measured dropsizes agreed reasonably well with the calcu- 
lated sizes and that a critical We /N /Ri ratio exists below which the 
derived equations were not applicable. 

Ihe experimental conditions and results are presented in Table 8 . The 
parameter P included in this table represents the calculated percentage 
of liquid-stream penetration into the gas jet (P = 100 represents 
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penetration to the center of the K as jet). Values in excess of 100 percent 
mean that the liquid streams possess more than the necessary velocity to 
penetrate to the center and therefore collide. 

The data obtained show that triplet ( 1 iqui d-gas-1 iquid) element atomiza- 
tion is a somewhat complex phenomenon. Typical data were plotted to 
illustrate the overall droplet mass distribution as a function of dropsize. 
Figure 97 shows the distributions obtained when 0 v .is > 2.0 (test 52), 

< 1.9 (test ’il), and ~ 1.9 (tests '.2 and h8 ) . With the largest values of 
<P, fairly broad, flat-shaped curves typically resulted, whereas with the 
smallest values of 0 much sharper curves were produced. In the intermedi- 
ate 0 range the distribution curves have a double-peak shape. These 
results are consistent with those which would be predicted since larger 
values of 0 are caused by reduction in the quantity of gas available for 
atomization. Because of this, a greater gas velocity gradient downstream 
of the injection point occurs. The broader range of relative gas-to- 
1 iquid velocity thus created results in a larger dropsize range. 

Three tests (V2, 'tl, and Vi) were conducted in which the We \jRo ratio 
was less than, or near, unity. The distribution obtained in run ',2 is 
shown in Fig. 97. It can be seen that a very broad dropsize range was 
produced under these conditions, and as shown in Table 8 , the experimental 
mean drops! zes were much larger than calculated. In view of these obser- 
vations, it seems apparent that a different mechanism of atomization 
occurred. 


The calculated value of the penetration parameter, P, was over 80 
percent for 10 of the 11 experiments. Within this range, the penetra- 
tion parameter should have negligible effects on the dropsize character- 
istics measured, However, one experiment (test was conducted in which 
the penetration parameter was reduced to percent, but the values of 
calculated and measured mean dropsizes remained similar. The rather 
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small influence of the decree of penetration on dropsize might be expected, 
since a -.1 e..ange in the penetration parameter only causes an 8— nercent 
change in the value of 0. This small change occurs because, with 50- 
percent penetration, 75-percent of the gas flowrate is involved in the 
liquid atomization. These results indicate that penetration will probably 
influence propellant distribution more significantly than it does the de- 
gree of atomization. 


In summary, the experimental data show that the gas-liquid atomization 
equation developed (Eq. 14) is applicable to the design of the gaseous 
fluorine-liquid lithium injector. Further, application of this equation 
is limited to cases in which the We ^Ke ratio is greater than unity. 
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THRUST CHAMBER CONCEPT AND 
DESIGN OF COMPONENTS 

The tripropellant thrust cheimber concept developed in this program was 
based on the premise that lithium/fluorine combustion should be com- 
pleted prior to hydrogen addition and subsequent achievement of thermal 
equilibrium. Thus, in essence, a two-stage system was employed, in the 
first of which hot LiF was formed (theoretical temperature = 9800 R, 
molecular weight = 26), followed by hydrogen addition resulting in 
reduced temperature (3500 to 5400 R) and molecular weight (5 to 9). 

Component designs were to be flexible to the extent that variations in 
upstream and downstream chamber volume could be made so that the com- 
bustion and mixing volumes required for high efficiency could be defined 
experimentally. Potentially, in the limit, the design would allow the 
determination of the basic necessity for a two-stage system as opposed 
to injection of all propellants in the same vicinity. 

4 

The basis upon which the Li/Fg injector was to be designed depended 
upon the results of calculations based on the analytical combustion ‘ 

model and on the simulated lithium atomization and propellant distri- 
bution criteria obtained in the cold-flow atomization study. In summary, j 

the combustion model results indicated that a maximum volume mean 1 

lithium drop size of about 20 microns was required for high performance. 

The cold-flow data showed that to achieve this drop size, gas atomi- 
zation of the liquid lithium would be required. Accordingly, the j 

lithium injector concept which was developed consisted of the injection j 

of lithium streams into high-velocity gas jets. These gas jets were to 
be obtained by use of a fluorine-rich gas generator upstream of the 

lithium injector. ' 
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A conceptual illustration of the thrust chamber assembly is shown in 
Fig. 98. A schematic scale drawing of the actual assembly is presented 
in Fig. 99. 

Nominal design operating conditions for the system were: 

Mixture ratio, F Q /Li = 2.74 

Hydrogen added, percent = 15 to 40 of total flowrate 
Chamber pressure, psia = 500 
Sea level thrust, lbf = 2000 






? 
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Design criteria and descriptions of the various components are 
discussed below. 


GAS GENERATOR DESIGN 

The gas generator concept arose from the need to supply hot (500 to 
1000 F), pressurized (400 to 700 psi) fluorine to the Li /F r) injector 
to obtain a high degree of lithium atomization. It was recognized that 
a convenient way of obtaining such hot gas was to react liquid fluorine 
with a very small amount of hydrogen, at mixture ratios which would pro- 
duce the desired gas temperatures. Figure 100 shows the variation of 

G ^2 combust i°n gas temperature in the mixture ratio range 200 to 500 
at chamber pressures of 400 to 1000 psia. Corresponding shifting- 
equilibrium c* values are shown in Fig. 101. Within the mixture ratio 
range of interest (500 to 500), the small amount of HF which would be 
present in the gas generator output would not significantly alter its 
properties from those of pure fluorine. The gas generator assembly is 


1 

I 





60 


I 




shown schematically in Fin;. 102. H consisted of separate hydrogen and 
fluorine injectors, two identical chamber sections, and a nozzle plate. 

Design details and the function of each of the components are discussed 
below. 


Hydrogen Injector 

Since the hydrogen flowrate into the gas generator was necessarily very 
small (0.005 to 0.009 lb/sec), use of hydrogen injection momentum to 
atomize the liquid fluorine, as in a triplet element, was not practical. 

The hydrogen was injected through nine 0.0135-inch, showerhead orifices 
drilled on a 2-inch basic diameter. Design pressure drop through these 
small orifices (square entrances) was such that flow of ambient-temperature 
hydrogen through them was sonic. The injector was fabricated from OFHC 
copper (Fig. 103) . 


Fluorine Injector 

Liquid fluorine was injected through 36 circumferential, self-impinging 
doublet elements arranged around the inner diameter of the stainless-steel 
ii.jector (Fig. I0>i) and forming radially directed sprays which intersected 
the hydrogen jets. The orifice diameter was 0.02b inch and the impingement 
angle was 60 degrees; design pressure drop at the nominal fluorine flowrate 
of 2./>i lb sec was 250 psi, and injection velocity was 110 ft/sec. Cold 
flow of the hydrogen-fluorine injector combination with water and nitrogen 
as propellant simulants, at relative momenta approximating those in a hot 

firing, is shown in a photograph (Fig. 105) which indicates a good degree 
of liquid atomization. 



Combustion Chamber 






The function of the combustion chamber was to permit interaction of the 
fluorine and hydrogen and to provide uniform gas properties at the 
chamber exit. Calculations made with the one-dimensional, vaporization 
rate-limited combustion model previously discussed showed that combusion 
would be complete within a 5-inch chamber length for 80-micron fluorine 
droplets, and w .thin 10 inches for 120-micron droplets. These droplet 
sizes are of the order of magnitude produced by the 0.026-inch impinging 
fluorine streams. Hence, a chamber consisting of two identical 5-inch 
sections was designed to permit use of either 5 or 10-inch chamber 
lengths. Since gas temperatures within the combustion chamber are rela- 
tively low, 1018 carbon steel was selected as the material of construction. 
To promote gas mixing within the chamber, perforated stainless-steel 
mixing plates were welded at the center of each section. Experimental 
temperature data showed that they effectively increased the degree of 
combustion gas uniformity. 

Nozzle Plate 

The stainless-steel nozzle plate (Fig. 106) was both a model of the in- 
jector into which the gas generator output was to be directed and an 
exit nozzle for control of chamber pressure during gas generator evalua- 
tion tests. For firings at 400 to 500-psia chamber pressures, the 
nozzle plate contained thirteen 0.185-inch orifices. A plate with seven 
0.213 orifices was used for tests at 700 psia. Figure 107 shows the gas 
generator components in order of assembly, with only one of the chamber 
sections* 
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t J T in um/fluorine injector 

The Li/jjg injector was the most important component of the thrust 
chamber because of its profound effect upon combustion efficiency. Since 
its design posed several unique problems, two distinct designs were 
coopdeted: (l) a triplet patter:- (two lithium streams impinging into a 
showerhead fluorine gas jet) and (2) a doublet pattern (one lithium stream 
injected into a showerhead gas jet). The design requirements, i.e., the 
necessity for gas atomization and for adequate liquid penetration into 
the gas jet, were the same in both cases. 


The major differences between the two injectors were the method and mate- 
rial of construction. The triplet, which was considered functionally 
superior in terms of providing a more uniform propellant distribution, 
was, by necessity, fairly complex and required machinable metal construc- 
tion. The durability of this design was unknown and it was limited to 
temperature levels which the selected metal (stainless steel) was capable 
of withstanding. The doublet design was based upon the assumption that 
the heat flux rates encountered would be too severe for common metals, 
and therefore a simplified design which permitted the use of a more 
thermally compatible material would be necessary. For this reason, the 
doublet design employed graphite as construction material and it was 
used for the major portion of the experimental program. 


Triplet Injection Pattern 

.Liquid lithium and gaseous fluorine were injected through triplet 
elements with the outer liquid lithium streams intersecting at a 
90 degree angle over a central fluorine gas orifice. As a result of 
data obtained during the lithium flow studies (described in detail in the 
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Liquid Lithium Flow Characterization section), 0.043-inch-diameter 
lithium orifices were selected for a seven-element injector; the design 
pressure drop at nominal flowrate (l.OO lb/sec) was 275 psi (Cp = 0.8). 


Results obtained in the metal atomization studies reported in a pre- 
ceding section indicated that the triplet element should be sized so 
that the outer liquid metal streams would have sufficient momentum to pene 
trate the central gas jet completely (i.e., to its center). Application 
of this criterion for optimum atomization showed that complete penetration 
is difficult to achieve at high gas velocities; a compromise is required 
between high gas velocities, which produce finely atomized metal droplets, 
and low gas velocities, which ensure complete penetration and hence effi- 
cient aerodynamic contact between gas and liquid. Thus, high gas velocity 
results in the production of finely atomized metal droplets for that 
portion of the liquid which penetrates the gas jet. The degree of pene- 
tration may be so low, however, that a substantial fraction of the liquid 
does not contact the gas and is therefore very poorly atomized. 
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Equation 26 is a continuity equation which expresses the F^Li injection J 

mixture ratio as the ratio of the weight flowrates through the single 1 

fluorine orifice and the two lithium orifices of the triplet injector 
element. Equation 20 was developed in the cold flow portion of this 
program; it describes the penetration of the liquid lithium streams into 
the central showerhead fluorine jet. Equation 14 is a modified form of 
the Volf e-Anderson equation (Ref. 26) which was verified in the cold-flow 
portion of this program; it expresses the lithium droplet size as a function 
of the injection conditions and propellant properties. Factors which must *1 

be considered in the use of these equations were discussed in the section i 

of this report which deals with liquid metal atomization. 

Results obtained with Eq.14, 20, and 26 are presented in Fig. 108A and 
108B, , which show the predicted extent of lithium stream penetration 
into the fluorine jet and the predicted volume mean droplet size of the 
atomized lithium, respectively, (it should be noted that 50-percent pene- | 

tration in Fig. 108 represents penetration to the gas jet centerline.) 

Lithium is injected at 0.0714 lb/sec/orifice through 0.043-inch orifices 
inclined at an angle of 45 degrees to the face. The effects of variable 
fluorine injection temperature and orifice diameter and of variable chamber 
pressure are indicated. It is apparent that larger fluorine orifices allow 
better jet penetration along with increasing droplet size; the increased 
degree of penetration is beneficial whereas the increased drop size is 
detrimental. Conversely, smaller fluorine orifice diameters produce 
smaller droplet sizes but increasingly poor jet penetration. The 
selected size for the design chamber pressure (500 psia) was 0.272 inch. 

This size represents the approximate optimum design condition for liquid 
penetration to the gas jet centerline and minimized drop size for the 
nominal gas generator design operating temperature of 1000 R. 


Triplet Material Selection. The most satisfactory materials for hot 
gaseous fluorine service are nickel, monel, copper, aluminum, and, to a 
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lesser degree, the 300-series austenitic stainless steels (Ref. 30) . 

For liquid lithium service, available data (Ref. 31 and 32) indicate 
that the best materials are pure iron, stainless and carbon steels, and 
certain refractory metals (molybdenum, tantalum, columbium); nickel and 
nickel alloys are probably satisfactory for temperatures under 600 F, 
while copper, aluminum, and their alloys have poor resistance to attack 
by molten lithium. 

Selection of the injector material involved a number of tradeoffs, not 
only because of conflicting compatibility properties but also on account 
of relative cost, availability, ease of machining, and high-temperature 
strength. Consideration of these factors led to the selection of stainless 
steel (321) us the basic injector material, with nickel as the best alter- 
nate. Materials bordering on state-of-the-art development were not con- 
sidered as candidates because of scarcity of compatibility data and generally 
excessive cost and delivery lead time. The relatively low thermal conduc- 
tivity of stainless steel was a recognized disadvantage in this application 
with regard to both preliminary heatup time and firing duration capability. 
For the latter reason, an injector design was developed which allows pro- 
tection of most of the injector face from the combustion gases during 
firing (described below). 


Triplet Injector Heating . To heat the injector to 500 F and keep it at 
this temperature prior to firing, heater rods (3-kilowatt total capacity) 
arranged in a radial pattern around the injector center were provided. 
Although some injector heating was obtained from the hot gaseous fluorine 
as well as from hot helium purge gas prior to the flow of lithium, these 
were auxiliary heat sources, not adequate in themselves to satisfy total 
injec'tor heating requirements. 

Injector Fa ce Thermal Protection . The extremely hot (^9800 R) 
combustion gases in the recirculation zone imposed a requirement for 
thermal protection of the injector face. Such protection was provided 
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by inclusion of a 0.75-inch layer of carbon cloth/phenolic resin abla- 
tive over the entire face of the injector except for minimum areas 
around the propellant orifices. The carbon cloth layers in the ablative 
were oriented parallel to the injector face. This ablative protector 
could be replaced after each firing series if necessary. 

Photographs of the triplet injector and its components are presented as 
follows: Pig. 109 shows the inlets to the seven fluorine gas orifices; 

Pig. HO is a face view of the triplet injector with the protective 
ablative face plate installed; Pig. Ill and 112 are face and back views 
of the ablative plate, respectively; and Fig. 113 shows the triplet 
injector with the face plate removed* 



Doublet Injection Pattern 

The lithium/fluorine injector is subjected to extremely severe chemical / 

and thermal environments. On the inlet side, it is exposed to fluorine * 

gas at temperatures up to 800 F and pressures up to 700 psi, as well as 

to liquid lithium at 500 P and 500 to 800 psi; at its face, combustion 

chamber temperature may be as high as 9800 R for firing durations up to 

8 seconds. No comon metal can withstand these conditions without J 

protective devices such as thermal resistive coatings or ablative shields, j 

and the adequacy of such protectors is questionable. Better results i 

might be obtained with refractory metals, but not enough is known of 

their compatibilities to justify their high cost and difficulty of 

fabrication* 


Of the nonmetal lie materials, possible candidates are graphite, ceramic 
oxides, or metal carbides. Graphite was selected as the material for 
the alternate injector design, for the following reasons: 
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1. Experience with graphite in combustion chambers and nozzles 
with fluorine and other fluorinated oxidizers has been generally 
satisfactory, although experimental conditions to date have 

not been as severe as those to which the lithium injector might 
be subjected. 

2. Various properties of graphite are favorable for this appli- 
cation (high thermal conductivity, good resistance to thermal 
shock, chemical inertness, ease of machining). 

3. The common ceramics and carbides have one or more disadvantages 
which, at present, rate them below graphite for use as injector 
material (lack of appropriate compatibility data, brittleness, 
variable porosity, difficulty of machining, high cost). 


Use of triplet elements in the first Li/Fg injector required fairly 
elaborate internal manifolding which could be eliminated by conversion 
to LLi GF^ doublets. Calculations of the degree of penetration of the 
fluorine gas jet by the liquid lithium stream and of the mean droplet 
sizes of the atomized lithium indicated that acceptable levels of both 
parameters could be obtained with a doublet injection element. Thus, 
for an 8-element injector with fluorine orifice diameter of 0.297 inch 
and lithium orifice diameter of 0.055 inch, under the injection con- 
ditions for firings at 500-psia chamber pressure, the calculated degree 
of gas jet penetration was approximately 85 percent, and calculated mean 
lithium droplet diameter was approximately 10 microns. 

Fluorine from the gas generator was injected through showerhead orifices. 
Liquid lithium was injected through individual tubes to each element 
through orifices inclined at 'i5 degrees to the injector face. The tubes 
were 321 stainless steel and were fed from a manifold outside the injec- 
tor body. This assembly is shown in Fig. 114. Figure il5 is a face view 
of the injector. 
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The external lithium manifold was close-coupled to the lithium/fluorine 
injector and was heated by Themwire heating cable wrapped around it. The 
lithium inlet tubes were not electrically heated within the injector body 
but were heated instead by passage of hot helium purge gas through them 
immediately before the lithium flow. Calculations showed that the inner 
surfaces of the tubes could be heated to about 500 F by flowing helium gas 
at 800 F through them for 4 seconds at normal purge gas pressure. This 
was confirmed experimentally by embedding a tube in a block of graphite, 
heating it by a 4-second flow of helium at 800 F, and then flowing liquid 
lithium through it without difficulty. Preheating of the tubes by con- 
duction from the heated manifold provides added assurance of temperature 
attainment . 


COMBUSTION CHAMBER 
Chamber Geometry 

Design criteria for the F^/hi combustion chamber were established from 
results of calculations based on the Lambiris-Combs one-dimensional, 
steady-state combustion model (Ref. 20) . Assembled results of these 
calculations are shown in Fig. 81, which present T)^ of the F /hi 
system as a function of L* and mean lithium droplet size. 

At 500-psia chamber pressure, lithium combustion occurs primarily in the 
first chamber section (item 7, Fig. 99). The actual residence time of 
the propellants in the first chamber section is nearly constant and inde- 
pendent of the amount of hydrogen injected or the threat size, provided 
the chamber pressure and the flowrates of fluorine and lithium are constant. 
Consequently, computation of effective L* for the first chamber section was 
based on the throat area which would give 500-psia chamber pressure when 
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fired at the nominal F^/hi ratio (2.7 4) and total propellant weight flowrate 
(3.74 lb/sec), with zero hydrogen injection. As described below, the throat 
area required was 1.75 sq in.; with a 5-inch chamber diameter, this corres- 
ponds to a contraction ratio of approximately 11 and an effective L* of 180 
inches . 

It is apparent from Fig. 81 that increasing the effective L* of the first 
chamber section beyond 18(. inches results in only minor gains in calcu- 
lated combustion efficiency, even with lithium droplet sizes as large 
as 50 microns, whereas decreasing L* below 180 inches has an increasingly 
adverse effect. Hence, the 180-inch value, corresponding to a chamber 
length of 15 inches, was selected for the first chamber section. The 
5-inch, second chamber section was designed to allow changes to be made 
in the overall L* during the experimental program by permitting use of 

5-, 15-, or 20-inch combustion chamber lengths; the latter is shown in 
Fig. 99 . 


Chamber Materials 

The F^/Li combustion chamber upstream of the hydrogen injector is sub- 
jected to extremely high- temperature combustion products. Heat transfer 
calculations with the Bartz equation (Ref. 33) indicated that uncooled 
topper chamber sections would be satisfactory for 10-second firings. 
However, Rocketdyne experience with fluorine and other fluorinated oxi- 
dizers has shown that heat flux near the injector face is frequently 
much higher than values calculated from the Bartz equation and may, in 
fact, approach that at the nozzle throat. It has also been found that 
under these conditions and in +hc- absence of water, ATJ graphite gives 
excellent service, nence, this material was chosen for the combustion 
chamber sections in the form of 1 . 5-inch- thick cylindrical liners. 
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Compressible, porous graphite washers were used to accommodate stresses 
associated with thermal expansion of the graphite liners in the axial 
direction and asbestos sheet in the radial direction. Steel pressure 
shells were used for both chamber sections. 

HYDROGEN INJECTOR 

The hydrogen injector should introduce gaseous (ambient) hydrogen into 
the flowing F^/hi combustion gases in such a manner that efficient mix- 
ing occurs in the chamber section downstream of the injector. This 
implies a requirement for maximum unformity of hydrogen distribution 
across the chamber, which may be satisfied by use of tranversely directed 
orifices arranged so that each one directly opposes another across the 
chamber. The impingement point of each pair of jets then becomes, in 
effect, a hydrogen injection point from which diffusion and turbulent 
mixing may begin. The principal advantage of this simple concept is its 
capability of providing uniform mixing without use of projections into 
the high- temperature gas stream. 

Jet Penetration 

The hydrogen injection concept was based on the premise that each of the 
gaseous hydrogen jets would be able to penetrate the F^/Li combustion 
gases and impinge upon its opposing jet. To verify the penetration 
capability of gas jets, cold-flow schlieren photographs were taken to 
observe the penetration characteristics of high-velocity ( 650 , 1300, and 
2600 ft sec), 0. 1 -inch-diameter , helium jets into a low-velocity ( 1 5*5 and 
110 ft sec), 2.0- inch-diameter gaseous nitrogen stream. 


The schlieren experiments were used as a design basis for the hydrogen 
injector. Since ii was uncertain whether jet breakup or jet turning is 
the limiting factor for gas jet penetration, methods based on both these 
phenomena were used to scale the He/ GN 0 cold-flow results to the H / 
(Li-F 2 ) hot f ring situation. If jet turning were the limiting factor, 
the cold-flow results indicated that all hydrogen jets would penetrate 
the Li/Fg combustion gas stream and form the patterns indicated in Fig. 
116. This corresponds to a condition in which the hydrogen mass injection 
rate is 15 percent or more of the total propellant flowrate. On the other 
hand, if jet breakup were the limiting factor, the hydrogen injection rate 
must be at least 22 percent of the total propellant weight flowrate to 
produce the same results. The percentage of hydrogen injection required 
for complete penetration, therefore, is between 15 and 22 percent. 


The hydrogen injector is expected to be an efficient mixing device at the 
higher hydrogen injection rates (22 to 40 percent); at the lower injection 
rates (15 to 22 percent) the mixing efficiency might possibly drop because 
the opposing hydrogen jets may not be ideally cylindrical but rather 
elongated and distorted when they impinge. Acceptance of possibly less 
satisfactory hydrogen jet penetration at the lower injection rates was 
the preferred choice over the following two alternatives: 

1. Use of higher hydrogen injection velocities. The resultant flows 
at the higher injection rates would raise the requisite hydrogen 
supply pressure to prohibitive levels and, in addition, would 
produce underexpanded, "bushy" hydrogen jets. 

2. Use of fewer hydrogen orifices, with increased diameter. This 
would provide better jet penetration capability, as indicated by 
Eq. 21. To maintain the same hydrogen injection velocity, a 
smaller number of orifices would be required. This would result 
in fewer hydrogen jet impingement points, which might degrade 
mixing efficiency. 
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Hydrogen Injector Design 

The hydrogen injector design consisted of a ring-shaped copper segment 
with the same inside diameter as the combustion chamber and containing 
64 transversely directed, 0.10-inch diameter orifices. The orifices 
were arranged in 32 directly opposing pairs (each of the 32 impingement 
points located approximately at the centroid of an area equal to l/32 of 
the total chamber cross section) on four diameters separated by 0.4 inch, 
and successively rotated 45 degrees, as indicated in Fig. 117. 

Estimation of the heat transfer rate to the hydrogen injector face and 
of the cooling capability of the hydrogen flow through the injector 
orifices indicated marginal face cooling in four injector areas where a 
relatively large gap between injection orifices exists. To improve this 
situation, eight 0.052-inch-diameter holes were provided for additional 
cooling on the upstream edge of the injector. 

A copper spacer similar to the hydrogen injector, fitted with pressure 
taps but without orifices, was fabricated for use between the combustion 
chamber and the nozzle section in the Fg/Li firings. 


MIXinG CHAMBER 

The function of the mixing chamber is to ensure that the hydrogen is uni- 
formly mixed with the F^Li combustion products and that thermal equili- 
brium is attained. The material of construction was OFHC copper. Copper 
is suitable here because the temperature decreases markedly (from —9800 to 
~4000 R) after hydrogen is added. 
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NOZZLE SECTION 
Selection of Throat Areas 


The experimental firing program plan included variations in hydrogen flow- 
rate over the range 15 to 40 percent of total flow. The planned variations 
in propellant flowrates in these firings required corresponding variations 
in throat area to maintain desired chamber pressures (Fig. 118). Ideally, 
a separate nozzle with the proper throat area would be constructed for each 
test condition. However, it is practical to utilize only a small number of 
nozzles to cover the range of test conditions. At some test conditions, 
therefore, chamber pressure or flowrate will be slightly different from 
the nominal values because the exact size nozzle is not used. 

The criteria for selection of the nozzle throat areas were maintenance of 
chamber pressure at 500 psia and of total fluorine plus lithium flowrate 
at 3.74 lb/sec, corresponding to 1.0 lb/sec lithium flowrate at the design 
Fg/ki mixture ratio (2.74). The 500— psia chamber pressure criterion conformed 
to program plan requirements. The constant lithium flowrate criterion re- 
sulted from the need to maintain high lithium injector pressure drop and 
injection velocity to promote a high degree of atomization by the gaseous 
fluorine. 

At a test condition consisting of a given Fg/Li ratio and hydrogen percent- 
age, assumption of constant characteristic exhaust velocity results in the 
following proportionality: 
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Furthermore, since the F 2 /Li ratio and percentage hydrogen are fixed for 
a given test condition, the total flowrate is proportional to the lithium 
flowrate, and the proportionality above becomes 
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A sufficient number of throats should be used such that the fraction 
P «Al,i would not vai T significantly from its nominal value. For the 
test conditions required, the maximum variation in P c /W Li was calculated 
as a function of the number of throat areas to be constructed. These 
results are presented in Fig. 119, which shows that if two throat areas 
were used, the variation in P^. f ro m the nominal value would be 
approximately ±24 percent. Distributed evenly between chamber pressure 
and lithium flowrate, this means that the entire range of conditions 
could be tested within approximately a ± 12-percent variation in chamber 
pressure about the nominal 500 psia, and a ± 12-percent variation in 
W Li about the nominal 1.0 lb/sec. If three throats were constructed, 
P c^Li can 8ta y wi *hin ±16.5 percent of the nominal value. A consider- 
able gain (7.5 percent) is obtained in changing from two to three throats. 
Kith four throats, P C /V L1 can stay within ±12.3 percent of the nominal 
value, a further gain of only 4.2 percent. 


Three nozzles were fabricated; throat dimensions were as follows: 


Nozzle 


Throat Diameter, 
inches 


Throat Area. 
in. 2 

Contraction Ratio 

1.72 

11.5 

2.38 

8.3 

3.30 

5.9 
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1.48 

1.74 

2.05 
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The range for each of the three nozzles is shown in Fig. 120, Each has 
a ±16.5 percent variation of P c /^ Li from the nominal value when used at 
the extremes of its range; when used at the middle portions of the ranges, 
the variation in P c /^ Li from nominal would be less and, in fact, would at 
one point be zero. 

Selection of Nozzle Material 

The criterion to be met in the selection of nozzle material was the 
requirement that test duration should be at least 4.0 seconds (preferably 
longer to allow greater test sequence flexibility). Candidate nozzle 
materials, on the basis of substantial Rocketdyne experience with flu- 
orinated oxidizers, were copper and graphite. 

A transient heat transfer analysis was carried out on the nozzle; the 
Bartz equation (Ref. 33) was used to calculate throat heat transfer 
coefficients, and the techniques given in Ref. 34 were used to obtain 
the transient solution of the time-dependent heat conduction equation. 

It was assumed for these calculations that initial temperature of the 
copper nozzle was 100 F and that maximum run duration corresponded to 
the time required for the inside surface temperature to reach 1250 F. 

Figure 121, which summarizes the results for the copper nozzle, shows 
that maximum run duration would be about 2.2 seconds. Since this is consider 
ably below the desired level, copper was rejected as a possible nozzle 
material. The variation of allowable operating time of the copper 
nozzle over the total range of hydrogen addition is quite small due to 
the fact that throat heat flux decreases by only about 35 percent as 
the amount of hydrogen increases from 0 to 40 percent and chamber tem- 
perature decreases correspondingly from 9800 to 3500 R. Heat flux 
remains relatively constant despite the wide variation in gas temperature 
primarily because of the compensating change in gas heat capacity 
(o. 40 B/lb-R to 1.93 B/lb-R). 
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Rejection of copper as the nozzle material pointed to the selection of 
graphite in its place. The graphite nozzle heat transfer analysis was 
the same as for the copper nozzle except that initial temperature was 
assumed to be 500 F and final inside surface temperature to be 4000 F, 

The first assumption is conservative, although hot gaseous fluorine from 
the fluorine gas generator would raise the nozzle temperature above 
ambient. Figure 122 presents the transient heat transfer analysis results 
for the cases of 0-, 10-, and 20-percent hydrogen addition. Hydrogen 
addition of 10 percent or more easily allows run durations in excess of 
4.0 seconds at chamber pressure of 500 psia; however, only short runs 
are possible when no hydrogen is added at 500-psia chamber pressure. 

Fig. 123, obtained from the same analysis, shows that lower initial 
graphite temperatures result in only relatively small increases in allow- 
able run duration, even when the nozzle is prechilled to -200 F. 

The nozzle section of the thrust chamber (Fig. 99) consisted of an 
ATJ graphite liner enclosed in a steel shell. A compressible, porous 
graphite washer and an asbestos outer liner were included, as in the 
chamber sections, to accommodate thermal stresses. The steel shell was 
tapered to retain and seal the liner. 
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EXPERIMENTAL EVALUATION OF FLUORINE-RICH 
GAS GENERATOR 

A fluorine-rich gas generator was designed to provide the high-velocity 
gases which accomplish liquid lithium atomization. The nominal design 
operating conditions were 700-psia chamber pressure and 500:1 F /H 
mixture ratio. Since this was a new type of gas generator, it ^con- 
sidered necessary to carry out a series of experiments to ensure its 
proper operation. 

Accordingly, sixteen experiments were conducted, with the following 
objectives: 

1. To demonstrate feasibility; 

2. To investigate the effects of varying operating conditions on 
hardware durability and gas temperature. The operating 
variables were to be Fg/Hg mixture ratio, chamber pressure, 
and duration; 

?. To investigate the effects of thrust chamber geometry and 

thermal condition. The geometric variables to be studied were 
chamber length and the presence of mixing plates. The thermal 
condition of the exit plate (which simulated the stainless- 
steel triplet injector) was varied by preheating this com- 
ponent. This variable was investigated primarily because 

of the dependency of fluorine-steel compatibility on metal 
temperature. 


The gas generator experimental conditions and performance results are 
sunnarized in Table 9 . Performance was monitored by measuring chamber 
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pressure and propellant flowrates and computing the resultant charac- 
teristic exhaust velocities and also by measurement of gas temperatures 
with thermocouples inserted into the combustion chamber through the 
mating flange just upstream of the exit plate. Three thermocouples 
were located in the same plane at differing radial positions: (l) center- 
line of the chamber (2 inches from the wall) ( 2 ) 1-1/4 inches from the 
wall, and ( 3 ) l/h inch from the vail. 


"he first three tests were shcrt-duration runs intended to check out the 
hardware and design concept. These firings conclusively established the 
feasibility of utilizing a fluorine gas generator to provide a hot fluid 
for atomizing liquid lithium jets. Combustion gas temperature was not 
measured, however, due to thermocouple burnouts. Prior to the second 
series of tests, a mixing plate was installed in the first chamber seg- 
ment to provide a more uniform radial gas temperature profile. 

The second test series (runs k to 6) was intended to verify long-duration 
capability of the gas generator ac design mixture ratio and to provide 
reliable measurements of c* performance and combustion gas temperature. 

A maximum duration of 10.6 seconds was attained and no hardware damage re- 
sulted. Performance measurements consistently exceeded theoretical values, 
but this may he attributed to imprecision of hydrogen flowrate measurements, 
reflected in the determination of mixture ratio. Measured gas temperature 
is the average of three thermocouple readings, and for tests 5 and 6 com- 
pares very favorably with theoretical predictions. These results indicate 
that the hydrogen reacted completely and that all of the liquid fluorine 
was vaporized. 

The first six gas generator tests employed a 10-inch long chamber. In 
the first three, no gas mixer was installed in the first 5-inch chamber 
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segment. Three additional tests (7 through 9) were carried out in which 
only the 5-inch chamber (containing the mixer) was employed to determine 
combustion efficiency in the short chamber; also, lower mixture ratios 
were used to obtain higher gas temperatures. No difficulties were experi 
enced in these 10-second firings, and there was no hardware damage. 
Performance efficiency remained at 100 percent within the precision lim- 
itation of the very low hydrogen flowrate measurements. In the 5-inch 
chamber, measured gas temperature downstream of the mixer showed approxi- 
mately a 20-percent spread from the lowest value (at the center of the 
chamber) to the highest (at a point l/k inch from the wall); this com- 
pares to an average 10-percent spread in the 10-inch chamber. In all 
cases, average gas temperature is close to theoretical. A mixer plate 
was then installed in the second chamber section also, and subsequent 
radial variations were less than 5 percent. 

The effects of chamber pressure and mixture ratio on hardware durability 
and gas temperature were determined during tests 10 through 16. The 
chamber pressure range was approximately 350 to 720 psia; the mixture 
ratio range was about 370 to 660 . In addition, during the last three 
tests (l4 through 16), the exit plate was preheated to about 500 F 
(simulating the triplet injector firing condition) to determine whether 
a fluorine-steel compatibility problem would occur. 

As indicated in Table 9 , the performance of the gas generator followed 
the theoretical trends in terms of changes in temperature with mixture 
ratio over the ranges of chamber pressure and mixture ratio covered and 
no hardware damage was observed. The final resulting plot of temperature 
vs mixture ratio is presented in Pig. 124 t which shows excellent agree- 

•ment between experimentally observed temperatures and the theoretical 
values. 
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These experiments demonstrated the gas generator chamber geometry and 
operating conditions required for the production of (essentially) gaseous 
fluorine at high temperatures and pressures over long durations, and 
hence its feasibility for the F^Li/H firings. 
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LIQUID LITHIUM FACILITY AND FLOW CHARACTERIZATION 


i-i. 
’ 1 



The objectives of this s. btask were first to design and build a facility 
for storing, heating, and delivering liquid lithium to the experimental 
thrust chambers, and then to conduct a series of lithium flow tests to 
check the operation of the various system components and to characterize 
the flow of liquid lithium through small orifices. 

LIQUID LITHIUM FACILITY 

Design, installation, and operational checkout of a liquid lithium sys- 
tem for use in experimental rocket motor firings constituted an impor- 
tant part of the overall program. Major system requirements were as 
follows : 



1. Capability for storage of at least 20 pounds of liquid lithium 

2. Capability for system heating, to melt solid lithium and to 
maintain the liquid at equilibrium temperature of 500F 

3* Capability for liquid lithium flow control and flowrate 
measurement 

Provision of a heated helium purge gas subsystem 
5. Capability of pressurization tc 1000 psi with filtered, 
moisture- and oxygen-free helium 

Austenitic stainless steel (30^) was selected as the basic structural 
material for the lithium system because it is compatible with hot liquid 
lithium, has suitable structural and thermal properties, and is readily 
available and machineable. The installed system, illustrated sche- 
matically in Fig 125, consists of the following major components* 



83 


Lithium Tank 

This is a 4-foot length of 6-inch, schedule 80, 304 stainless-steel 
pipe, with welded, ASA, ring-sealed, flanged ends. The tank volume is 
approximately 6.7 gallons. It is ASME coded and approved for 1145 psi 
at 600 F. A boss in the tank wall near the lower flange accommodates 
an immersion thermocouple. The upper flange contains a single port to 
which both the pressurization and vent systems are joined, and the 
lower flange has a single outlet port. 


Valves 

The three valves (lithium main, purge, and vent) which might contact 
liquid lithium are special, 1-inch, bellows-sealed , globe valves con- 
structed entirely of 300-series stainless steels*. The valve seat is 
machined directly into the body to eliminate the need for a valve seat 
insert. The other two system valves (pressurizing valve and freeze 
water valve) are conventional Annin globe valves. 


Flow Line 

The 1-inch line (304 stainless steel, schedule 10) between the lithium 
tank and main valve has welded joints to eliminate leakage hazards, and 
is U-shaped to compensate for thermal changes. An in-line, stainless- 
steel screen filter (0.018-inch mesh) Is installed in the main flow- 
line downstream of the main valve and purge line connection point. 


♦Manufactured by Wm. Powell Co., Cincinnati, Ohio 
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Flowmeter 


Lithium flowrate is measured by a magnetic flowmeter (MSA Model FM-4) 
welded into the flow line. The flowmeter consists essentially of a 
l— inch, schedule 10, 304 stainless— steel pipe within a permanent 
890-gauss magnet. The d-c voltage generated by flow of lithium through 
the line is a direct measure of the flowrate. Suitable heating elements, 
a chrome 1-alumel thermocouple, and complete insulation are integral 
parts of the flowmeter unit. 


The original design of the lithium system included a turbine flowmeter 
in the line between the magnetic flowmeter and the main valve. During 
preliminary tests, however, reliable flowrate measurements could not 
be obtained with the turbine meter, hence it was not used in the present 
investigation. 


Freeze Jacket 

An emergency "freeze" jacket around an 8-inch section of the line, with 
inlet and outlet ports for water flow, is included in the lithium line 
downstream of the magnetic flowmeter. In the event of line or valve 
failure downstream of the jacket during lithium flow, a plug of solid 
lithium can be formed in the line by water-cooling of the jacketed 
section. 


Helium Purification System 


To avoid introduction of significant amounts of impurities (such os 
water, oxygen, or particulate matter) into the lithium system with the 
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pressurant gas, a dryer (artificial Zeolite molecular sieve), a de- 
oxygenator (heated copper turnings), and a 10-micron filter are included 
in the helium supply line. 


Check Valves 


Specially designed stainless-steel ball check valves incorporating Inconel 
tension springs downstream of possible lithium flow are used in the 
lithium system. 


Purge Gas System 


The purge gas heater consists of a 4-foot length of 2-inch-diameter 
schedule 40, stainless-steel pipe filled with 0.25-inch steel balls. 
It supplies 500 F minimum helium for at least 10 seconds at nominal 
flowrates when initially heated to about 1000 F. 


Heating System 

Heating of the lithium system is accomplished by means of externally 
mounted electrical resistance heaters. Radiant-type heaters are used 
for the tank (Lindberg "Hevi-Duty", Model 50731) and for the purge 
gas line (Lindberg Model 50031). The lower tank flange, flow line and 
main valve are heated by wrap-around cable (chromalox "Thermwire" , 

Type TVM-20) embedded in "Thermon" heat conducting cement, and the 
lower tank flange is also heated by "Calrod" cartridge units. The 
helium de-oxygenator is heated by two 500-watt Watlow "Firerod' 1 
cartridge units. 
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For control purposes the heater electrical system is divided into four 
sections: one 230-volt, 3-phase, 12 kw source and three 120-volt, 

single-phase, 3 kw sources. The power from each source is varied 
independently by means of manually controlled variable transformers. 
Heatup rates are reproduced by current control through each transformer, 
The power sources terminate in plug-in connectors at the test stand to 
permit convenient tie-in of the heating units. 

Thermocouples are used to monitor temperatures at various points of 
the lithium system. These include two immersed in the liquid (one in 
the tank and one immediately upstream of the injector) and nine 
externally welded to the tank body, lower tank flange, magnetic 
flowmeter, flow line (2), main valve, purge gas heater, purge gas line, 
and delivery line to the injector. Two of the thermocouples are used 
for automatic temperature control of the magnetic flowmeter (Foxboro 
"Rotax" controller) and the main valve (Wheelco controller). 

The heated portion of the lithium system is insulated with fiber- 
glass insulation overwrapped with aluminum foil. The entire system 
is mounted on a frame to make it portable for cleaning or storage. 

The system, with upper tank flange removed, is shown in Fig, 12b and 

127. 


Lithium System Loading and Operation 


Lithium is purchased in the form of 1-pound or 2.7-pound cylindrical 
ingots individually packed in cans in an inert atmosphere. The cans 
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are opened and the tank is loaded with the desired number of ingots in 
an argon atmosphere. A small amount of metallic sodium (0.5 weight percent) 
is added between the ingots to improve the wetting capabilities of the 

lithium (proper line wetting is essential for accurate operation of the 
magnetic flowmeter) . 


System heating begins at the tank and the lower tank flange. After these 
reach approximately 300 F, the line and the main valve are successively 
heated. The purge gas heaters are activated simultaneously with those of 
of the main valve. When the lithium system reaches approximately 500 F 
and the purge gas heater reaches about 1000 F, the power inputs are re- 
duced to maintenance levels. Starting at ambient temperature, about 3 
hours are required to accomplish the heating process. The helium line 
de-oxygenator is heated from ambient to 500 F in about 30 minutes. 



System pressurization and pressure regulation are carried out with 
motorized loaders and dome-type regulators in the customary manner 


LIQUID LITHIUM FLOW TESTS 

A series of lithium flow tests was conducted with two primary objectives: 
(1) to make operational checks of the various lithium system components 
and procedures, and (2) to study the flow of lithium through orifices 
whose diameters and discharge inclinations typified those of a hot-firing 
injector. 
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The experimental series consisted of the following eight tests: 


Test Number 

1 and 2 


3 through 7 
8 


Description 

Flow at tank pressure and flowrate equal to those 
of anticipated hot firings 

Flow through single orifices into argon-filled 
vessels 

Flow through a single, self-impinging doublet 
element into argon-filled vessel 


Lithium System Checkout 


The lithium system loading and operational procedures were found 
to be generally satisfactory. Details of the various subsystem and 
component functions are described below. 


Valves . Operation of the bellows valves was satisfactory although some 
installation difficulty was experienced in obtaining leak-free seals at 
the serrated rings between the upper and lower valve bodies. Opening 
time of the lithium main valve when pressurized to 750 psiwasabout 85 
msec; closing time was about 40 msec. 


Although the main valve functioned satisfactorily, it was found that if 
there was an extended time period (over ten days) between valve actuations 
unpredictable valve behavior occurred. Under these circumstances, exces- 
sive valve opening time delays (l second or greater) were encountered, 
especially during the first attempt at valve opening. This problem was 
circumvented by system heatup on the day preceding the actual test 
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firing, and actuating the valve until the opening and closing times 
were as noted above. 


Heaters. The heaters and heater controls functioned satisfactorily and 
predictably. 

Fl owmeters . Flowrates indicated by the magnetic flowmeter were in good 
agreement with rates calculated from orifice C D values obtained from 
water calibrations and rates estimated from catch tank weight increases. 

The turbine flowmeter produced excellent water calibration results and 
appeared to be functioning properly when it was installed in the lithium 
line. During the flow tests, however, its behavior was very erratic; the 
output would usually peak at the start of the flow, then immediately re- 
turn to zero. Although the meter was not used, its removal from the system 
was impractical and it was left in place. 


P urge Gas System . The purge gas heater functioned satisfactorily. 

Sy stem Operation . Procedures for lithium loading, system pressurization 
and conducting purge gas and lithium flow tests were found to be 
satisfactory. 


Lithium Flow Test Results 


A three-step test procedure was initiated following system warm-up to 525 
±25 F (with purge gas heater at 1000 ±50 F). This consisted of the fol 

lowing: purge gas flow (2 to 5 seconds), lithium flow ( 3 to 9 seconds) 

and purge gas cycling. ’ 
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Te_ 3 ts No. 1 and 2 . These tests were made at lithium flowrate (l lb/sec) 
and tank pressurization level (750 psi) equivalent to those expected 
in the hot firings. Test objectives were to check valve operation and 
lithium flow under these conditions. A l/2-inch line terminating in a 
suitably-sized orifice was positioned downstream of the main valve. In 
Test No. 1, the lithium flow was directed into a covered, 25-gallon steel 
drum in which an argon atmosphere was maintained; in Test No. 2, the 
lithium was allowed to flow out to atmosphere (igniting after emergence 
from the line). No anomalies were encountered in these tests, and it 
was concluded that the system was capable of furnishing lithium at the 
temperature and flowrate needed for the experimental firings. 


Tests No. 5 a . n . d . 7 • These tests were made to obtain photographic 
recording of liquid lithium flow through orifices of various sizes and 
orientations, primarily to aid in providing design criteria for the 
lithium injector. Stainless-steel injector segments were used, each of 
which contained three identical nonimpinging orifices (Pig. 128). Cart- 
ridge heating elements (500 w) were fastened to the test segments to heat 
them to 500 F. The injector segment being tested was placed in a covered 
o x 12 x 26-inch, steel catch tank, between two Lucite windows, one for 
backlighting and one for photography (Fastax camera, 800 to 1000 frames 
per second, Kodak Ektachrome EF-7242 film). A stream of argon was main- 
tained through the catch tank, prior to and during the lithium flows. 


At a 300-psi pressure differential, it was found that steady flow at 
approximately the calculated rates could be maintained through 0.043-inch 
orifices but that sustained steady flow was not obtained through 0.031- 
inch orifices. This may have been due to blockage by solid particles of 
lithium oxide or nitride, since two of the three orifices were found 
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closed after the test. With the 0.043-inch orifices, discharge angles 
of 43 and 90 degrees with respect to the segment face gave free-flowing 
streams, while a 30 degree angle produced noticeable adherence of lithium 
to the face. A typical lithium flow test (Segment No. 1 3 , d = 0.043 in, 
ot= 90 degrees) is shown in Fig 129. 

As a result of these tests und of lithium injection requirements for 
minimum orifice sizes consistent with fuil. non-adhering flow at AP of 
about 300 psi , orifices of 0.043-in diameter and 45 degrees inclination 
were selected for the first lithium injector. 


T est No. 8 . The final lithium flow test of the series was made with n 
single, self-impinging doublet element (orifice diameter = 0.043 inch, 
impingement angle = 90 degrees). The purpose of this test was to obtlin 
photographic visualization of the spray formed by the lithium streams, 
for qualitative comparison with that typical of impinging Cerrosafe 
streams, since experimental results derived from the latter were the 
bases for important lithium injector design criterio. The lithium spray 
was directed into an argon-filled vessel fitted with Lucite windows, and 
strobe-lighted motion pictures (1000 frames per second, l microsecond 
flash per frame) were taken. Figure 130 shows enlarged single frames of 
both lithium and Cerrosafe sprays. Similar primary atomization into 
ligaments and large drops is evident, lending support to the application 
of Cerrosafe droplet size correlations to lithium atomization. 
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task II 

EXPERIMENTAL FACILITIES AND PROCEDURES 

TEST FACILITY' 

Since a detailed deacription of the liquid lithium ayatem has been given 
previously, only the fluorine and hydrogen facilities will be discussed 
m the present section. The experimental firings were carried out at the 
Propulsion Research Area of SSFL. A schematic diagram of the fluorine 
system, which is rated at 3000 psi, is shown in Fig. 131. Fluorine gas, 
supplied from a manifolded system of fifteen, 6 -pound, GF„ cylinders, is 
liquefied in an LN^-. jacketed condenser and fed to a jacketed, 13 -gallon, 
cylindrical run tank. A dip-leg tank outlet is used, leading to a jacketed 
shutoff valve and delivery line. A jacketed, three-way, main valve and 
a jacketed line between the valve and the injector ensure delivery of 
liquid fluorine to the oxidizer manifold. 

The hydrogen system (Fig. 130 consists of two branches, one to the gas 
generator (with flowrates of 0.003 to 0.009 lb/sec) and the other to the 
main combustion chamber (0.7 to 2.3 Ib/.ec). Each branch includes three 
parallel valves and venturi meters. Dy successive operation of these 
valves, several discrete hydrogen flowrates, both to the gas generator and 
to the main chamber, con be obtained in the course of a single firing. 

This capability of obtaining several data points per run proved to be of 
substantial value in the course of the experimental program. 


INST MENTATION 

Instrumentation locations are schematically indicated in Fig. 133 . The 
port t ntlor transducers used for the various types of measurement „ r e de- 
scribed below. 
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Thrust 


The thrust chamber mount was supported on flexures which allowed free move- 
ment parallel to the engine axis, restrained in the thrust direction by a 
Baldwin-Lima-Hamilton double-bridge load cell (Model U-382). 

Pressure 

All pressures were measured with bonded strain gage transducers (Taber 
Teledyne Series 206). Chamber pressures were determined at several axial 
locations both in the gas generator and in the main combustor. At each 
location 2, 3, or 4 circumferential pressure taps were manifolded together, 
and the manifold pressure was measured. The other pressure measurements 
shown in Fig. 133 were made with close-coupled transducers at each location. 

Flowrate 

Hydrogen . Hydrogen flowrates were measured with sonic venturi meters 
(Flow-Dyne Model Series V-16). 

~° rit1 * - Tvo turbine flowmeters in series (Fischer-Porter Model RF-1-50) 
were used to measure volumetric fluorine flowrates. 

Temperature 

Fluorine. Liquid fluorine temperature was measured by two shielded, plat- 
inum resistance bulbs (Rosemount Model 1 7 6) immersed in the liquid stream, 
one upstream of the first flowmeter and the other downstream of the second 


Hydroieu. Temperatures in the venturi plena were measured with iron- 
constautan thermocouples. 


9 * 


L 


► 


i 


f 



ti 

!!/ 

U 

I 


ff 

f; 



V 

I 


i 


I 

1 






C hamber . Temperatures in the gas generator chamber and in the wall of 
the combustion chamber were measured with chromel-alumel thermocouples. 

Data Recording 

Pressure, temperature, and flow measurements were recorded on tape during 
each firing by means of a Beckman Model 210 Data Acquisition and Recording 
System. This system acquires analog data from the transducers, which it 
converts to digital form in binary-coded decimal format. The latter are 
recorded on tapes for computer processing. 

The Beckman Data Acquisition Unit sequentially samples the input channels 
at a rate of 5625 samples per second. Programmed computer output consists 
of tables of time vs parameter value (in engineering units), printed out 
as either the instantaneous values at approximately 10-millisecond inter- 
vals during the firing or as average values over prespecified time slices, 
together with calibration factors, prerun and postrun zero readings, and 
related data. The same computed results are also machine plotted and dis- 
played as CRT outputs on appropriately scaled and labeled grids. < 

Primary data recording for these firings was on the Beckman 210 system. 1 

In addition, the following auxiliary recording systems were employed: 

I 

1. An 8-channel Brush Mark 200 recorder was used in conjunction with 
the Beckman unit to establish time intervals for data reduction 
and for quick-look information at the most important parameters. 

2. A CEC, 56-channel, direct-reading oscillograph was useu as backup 
for the Beckman 210 System for indication of possible oscillatory 
combustion and for precise recording of firing-event timing. 

5. Direct-inking graphic recorders were used to set prerun propel- 
lant supply pressures, for monitoring fluorine system chilldown 
and lithium system heating, to provide quick-look information, 
and as secondary backup to the Beckman and oscillograph recorders. 
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CAL I BRAT r ON PROC FDURES 


Transducer calibrations were used not only to obtain appropriate factors 
for test data reduction, but also to develop statistical histories of each 
transducer so that estimates of short- and long-term deviations could be 
made and probable error bands calculated (see Appendixes B and C for de- 
tailed discussions). The calibration methods used for the various types 
of transducers are described below. 


Thrust 


The thrust-measuring load cell was calibrated in-place. A permanently 
mounted, manually operated, hydraulic force cell was employed which de- 
flected the load cell exactly as did the engine, through a yoke-tension 
rod system. Known loads were applied to the force cell through a Morehouse 
compression-type, temperature-compensated, proving ring calibrated by the 
National Bureau of Standards (NBS). 


Pressure 


Pressure transducers were calibrated by mounting them on stand manifolds 
in which pressures were read with high-precision lieise Bourdon-tube gages. 
The latter were calibrated periodically on Ruska dead-weight testers with 
weights traceable to NBS standards. 


Flowrate 


Fl uorine . Calibrations of the turbine flowmeters to obtain 
rates ‘as functions of rotational speeds were made with wale 
of these cycles-per-gal Ion factors to liquid fluorine usage 
application of corrections which allow for the differences 


vo 1 ume f ] ow- 
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and viscosity between water and LF^. The temperature correction (70 to 
-310 F), which is a function of meter material and not of meter size, has 
been estimated as 1.005 (Ref. 35) and 1.009 (Ref. 36), average: 1 . 007 . 

The viscosity correction, which is a boundary layer phenomenon and there- 
fore depends on flowmeter size, was estimated as 0.992 ( 1-inch meter, 

Ref .35), so that the net correction applicable to the water calibration 
factors was (1.007 x 0.992), or 0.999. This was within the readability 
limits of meter output and was therefore considered negligible. Hence, 
volumetric flow factors determined with ambient— temperature water were 
used for liquid fluorine without correction. 


Hydrogen. The sonic venturi meters were calibrated with hydrogen by the 
manufacturer.* Mass flowrate was determined from the following equation 



The calibrations established curves which gave values of K as functions 
of gas pressure and temperature at the venturi inlet. The flow coeffi- 
cient curves were calculated from the ASME equations (Ref. 37). 

Calibration procedures for the pressure and temperature sensors used in 

conjunction with the venturi meters are described in other parts of this 
section. 


*F1o.-Dyne Engineering, Inc., 3701 K. Vickery St., Fort Worth, Te.e. 76107 
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Lithium. The magnetic flowmeter was calibrated for lithium by the manu- 
facturer* to establish the variation of the flowrate factor (gpm/mv output) 
with temperature. A Leeds & Northrup Type K potentiometer was used to 
calibrate the amplification-recording circuits. 


Temperature 

R esistance Thermometers . Resistance of the platinum thermometers used in 
the LF 2 line was converted to millivolt output by a triple-bridge system. 
This was calibrated by substituting a decade resistance box for the sen- 
sor, and setting it at various resistances corresponding to a temperature- 
resistance calibration supplied by the manufacturer** for each instrument. 
These precision platinum resistance sensors have no significant calibra- 
tion drift. They were checked upon receipt by immersion in liquid nitro- 
gen and liquid oxygen; temperature readings were correct within the limits 
of readability. 


Thermocouples . Thermocouples were used on the basis of the standard NBS 

millivolt/temperature tables. Thermocouple recorders were electrically 
calibrated . 


Calibration Frequency 


The thrust load cell was calibrated immediately before every firing. 
Pressure transducers were calibrated at intervals dictated by the Measure- 
ment Analysis Program (Appendix B). The magnetic flowmeter, fluorine 
flowmeters, venturi meters, and resistance thermometers were calibrated 
just before the start of the experimental firings. 


*MSA Research Corp. , Fvans City, Pa. 

♦Mtosemount Engineering Co., 4900 West 7fith Street, 


Minneapolis 24, Minn. 
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FIRING PROCEDURES 


Fluorine System Pas si vat i 


on 


Prior to assembly, fluorine system components were thoroughly cleaned in 
accordance with standard prescribed procedures (Ref. 30). Passivation of 
the assembled system (to main oxidizer valve) by provision of protective 
fluoride films on exposed surfaces was carried out as follows: low- 

pressure gaseous fluorine was introduced into the system and maintained 
for successive 15-minute periods at 5, 10, and 15 psi; finally, 20 psi 
was maintained for several hours. 


The feed line-thrust chamber system downstream of the main valve was 
passivated immediately before firing by flowing fluorine through the sys- 
tem for several short intervals of time. 

Hun Procedure 

Prerun and postrun zeros were taken with the fluorine and lithium inlet 
lines at run temperatures (-310 and 600 F, respectively) to avoid line 
temperature effects on thrust zero-point readings. 

Firings were controlled by an automatic el-ctronic timer, the sequencing 
is indicated in Fig. 134. For the gas generator both startup and shutoff 
were oxidizer-rich, while in the main combustion chamber there was a hydro- 
gen lead (over the lithium) at startup and a corresponding hydrogen lag 
at shutoff. The purpose of the preliminary hot helium purge was to heat 
the lithium injector tubes. The three step. i„ the main hydrogen flow 

indicate operation of the three main valve, for attainment of three dis- 
crete run levels. 
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EXPERIMENTAL EVALUATION OF 
LITHIUM- FLUORINE-HYDROGEN COMBUSTION 


The original Task II experimental program plan consisted of two phases: 
Phase I, involving F^/Li injector design characterization and Phase II, 
involving evaluation of tripropellant performance. 

The objectives of Phase I were to confirm that performance efficiency was 
a function of the degree of lithium atomization and to define injector 
design criteria in terms of F 0 gas jet momentum and orifice size under 
actual hot-firing conditions. These data were to be used for final design 
of an injector with demonstrated high performance so that if tripropellant 
inefficiencies were observed, the losses would not be erroneouslv attri- 
buted to incomplete formation of lithium fluoride. 


The objective of Phase II was to evaluate tripropellant performance as a 
function of chamber volume, percent hydrogen added, and hydrogen injection 
station. These tests were to be conducted at optimum F 0 /Li mixture ratio 
(2.7'i) and nominal chamber pressure of 500 psia. Thus, the primary objec- 
tive of Phase II was to determine the efficiency of the hydrogen-combustion 
gas mixing process as a function of the above-mentioned variables. 


However, because of operational difficulties and initial hardware damage, 
no significant data were obtained in Phase I and it was necessary to pro- 
ceed directly to Phase II without the injector design information desired 
from Phase I. However, the lithium/fluorine injector design that was em- 
ployed, which was based entirely upon cold-flow and analytical results, 
proved to be satisfactory. It produced nearly theoretical F^/Li combustion 
efficiency. 


Each of the eleven firings conducted is separately discussed in the sectic 
below. This is followed by a general discussion and graphical display of 
+ he experimental results. 
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EXPERIMENTAL RESULTS 


Individual Test Descriptions. Phase I 


Run No. l. The objectives of this test were to check out the facility 
and to evaluate the F 2 /Li triplet injector design. Target conditions and 
parameters for this firing were as follows: 


Chamber Length, inches = 
(injector face to start 
of nozzle convergence) 


W Ip , lb/sec 


2 

V|j^ , lb/sec 

= 

W GH , lb/sec 

r s 

2 

MR (F 2 /1I 2 ) 

= 

MR (Fg/Li) 

= 

T R 

GG’ 

= 

P GG’ psia 

= 

P £ , psia 

= 

F, lbf 

_ 


17 

2.74 

1.00 

0.0055 

495 

2.74 

960 

370 

235 (90 percent 77 „) 
c* 

1045 (90 percent 7?^. ) 


The planned firing sequence was to consist of: (l) gas generator operation 
for 3 seconds, (2) combined gas generator and F Q /Li firing for 2 seconds, 
and (3) an additional 2.5 seconds of gas generator operation only. Both 
the lithium purge (hot helium) and a GN 2 purge through the chamber pres- 
sure taps were to be on throughout the firing, except during portion (2). 


Prior to the firing, the fluorine manifold and injector were passivated 
by pressurizing the LF 2 tank and opening the fluorine main valve for 
three 1-second intervals. The firing was initiated with an oxidizer-rich 
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startup of t lie gas generator; fluorine flowrate (3.2 lb/sec), hydrogen 
flowrate (0.006 lb/sec), gas generator chamber pressure (355 psia), and 
temperature (910 R) were close to targeted values when steady-state con- 
ditions were reached (in seconds). After 3 seconds of gas generator 
operation, during which a hot helium purge was maintained through the I 

lithium system, the lithium purge valve was closed. The lithium main 
valve was programmed to open simultaneously. However, there was no indi- 
cation that it did open; the lower valve stem-travel micrcswitch was not 
actuated nor was any lithium flow shown by the magnetic flowmeter. For 
2 seconds, therefore, the lithium purge was off and the gas generator was 
in operation without mainstage combustion. At the end of this time the 
lithium purge valve was opened (with the programmed simultaneous closing 
of the lithium main valve), whereupon a burnthrough occurred in the lithium 
injector and the test was manually terminated. j 

Significant damage from the burnthrough was confined to the lithium in- I 

jector, which was irreparable. Examination of the lithium injector showed 
that about one-third of the face was burned completely through to the back, 
and that reaction probably began in the manifold. On the part that was 
not burned through, the fluorine orifices were eroded into cone-shaped 
openings with bases at the face, indicating that this burning progressed 
from the gas side into the main injector body. When the lithium delivery 
line was removed^, traces of lithium nitride were found near its exit, 
indicating that a small amount of lithium had passed through the lithium 
main valve, even though the travel microswitches on the valve stem had 

not been actuated and the magnetic flowmeter showed no discernible lithium 
flow. 



The following appears to be the most probable course of events preceding 
the injector burnthrough: when the lithium purge was shut off and the 
lithium main valve did not simultaneously open, gaseous fluorine diffused 
into the lithium injector manifold; the lithium valve opening was delayed 
until just before the signal to close was given so that the stem was barely 
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lilted from its seat; the very small amount of lithium which passed through 
the valve was atomized by the helium purge which immediately followed it 
and was carried into the lithium manifold, where ignition with the fluorine 
already present resulted in the burnthrough. 

Several postrun checks of the lithium valve were made to investigate its 
operation since its physical removal for inspection was impractical. In 
the earlier lithium flow tests, time elapsed from open signal impulse until 
actuation of the lower valve stem microswitch (indicating start of opening) 
was -120 milliseconds; time between lover and upper microswitch actuations 
(considered to be the valve opening time) was approximately 85 milliseconds. 
This compares with the delay of nearly 2 seconds observed in the F 2 /Li 
firing. To investigate the timing of the valve, it was cycled a number 
of times on the day following the firing. The timings observed during 
these cyclings indicated that several cyclings of the valve were needed 
to establish normal operational timing after a period of inactivity, 
possibly because of initial delays in the functioning of the bellows 
around the valve stem. They also pointed to a requirement for estab- 
lishment of normal valve timing by cycling the valve immediately prior 
to a firing. 


Run_No. 2. The following changes were made in the test procedures to 
avoid possible recurrence of the lithium valve malfunction and injector 
burnthrough: 

1. After warmup of the lithium system to 500 F, the lithium main 
valve was cycled three times (with lithium flow to atmosphere), 
in the last two of which the time elapsed from valve signal to 
valve liftoff (130 milliseconds) and from valve liftoff to open 
position (85 milliseconds) were in the normal range. 

2. Temperatures of the preheated lithium injector and of the helium 
purge gas were reduced from 550 to 450 F and from 1000 to 700 F, 
respectively, to lessen the possibility of ignition of the in- 
jector body by fluori r *. 
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Closing of the lithium purge valve was slaved to actuation of 
the upper microswitch on the stem of the lithium main valve so 
that purge pressure would be maintained until the main valve was 
fully open. 

A sampler circuit was installed to check the lord cell output at 
a progrumned instant in the firing so that the run would be cut 
if thrust were below a preset level. 

5. Automatic run cut could also be triggered by discontinuity in a 
fine wire wound around the oxidizer system and injector such as 
would be caused by a serious fluorine leak. 


Target conditions for the second firing were the same as for run No. 1. 

Decouse of contact malfunction in one of the channels of the electronic 
timer used to sequence firing functions, two attempted firings preceded 
the actual test. The timer malfunction in the two aborted runs resulted 
in premature closing of the fluorine valve during gas generator operation 
prior to opening of the lithium valve. This was followed by automatic 
cutoff by the thrust sampler circuit 0.3 second after opening of the lithium 
valve. Run No. 2 was carried out following rectification of the timer 
malfunction and sequencing occurred as programmed. This firing was man- 
ually terminated after 1.5 seconds of lithium flow when failure of a check 

valve in the fluorine purge system released some fluorine into the test 
pit. 


Postfiring examination of hardware, test data, and motion pictures showed 
that the center portion of the lithium injector was burned out, and that 
this burnout occurred during the first part of the run before lithium 
valve opening. Consequently, no usable performance data were obtained. 
The posttest examination also indicated that the premature closing of the 
fluorine valve in the preceding two firing attempts resulted in lugh- 
temperoture transients which subjected the gas generator chamber and the 
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li tli tun injector to substantial overheating. Although no external hard- 
vare damage was visible after these attempts, close inspection was not pos- 
stble because of t!ie presence of lithium in the chamber and nozzle as well 
as on the test pit floor. Hence, the condition of the lithium injector 
after the aborted tests is not known. 


Hun_No_. 2 - Test objectives and target run conditions for the third firing 
were the same as for the first two. In this firing, however, the graphite- 
core, doublet-pattern lithium injector was used because a triplet was not 
available. 

Prerun preparations included the customary lithium system heatup and cy- 
cling of the lithium main valve to ensure proper functional timing. It 
was observed at test start, during operation of the gas generator, that 
the fluorine flowrate was significantly below the expected level. The 
run was therefore manually cut after 2.1 seconds of F., II., firing (prior 
to opening of the lithium main valve). No hordwure damage occurted, 

t'xarni nation of the liquid nitrogen system used to chill the fluorine sys- 
tem uud a series of blowdowns using liquid oxygen as fluorine simulant 
indicated that the liquid nitrogen flowrate through the ja -ket around the 
fluorine line between the tank and main valve was not adequute on excep- 
tionally hot days to prevent partial gnssi fication of fluorine within the 
line. Suitable modifications were accordingly made in the liquid nitrogen 
system to permit attainment of a sufficient degree of line chill even 
during hot-weatner conditions. 

SHE No » !* • This was the first firing in which F Li combustion data were 
obtained. Target run conditions and test sequence were the same a- in 
previous runs; tie graphite-core, doublet-pattern injector was used. 
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The firing was cut by the thrust sampler circuit at the 3.1-second point, 
after 0.8 second of lithium flow, because thrust had not reached the preset 
minimum level. This was due to low lithium flowrate, as shown i„ the tab- 
ulation below. The table gives expected and observed levels of operating 
parameters during both the preliminary gas generator operation and the 
*2 Li firing. The agreement between targeted and experimental values is 

good in all cases except for lithium flowrate and, hence, for the Li 
mixture ratio. 


Parameter 


Liquid Fluorine Flowrate, 
lb/ see 

Gaseous Hydrogen Flowrate, 
lb/sec 

F.,, II ( Mixture Ratio 

Liquid Lithium Flowrate, 
lb/sec 

F <v /Li Mixture Ratio 


Chamber Pressure, psia 

Thrust, pounds 

experimental Character- 
istic Velocity, ft/sec 

Theoretical Character- 
istic Velocity, ft/sec 

Cncorrected Character- 
istic Velocity efficiency, 
percent 


Gas Generator Alone 

1 Gas Generator and 
1 Lithium Flow 

Target 

Fxperi mental 

Target 

experimental 

*-*.97 

2.85 

r) -» t 

2.60 

0.0033 

0 0035 

ir* 

ir* 

C 

o 

o 

0.0055 

537 

528 

593 

582 



1.00 

0.25 



2.75 

10.5 

290 

303 

315 

311 

o 

X 

830 

900 




235 

109 



1053 

500 



5090 



55 

00 
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The graphite injector body was intact after the firing, and the fluorine 
.ri fires were not significantly enlarged. Six of the eight lithium de- 
livi.y tubes in the injector were found to be plugged with solidified 
lithium after the firing. Analysis of the contents of the plugged tubes 

showed n significant degree of lithium contamination (epproxl raately 5 
per~ent ). 
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To remove the contaminants, which had accumulated during repeated contact 
of molten lithium with impurities, the lithium system was emptied and 
thoroughly cleaned. In addition, a series of tests was carried out to 
determine whether the blockage resulted from clogging of the tubes by con- 
taminants in the lithium or from inadequate heating of the tubes prior to 
the firing. Blowdowns with both hot helium and lithium were made through 
the ring manifold and lithium injector tubes. The manifold was preheated 
electrically, while the delivery tubes were preheated prior to lithium 
flow with a heated helium purge. These experiments indicated that the 
tubes were adequately heated by the hot purge gas under firing conditions 
( 1-second flow of 1000 F helium). At this time, an in-line, stainless- 
steel, screen filter (0,018-inch mesh) was installed in the lithium delivery 
line to remove any residual impurities. 

Although lithium flow was delivered through only two of the eight injector 
tubes, which resulted in poor propellant distribution, I., Ii combustion 
efficiency was over percent, indicating promise of high combustion ef- 
ficiency under design operating conditions. 


Individual Test Descriptions. Phase f r 

- ^ - n Xo - 1 - At this P Mint tl><* program, the decision was made to proceed 
to the Phase II tripropellant firings, in which the graphite doublet in- 
jector would he used. It was recognized that the previous failures of the 
triplet injector were attributable to facility malfunctions rather than 
to its design. However, the doublet injector appeared to be capable of 
performing quite adequately and was apparently more durable. Target test 
conditions for run No. *; were as follows: 

l’ c ">00 psin 

F, I.i Ml 2.7* 

H.j Addition, (percent of total flow) - 40 
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The thrust chamber configuration employed (Fig. 11")) was that shown sche- 
matically in Fig. 9*), i.e., go s generator chamber, doublet injector, 20-inch ' , 

combustion chamber, hydrogen injector, 10 -inch mixing chamber, and nozzle. 


During the first second of gas generator operation, the fluorine flowrate 
was nearly at the targeted value, but it then declined to about half the 
targeted value for the remainder of the firing. Posttest examination of 
run data and procedures indicated inadequate chilling of the fluorine in- 
jector, with consequent tvo-phase flow through the orifices, as the cause 
of the low fluorine flowrate. The insufficient chilling was attributed 
to the GN,, purge maintained through the fluorine injector during the chill- 
down period. The pretest purging procedure was modified for later firings. 

The lithium flowrate was at the targeted level during the first 80 milli- 
seconds of mainstage, alter which it increased significantly to a point 
substantially higher than the calibration range of the magnetic flowmeter. 
This increase in lithium flowrate was due to burning of the orifice tips 
of the lithium delivery tubes and to low chamber pressure resulting from 
the decreased fluorine flowrate. In spite of the unfavorable degree of 
lithium atomisation and the poor distribution resulting from the very low 

nixture ratio ( 1 . 0 )), corrected c* efficiency was approximately 9 /, 
percent. 

Hardware damage resulting from this firing was relatively minor and in- 
cluded burning of the inner core of the main hydrogen injector (copper) 
and hydrogen mixing chamber (copper), erosion of the graphite core of the 
lithium injector, and burning of the stainless-steel gas generator chamber 
wall immediately upstream of the lithium injector. 
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R - n - No - *? • Targeted parameters for this firing were the same as those for 
Run No. ~). Hard ..are Modifications prior to the test included substitution 
of a graphite (for copper) hydrogen injector core (Fig. 1 %), a graphite liner 
in the copper mixing section, and a graphite mixer plate and liner in the 
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downstream half of the second gas generator chamber (nearest the lithium 
injector). In addition, the lithium injection tubes v/ere recessed about 
l/8 inch from the face of the injector to provide a measure of protection 
from the hot chamber gases. 

The test was completed successfully. Observed parameter values are com- 
pared to the corresponding targets in the following tabulation. Details 
ol the computational procedures and corrections applied are given in 
Appendix A. 


Run No. 6 

Target 

Actual 

lb/ sec 

2.52 

2.59 

W LLi’ lb/aec 

0.94 

1.08 

W G|[ (main), lb/sec 

2.26 

2.15 

Wqjj (gas generator), lb/sec 

0.0055 

0.0055 

F 0 /Li MR 

2.68 

2.40 

II (> , percent 

40 

37.0 

Gas Generator Pressure, psia 

535 

489 

Chamber Pressure, psia 

490 

449 

Thrust, pounds 

2360 

2185 

r l c +(co l rr), percent 


100.0 

TJj (meas), percent 
s 


91.5 


Posttest hardware examination showed that the orificed ends of the lithium 
injection tubes were burned away and that there was substantial erosion 
of the stainless-steel portion of the gas generator upstream of the graphite 
insert. The burning of stainless steel adjacent to graphite had been gen- 
erally observed during these firings. This suggested initiation of 
fluorine-steel reaction by the heat of surface reactions of fluorine with 
graphite. To minimize the possibility of stainless steel burning in sub- 
sequent firings, nickel liners were installed in the first gas generator 
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cuanber to act as a buffer between the graphite in the downstream chamber 
and the stainless-steel fluorine injector. 

him No. 7. Thrust chamber geometry and target run conditions for Run So. 7 
'ere the same as for Run No. 6 except for the percent hydrogen added. In T 

this case, two hydrogen flowrates were planned during the test, by use of j 

the parallel hydrogen feed system previously described. In addition, the j 

lithium injection tubes were recessed further from the injector lace (about ! 

1/2 inch) to increase their service time beyond 1 second | 

Both of the planned steps were achieved and c* efficiency remained high: 4 

step 1, II 2 = 2 k percent, corrected c* efficiency = 100 percent; step 2, 

H o = 20 percent, corrected c* efficiency p)l percent. 

RunJlo._8. The test objective of Run No. 8 was to determine the influence 

of hydrogen injection station on performance. Accordingly, the hydrogen 

injector was moved forward 11 inches. The resulting chamber configuration 

consisted of (l) the gas-generator plus lithium injector, (2) a 1-inch- ^ 

long F^/Li combustion zone (previous length was 20 inches), (l) the hydro- * 

gen injector, and ('») a 27-inch-long mixing chamber (previous length was 

12 inches). The target run conditions included three discrete hydrogen 

flowrates, corresponding to 11, 21, and 11 percent of total propellant 

flowrate. j 

4 

1 

The test was successfully completed and data points for all three hydrogen 5 

flow conditions were obtained. Performance remained liigli (rj >07 percent) 

— ‘ - n - No - 2 - The test objective of Run No. 0 was to determine the effect 
of chamber length on performance. To accomplish this, the downstream mixing 
chamber length used in Run No. 8 was reduced from 27 to 17 inches. The 
target run conditions and remaining thrust chamber geometry remained un- 
changed. The test was successfully completed; three hydrogen flow conditions 


111 




were evaluated * However, the measured lithium flowrates indicated that 
four of the eight injector tubes had plugged either at or prior to test 
start, because the indicated flowrate, although constant, was only about 
00 percent of the target value. In spite of the resulting high F 0 /hi 
mixture ratio and attendant nonuniform mixture ratio distribution (i.e. , 

lithium was not injected into four of the GF jets), performance remained 
high. 


Hun s No. 10 and 11 . Subsequent flow checks and two additional attempted 
firings (Runs No. 10 and 11 ) revealed that the lithium main value was the 
cause of the difficulty. During Run No. 10, the main valve failed to 
open properly (the opening time was about 750 milliseconds compared to 
150 milliseconds observed in previous tests) and all injection tubes were 
plugged. Continued investigation revealed that not only was valve opera- 
tion unpredictable, but significant leakage was occurring across the valve 
seat and becoming more severe with time. Initially, the valve only leaked 
under pressure; later, however, valve leakage was observed as soon as the 
lithium was liquified, even without tank pressurization. 

Since the lithium injecter tubes were heated by hot purge gas, the pres- 
ence of liquid lithium in the run line downstream of the main valve would 
result in tube plugging (lithium would be carried into the injector tubes 
before they could be adequately heated). No additional tests were attempted 

Data Summar y 

Performance levels were established on the basis of c* efficiencies ob- 
tained from measurements of chamber pressure and thrust with applicable 
corrections for energy losses, throat area changes, and departures from 
ideal, invisod. one-dimensional flow. Details of the procedures used 

for these computations and for estimation ol the correction factors are 
given in Appendix A. 
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A measurement error analysis was carried out (Appendix C) to obtain esti- 
mates of the^ uncertainty intervals associated with the determinations of 

L ,7 c*Jp c and [Vi* This anal ysis indicated the following typical error 
bands at the 2 a (95 percent) confidence level: 

F° r fu c# ^ : -1.0 percent 
c 

p c *J F ■’ 11 -3 percent 

A summary of the experimental performance data is presented in Table 10. 

Note that ^ (mens) is the observed specific impulse (corrected to vacuum); 

multiplication of this efficiency by <P TQTAL gives c* efficiency based on 
thrust, 

i 

The values of c* (theo) and I s (theo) were obtained from theoretical per 
formance calculations carried out for each individual set of run condi- 
tions. This avoided errors arising from the triple interpolations that 
would otherwise be required (for P c> MR, ’d percent Hg). 

Th. good agreenen' b.t.o,„[ , c# ] p dddL.t i. 8 h..„ i„ Fig. 1,7, id 
which only minor departures from tSe ideal 45-degree line are evident. 
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DISCUSSION OF EXPERIMENTAL RESULTS 

The nine separate experimental conditions cov?red in runs No. 6 through 
9 provided a total of 18 performance data points, nine each for c* efficiency 
based on chamber pressure and on thrust. Correlations of these data with 
analytical predictions of the degree of lithium atomization, fraction of 
lithium vaporized and burned, and uniformity of propellant distribution 
are presented in this section. 

The lundamental requisite for the practical feasibility of the F^/Li/H 
tripropellant combination is that all, or nearly all, of the injected 
lithium be vaporized and burned. A series of curves showing predicted 
c* efficiency for varying percentages of lithium vaporized as a function 
ol percent hydrogen, at nominal chamber pressure and F^/Li mixture ratio, 
is presented in Fig. 138. They were calculated by the methods discussed 
in preceding sections of this report. Superimposed on these curves are 
the experimental Ly observed c* efficiencies. It is apparent that in all 

cases the injected lithium is essentially completely reacted with the 
fluorine, and the combustion products are uniformly mixed with the hydro- 
gen; hence the fundamental feasibility of the Fg/Li/Hg system is 
established. 

The results illustrated in Fig. 138 also point up the fact that separate 
determination of Fg/L* combustion efficiency, originally planned as Phase I 
of the experimental firing program, was not essential since overall ef- 
ficiencies on the order ol 100 percent were obtained. Consequently, 
bypassing of the Fg Li firings was not detrimental to achievement of pro- 
gram objectives. 

The nine experimental firing conditions cover a range of chamber pressure 
(348* to psiu) within which the predicted effect of chamber pressure 
on c* efficiency is negligibly small (Fig. 88 ). Further, the variation 
of Fg/Li mixture ratio included in these firings ( 2 .'i to 5 . 2 ) also has 
negligible predicted effect on c* efficiency as shown in Fig. 139 . For 
these reasons, the significant variables in the tripropellant firings 
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conducted were: (l) percentage of hydrogen, (2) point of hydrogen addi- 

tion. and ( 3 ) chamber length. 

The cylindrical portion of the experimental thrust chambers (i.e., the 
length between the lithium injector and start of nozzle convergence) 
included the F^/Li combustion chamber section, the 2-inch-wide hydrogen 
injector, end the mixing chamber between the hydrogen injector and the 
nozzle. (For convenience, the hydrogen is considered to be added at a 
single axial point, and the width of the hydrogen injector is included 
in the length of the mixing chamber section). The first six tests used 
a 32— inch chamber. In the first three, hydrogen was added at 20 inches 
from the lithium injector (giving a 20-inch combustion chamber); 

in the second three, hydrogen was added at 5 inches from the lithium 
injector (giving a 5-inch F^/Li combustion chamber). The final three 
tests employed a 22-inch chamber with hydrogen addition at 5 inches. 

Each of these three chamber configurations will be separately discussed 
in the following paragraphs in terms of the degree of lithium atomization 
and hydrogen mixing efficiency. 

The results of the first three tests are plotted in Fig. 140 together 
with the c* efficiencies estimated for lithium mean droplet sizes of 30 
and 70 microns. Corrected c* efficiency was essentially 100 percent 
over the hydrogen addition range of 19 to 37 percent of total propellant 
flowrate. These results indicate that a high degree of lithium atomi- 
zation was obtained (< 30-micron mean droplet size) in all three cases. 

Figure 141 presents the results of the second set of three tests, with 
the hydrogen addition point moved 15 inches upstream. Again, c* effi- 
ciencies of nearly 100 percent were obtained except at the lowest per- 
centage of hydrogen, where c* efficiency was 95-5 percent. The decrease 
in performance at this point is due to incomplete peuetration of the 
F L /Li combustion gas by the hydrogen, as discussed in the design portion 
of this report, with resulting poor distribution. 
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Results of the third set of firings, in which the 5-inch F^/Li combus- 
tion chamber length was retained while the mixing chamber length was 
decreased to 17 inches, are shown in Fig. 142. Again, c* efficiencies 
of essentially 100 percent were obtained except at the lowest percentage 
of hydrogen, where incomplete penetration occurred. The difference in 
c* efficiency at the 19-percent hydrogen level between the 20- and 5-inch 
F^/Li combustion chamber length configurations (Fig. 140 and 142) indi- 
cates that completeness of penetration of the hydrogen at this level is 
marginal and, in addition, that the 20-inch chamber ensured completeness 
of the F^Li reaction and consequent greater combustion gas uniformity at 
the hydrogen injector than did the 5-inch F^/Li chamber. 

Figure 143 shows the effect of F^/Li combustion chamber length on pre- 
dicted c* efficiency (with three total chamber lengths) for a lithium 
mean droplet size of 30 microns and 25-percent hydrogen addition at 
500-psia chamber pressure and 2.74 F^Li mixture ratio. Observed c* 
efficiencies with the three thrust chamber configurations employed are 
also shown. With actual lithium mean droplet sizes indicated to be less 
than 30 microns, predicted c* efficiency at zero F^/Li combustion chamber 
length (i.e., all three propellants injected at the same axial point) 
is over 95 percent even with the 12-inch overall chamber length. 


CONCLUDING REMARKS 


The experimental results of this program clearly demonstrate that essentially 
complete combustion can be obtained with the Fg/Li/lI (} tripropellant combi- 
nation, and confirm the validity of the criteria upon which the designs of 
the thrust chamber componejtn were based. 

Program experience showed that no major difficulties are involved in the 
simultaneous use of a cryogenic and a heated liquid-metal propellant, pro- 
vided careful attention is paid to operating procedures. It was also 
found that adequate material compatibility for the Fg/Li/Hg system may 
be obtained without recourse to exotic materials. 

Although combustion efficiencies were high in the thrust chamber config- 
urations er.ployed, several simplifying modifications remain to be experi- 
mentally investigated, such as: 

1. Use of smaller volume chambers (The results of the present program 
indicate that high c* efficiencies may be obtained with signifi- 
cantly shorter chambers than were employed because of the small 
lithium droplet sizes produced by gas-augmented atomization.) 

2. Use of nonreactive gas for lithium atomization; i.e., use ot a 
fuel-rich, instead of an oxidizer-rich, gas generator (This would 
minimize combustion at vhe injector face and the hazards of pre- 
mature fluorine-lithium interaction.) 

3. Avoidance of difficulties inherent in the use of a liquid-metal 
valve by employing alternate methods of lithium flow control 
(e.g. , rupture disks) 

With the attainability of high combustion efficiency demonstrated, the 
next important step in determining the practical feasibility of the 
F^Li/Hg combination is an experimental investigation of nozzle efficiency 
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because to achieve the high theoretical specific impulse of the tripro- 
pellant system , both c* and C p efficiencies must be high. Such studies 
of the nozzle flow process are particularly essential when, as in the 
case of F^/Li/Hg, expansion of a two-phase system is involved. 
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NOMENCLATURE 


7T 

burned 


injected 


total 


transfer number = C (T - T )/A H 

p B A' v 

characteristic velocity- 

theoretical c* at the mixture ratio corresponding to the 
burned propellants 

theoretical c* at the mixture ratio corresponding to the 
total amount of all propellants injected 


- gas-phase heat capacity 
= drag coefficient or discharge coefficient 
= diameter 

~ volume mean droplet diameter 
= measured droplet image size of i th droplet 
= generalized mean droplet size 
= binary diffusion coefficient of oxidizer 
= dimensionless molar flux of oxidizer 


f».toich/ U ][/ P °" 


p - p (l - a M , . . 

0 stoich- 


total correction factor applicable to measured P 


= thrust 


- latent heat of vaporization 
= flow coefficient for sonic venturi meter 
= burning rate constant corrected for convection 
= stagnation burning rate constant 
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L* 

= 

characteristic chamber length 

zrr 


_ _ _ 

Le 

— 

Lewis number = X / (D P cVrt) 

M 

= 

total number of droplets counted (Eq. 22) 

^stoich 

= 

moles of gaseous product per mole of oxidizer 

N 

= 

number of droplets 

P 

= 

exponent expressing nature of mean droplet size (Eq. 22 

P 

= 

pressure 

P 

c 


chamber pressure 

P 

0 

= 

partial pressure of oxidizer 

Pr 

= 

Prandtl number 

q 

= 

exponent expressing nature of mean droplet size (Eq. 22 

fl 

= 

gas constant 

Re 

= 

Reynolds number = (AV)/jli^ 

t 

= 

time 

fc b 


breakup time (Eq. 18) 

b f 

= 

flight time (Eq. 17) 

P 

= 

temperature 

J 

= 

velocity 

\V 

= 

relative gas-liquid velocity 

1 

S 

= 

molar vaporization (burning) rate of droplet 

r 

= 

dimensionless droplet burning rate = w "c /2 rrX D 

s p A 

r 

= 

weight flowrate 

burned 


total weight flowrate of all propellants injected minus 
weight flowrate of propellants remaining unburned at the 
throat 

injected 

= 

total weight flowrate of all propellants injected 



We 

X 

P 

y 

o , c 

a 

0 

V 

*1 

s 

^stoich 

(v) F 
(v) „ 

c 

0 

A 

M 

P 

O’ 

0 

^total 


= Weber number = D T O (A V) /cr 

L g ' ' L 

= penetration distance of liquid stream into gas jet 
= weight fraction oxidizer in the gaseous phase 
= fraction of product vaporized 
- diffusivity ratio = TT X/(f T 

p 7 p 

= characteristic velocity efficiency 
= specific impulse efficiency 

= stoichiometric coefficient (moles of lithium per mole of 
fluorine) 

“ characteristic velocity efficiency based on measurement of 
thrust 

= characteristic velocity efficiency based on measurement of 
chamber pressure 

= angle between liquid stream and injector face 
= thermal conductivity 
= viscosity 
= density 
= surface tension 
= defined in Eq. 14 

= total correction factor applicable to measured thrust 


Subscripts 


A, B, C, - locations A, B, and C, respectively, of Fig. 77 
D = droplet 

g = gas 


GG 


gas generator 



L = 1 

t - throat 

Superscripts 

= = value in zone BC (Fig. 77) 

= value in zone AB (Fig. 77) 

o = initial condition 


TABLE 1 


^2 ' 153 R ) ' /hLi ( 96 ° R ) pitted theoretical performance parameters 


i Plotted 

Ordinate . Atbc^sss j Expansion Parameter 


Constant 

Parameters 





(153 Bl /LLi (960 R) /GHg (537 R) PLOTTED THEORETICAL PERFORMANCE PARAMETERS 
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TABLE 4 


DATA USED TO CALCULATE LITHIUM DROPLET BURNING RATE 
1 

Term 

Value 

Comments 

r A 

3900 R 

Boiling point of lithium 


r 

at 500 psia (itef. 59) 

X. 

1.9 X 10 Btu 

aec-f t-F 

Thermal conductivity of 
gas in zone AB (Fig. 77) 

— 

P 

, p Btu 
lb-F 

Heat capacity of F Q -Li 

- 


combustion products 


21.2 x 10' 6 - Btu 

sec-f t-F 

Thermal conductivity of 
F 0 -Li combustion products 

^stoich 

2.0 

Moles of Li/ mole F Q 

^stoich 

2.44 

Moles of gaseous product 
per mole of oxidizer 

i 

o, c 

0.25 

Assumed constant at this 
intermediate value 

p 

0.7 Btu/lb-R 

Heat capacity of gas 
in zone AB (Fig. 77) 


7,600 Btu/lb 

Heat of vaporization of 



lithium at 5900 \i (Ref. >yj 

e 

1.0 

Lewis number 


33 





























TABLE 7 

PHYSICAL PROPERTIES OF CERROSAFE* AND LITHIUM 


Property 

Specific Gravity at 68 F 
Melting Range, P 

Liquid Surface Tension at Melting Point, 
dyne/ cm 

Volume Change (Liquid to Solid), percent 

Volume Change (Growth After Solidification), 
percent 

Composition , percent: 

Bismuth 

Lead 

Tin 

Cadmium 


LITHIUM 

Cerrosaf e 

Lithium 

9.4 

0.55 

0 

ON 

1 

o 

■sO 

4 

357 

400 

395 

-2.0 

— 

0.5 

— 


^Manufactured by Cerro Sales Corporation, New York, N. Y. 



TABLE 8 

DROPLET DISTRIBUTION DATA, ATOMIZATION OF CERROSAFE WITH NITROGEN 
IN A LIQUID-GAS -LIQUID TRIPLET ELEMENT 
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DATA SUMMARY 



Append 










Mixture Ratio 


Figure 1 


Variation of c* (shifting Equilibrium) with 
Mixture Ratio at Indicated Chamber Pressures 
I<F n (153 R)/LLi (960 R). 








Chamber Te/nnerature 


✓ V ‘ 



Figure 2. Variation of Chamber Temperature with Mixture Ratio 

at Indicated Chamber Pressures, LF 0 (153 R)/LLi (%0 R). 










Vacuum Specific Impulse, Sec 



Figure 3. Variation of Vacuum Specific Impulse (Shifting 
equilibrium) with Fixture Ratio at Indicated 
Chamber Pressures and Expansion Ratios, 

LF 2 (153 R)/LLi (960 It). 
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Coefficient 







Characteristic Exhaust Velocity 
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Figure 5 . Variation of c* with Mixture Ratio at P = 300 
psia, LF 2 (153 R)/U,i (960 R) , c 

Curve Is Shifting Equilibrium, Ref. 14 
Curve 2: Shifting Fouilibrium, Present, Studv 

” ' Fro?en Composition, Ref. H 

Fro?en Comnosition, Present Study- 


Curve 1 
Curve )i 
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Specific Impulse. Sec 



Thrust Coefficient 
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Figtire 56 . Vacuum Specific Impulse (Frozen Composition) as Function of Fxoansion 
Ratio at Indicated Percent ffydrogen, P = F00 osia, F„/T,i M.R. = ?. 7 ) ; 
LF 2 (153 P)/SLi(36 R)/L^(36 R). c 2 
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Vacuum Thrust Coefficient (Shifting Enuillbrium) 
as function of E^/U fixture Ratio at Indicated 
Percent Hydrogen, £ = 900 osla, e= 900. 
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Figure 68 . 


Variation of Vacuum Thrust Coefficient (Shifting 
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Mechanism A - Small Droolets 
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Pifrure 76. Ignition and Combustion Mechanisms of Liouid Lithium i 


Figure 77. Simplified Schematic of Lithium Droplet Burning in 
lluonne Atmosphere by Vaporization Rate-Limited 
Mechanism. (Lithium vapor diffuses away from the 
liquid droplet (A) to the flame front (B) where 
it reacts with fluorine. The LiF then diffuses 
outward.) 




PERCENT Li VAPORIZED AND REACTED 


re 7> . Calculated Variation of e» Efficiency 
vith Percent Lithium ^aoorized, at 
Indicated Percentages of Hydrogen. 

All Vanorized Prooellant.s Assumed to 
be ComDletely Reacted. 
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Figure 83. Calculated Variation of c* Efficiency with Percentage 
Hvdrogen, Added at ^-inches from F«/Li Injector Face 
in 10-inch Chamber, for Indicated Volume Mean Lithium 
Droplet Diameters (in Microns). 



























Figure 91 . Schematic of Apcaratus Used for Jet Penetration Studies 
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; Figure 95. Atomization of Liquid Cerrosafe by GN 

I 1° Triplet Element. ^ 

Liquid Orifices: D « 0.06()-inch ' 

Gas Orifice: D - 1.50-incb, | 

V - 1,090 ft / sec 





Figure 96. Photomicrograph of Solidified 
Cerrosafe Farticles Collected 
During Test No. 49, 100X. 
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Droplet Diameter, D , Microns 

ire 97. Experimental Droplet Distributions, Atomization of Liquid Cerrosafe. 








Figure 98. Coneeotual Schenatic of If ? /LLt/aH ? Thrust CThamber . 




Thrust Chamber Schematic. 








Figure 100. Variation of Chamber Temperature with 

Mixture Ratio at Indicated Chamber Pressures 
LP 2 (153 R)/GB 2 (537 R). 
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Fif-jure 102 . Gas Generator Assembly, Schematic 
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Figure 107 . Gu Generator Componenta 





ter. Microns 



A. Penetration of IFf Jet by LLl Stream as Function of 

r iF ? Orifice Diameter at Indicated P and T« . 
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8td«, Shovln* Inlet* to GF^ Orifi c#«. 
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5AAV.-4 2'. 67-81 

P, * urp 110. Triplet-Pattern Lithium Pluorine In|ector, Pace Vi»v vlth 
Ablative Shield in Place. 
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5AA3V'«, 21/67-SIG 

I 

i 

Figure 111. Ablative Plate for Triplet-Pattern Lithium/. luorine 
Injector, Fare View. 






5AA33-V'2l/67-SID 

gure 113 . Lithiua/Fluorine Injector, Face View 

Without Ablative Shield. 
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Figure K4. 


?!?5 U r t : P *“' r “ r Uh Wl’I»orio. Injector and Lithiu. Hani- 
fold, Gn> Generator Side, Showing Fluorine Orifice Inlets. 







Figure 113. 




’ Lithiu "Aluorine Injector, Face Viev 
Without Ablative Shield. 





Figure 114. 


U /—OlA 


^?H le r t r Pa r U "“ !' Uh ‘ u "/ Flu ' jri " e Injector and Lithiu. M,„i. 
fold, Gao Generator Side, Showing Fluorine Orifice Inlets. 
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Figure 115. 


Doublet- Pat tern Lithium/Fluorine Injector, Face View Showing 
Showerhead GF^ Orifices and Outlets of Lithium Delivery Tubes. 
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Impingement Points of 32 Pairs of Opposing Hydrogen 
Jets, Schematically Projected on a Single Plane 



. °chematic of One of Four Sets of Eight Pairs of 
Impinging Hydrogen Jets, Each Set in a Plane 
Successively Separated by 0.4 inch, and Success- 
ively Rotated 45 Degrees. 

Figure 117. Hydrogen Injector Orifice Arrangement 









Figure 120. Operating Ranges of Three Norzles. 
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Figure 121. Maximum Run Duration of Cooper No??le as Function 
of Hydrogen Percentage. 

P = 500 osia 
c 

F ? /Li M.R. - 2.7U 

Initial Copper Temper e a 100 F 

Maximum Surface Temne. .cure at Run Termination = 
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Figure 125. Liquid Lithium System, Schematic Flow Diagram 
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Figure 129. 


Liquid Lithium Flowing Through 0. 043-Inch Showerhead 
Orifices at 0.07 lb/sec/orifice, AP = 250 psi. 
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Figure 130, Comparison of Liquid Metal Flows from 
Impinging Doublet Elements. 


A. Cerrosafe, 0.060— inch orifices, ^ P = 100 psi, 

Field of View: 2" x 2". 

B. Lithium, 0.043-inch orifices, = 250 psi, 

Field of View: 1" x 1". 
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Figure 132. 


Schematic of Hydrogen 
in Both Main and Gas 


Flow System, Showing Parallel 
Generator Supply Branches. 
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Circuits 
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INSTRUMENTATION 


FLUORINE -HY0R06EN SYSTEMS 


(7) FLUORINE TANK PRESSURE 
(?) FLUORINE TANK TEMPERATURE 
(?) FLUORINE SUPPLY TEMPERATURE 
(?) FLUORINE flowrate 
(?) HYDROGEN SUPPLY PRESCJRE 
(?) VENTURI INLET PRESSURE 
(?) VENTURI INLET TEMPERATURE 
(?) HYDROGEN MANIFOLD PRESSURE 
(?) FLUORINE MANIFOLD PRESSURE 
(lO) LOAD CELL 
(Tj) G.G. CHAMBER PRESSURE 
(f|) 6.G. CHAMBER TEMPERATURE 
(i?) CHAMBER WALL TEMPERATURE 
(w) CHAMBER PRESSURE 


Figure 133. Schematic of Fluorine and Hydrogen Flow Systems, Showing 
Instrumentation Locations. 
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Li SYSTEM PURGE (HOT HELIUM) 



TIME, SEC 

Figure 13^. Typical Firing Sequence, Showing Timing of Propellant Flows and of Lithium Purg 










5AA33-9/6/67-SI 


Figure 13t), Hydrogen Injector for Main Hydrogen Flow. This Model Has A 
Graphite Core, in Which the Orifices Are Located, and a 
Copper M' r. if old. 
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L *20 INCHES 
Ha ADDED AT 5 INCHES 
H a ADOED * 25 PERCENT 

Djo s 30 MICRONS 


Figure 1.39. Predicted Variation of c* Efficiency with F^Li 
Mixture Ratio at Indicated Conditions. 
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Figure 1 M> . Experimental Data (Run No. «» , Superposed on Predicted 
Indicated Lithium Mean Droplet Diameters. 
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Figure V»3. Experimental Data Superposed on Curves of Predicted c* 
Efficiency as a Function of Length of Fjj/Li Combustion 
Chamber for Three Overall Chamber Lengths, at Indicated 
Conditions, 






APPENDIX A 


CALCULATION OF CORRECTED C* EFFICIENCIES 

INTRODUCTION 

The performance indices used in the present program were corrected c* 
efficiencies obtained independently from measurements of chamber pressure 
and of thrust. Details of the computational procedures and of the correc- 
tions applied are given in this appendix. 


CORRECTED CHARACTERISTIC VELOCITY EFFICIENCY BASED ON CHAMBER PRESSURE 


Characteristic velocity efficiency based on chamber pressure is defined 
by the following equation: 




( P o>n <Veff Sc 

(°*) theo 


(A-l) 


where 


(P ) 

C O 

( A t^eff 


stagnation pressure at the throat, psia 

. 2 

effective thermodynamic throat area, in. 

2 

conversion factor (32. 17** lbm-f t/lbf-sec ) 

total propellant weight flowrate, lbm/sec 

theoretical characteristic velocity based on shifting 
equilibrium, ft /sec 


Values calculated from Eq . A-l are referred to as "corrected" c* effici- 
encies because the factors involved are not measured directly, but are 
obtained by application of suitable corrections to measured parameters. 
Thus, stagnation pressure at the throat is obtained from measured static 


A-l 



pressure at, start of nozzle convergence by assumption of isentropic ex- 
pansion, effective throat area is estimated from measured geometric area 
by allowing for radius changes during firing and for nonunitv discharge 
coefficient, and chamber pressure is corrected to allow for energy losses 
from combustion gases to the chamber wail by heat transfer and friction. 
Equation A-l may therefore be written as follows: 


V, 
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P C A t K c f p f TR f I)IS f FR f IIL 


f^iheo 


(A-2) 


where 


A. 


w r 


T 

(c*) 

f p 

f TR 

f 

MS 

‘fr 

f HL 


thoo 


= measured static pressure at start of nozzle convergence, 
psia 

2 

- measured geometric throat area, in. ~ 
conversion factor (32. 1?** Ibm-ft/lbf-sec ) 

= total propellant flowrate, lb/sec 

= theoretical c* based on shifting equilibrium 
calculations, ft sec 

- factor correcting observed static pressure to 
throat stagnation pressure 

= factor correcting for change in throat radius 

= factor correcting throat area for effective discharge 
coefficient 

= factor correcting measured chamber pressure for fric- 
tional drag of combustion gases at chamber wall 

- factor correcting measured chamber pressure for heat 
losses from combustion gases to chamber wall. 


Methods of estimation of the various correction factors are described in 
the following paragraphs. 


Pressure Correction (fp) 


Measured static pressure at start of nozzle convergence was converted to 

stagnation pressure at the throat by assumption of no combustion in the 

A-2 


nozzle and application of the isentropic flow equations. The shifting 
equilibrium specific heat ratio (y) ranged from 1.20 to 1.26 for the 
experimental firings reported herein. At the applicable contraction 
ratio (5.95), the corresponding stagnation static pressure ratios do not 
vary significantly; the same correction factor (f ]} = 1.006) was therefore 
used for all of the firings. 


Throat Radius Correction 


( f TR^ 


Temperature gradients produced in an uncooled nozzle wall by flow o» hot 
combustion gases result in thermal stresses which affect throat radius. 
Consequently, the geometric throat diameter measured in an ambient- 
temperature nozzle is not the same as that which exists during firing. 

When firing begins, thermal penetration of the nozzle wall is small with 
respect to the wall thickness, hence the outer wall diameter is unchanged. 
The inner-wall material will therefore expand toward the center, resulting 
in a decrease in throat diameter. As heat penetrates throughout the 
nozzle wall, the outer diameter will also increase, allowing outward ex- 
pansion of the inner portion and consequent increase in throat diameter. 
Hence, throat area during firing is a function of time, as well as of the 
physical properties of the throat material and the temperature and pressure 
of the combustion gases. 

Calculations of throat radius change start with the transient thermal stress 
equations for a hollow cylinder (Ref. A-l ) . Various analyses based on 
these equations have been made. One which is applicable to relatively 
short-duration firings (Ref. A— 2) assumes zero plastic stress, zero stress 
in the axial direction, a quadratic temperature distribution through the 
wall, circular throat and nozzle, and zero temperature and radius changes 
at the outer wall. The result is expressed as change in throat radius: 
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A-3 


where 


R. = inner radius of wall 
1 

R = outer radius of wall 
o 

OL = coefficient of linear expansion of wall material 
H = Poisson’s ratio of wall material 
6. = temperature at inner wall 
0 q = temperature at outer wall 

Although graphite is not strictly an isotropic material, directional 
property variations are not large enough to significantly prevent appli- 
cation of Eq. A-3 . Inner-wall temperatures were estimated by the method 
of Ref. A-3. The throat area correction factor due to thermal shrinkage 
varied from 0.993 (at 40-percent addition) to 0.989 (at 15-percent Hg 
addition) . 


Throat Discharge Coefficient Correction (^ DI g) 

The discharge coefficient is defined as the ratio of actual flowrate 
through the throat to the theoretical maximum based on the geometric 
throat area and one-dimensional, inviscid flow. Values of the discharge 
coefficient may be estimated either analytically or from correlations of 
the results of experimental studies of gas flow through nozzles. In the 
present program the ratio of the upstream wall radius at the throat to 
the throat radius was large (4:l) so that the discharge coefficient was 
very nearly unity and either of the two methods could be employed. 

In a recent critical study of available theoretical analyses of nozzle 
flow (Rei. A-4), values of the throat discharge coefficient were presented 
as functions of nozzle wall/throat radius ratio and throat Reynolds number^ 
For the applicable radius ratio (4) and Reynolds number (l x 10 to 2 x 10 ), 


A-4 


the indicated value is 0.995. On the other hand, from correlations of 
experimental conical nozzle discharge coefficients obtained by various 
investigators (Ref. A-5) and the geometry of the present nozzle, a value 
of 0.993 is obtained. Agreement between the two values is good and the 
average was used as the correction factor: ” 0-994. 


Frictional Drag Correction (fp^) 

Calculations of c* based on chamber pressure are concerned with chamber 
phenomena up to the nozzle throat. Drag forces to this point are gen- 
erally small enough to be considered negligible, particularly when the 
contraction ratio is fairly large as in the present case. Hence, the 

factor f__, was taken to be unity. 
rK 


Heat Loss Correction (fj^) 

Heat transfer from the combustion gases to the walls of an uncooled thrust 
chamber results in loss of enthalpy and thus decreases chamber pressure 
and thrust. This enthalpy loss is substantially reduced in an ablative 
chamber and is effectively recovered in a chamber cooled regeneratively 
by one of the propellants , whose initial enthalpy is raised by the heat 
absorbed. To obtain a true indication of performance efficiency in an 
uncooled chamber, measured chamber pressure must be corrected by a factor 
which accounts for heat loss to the walls and to the injector face. 

The three modes of heat transfer considered in estimating the heat loss 
rate to the thrust chamber walls were radiation, convection, and condensa- 
tion. These heat transfer mechanisms were assumed to be uncoupled; other 
assumptions regarding the heat transfer modes were based on posttest ob- 
servations of engine hardware. 

Three separate regions are distinguishable in the experimental thrust 
chamber: (A) the gas generator section, (B ) the lithium/fluorine combustion 


section upstream of the hydrogen addition point, and (c) the chamber 
section downstream of the hydrogen injector. 

In region A temperatures am! flowrates are low enough so that heat losses 
to the chamber walls are negligibly small. 

Combustion of lithium and fluorine in region B produces gas temperatures 
on the order of 9000 to 9800 R. In this region radiation from the gas 
cloud is by far the dominant mode of heat transfer. The meager emission 
data available for LiF and HI .Indicate that for the chamber geometry and 
pressure involved, the gas radiates as a continuum black body and that beam 
length is much shorter than any chamber dimension. Assuming unity view 
factor and black body emission, estimated radiation heat flux is approx- 
imately 26 Btu/in . 2 /sec . 

Confirmation of the order of magnitude of the heat flux in region B was 
obtained in the first Fg/Li/H,, firing, in which an all-copper hydrogen 
injector was used. Postrun, the core of the injector was found to be 
burned away (graphite cores were therefore used in all subsequent firings). 
An estimate of the heat flux required to burn the copper core was made 
on the basis^of its geometry, the cooling capability of the hydrogen flow- 
ing through it, and the time to failure. The value thus calculated was 
29 Btu/in. 2 /sec, which agrees with the estimate of radiation heat flux. 
Apparently total heat flux in region B (radiation plus convection) re- 
mains relatively constant at the indicated high level. 

Heat flax in region C was estimated on the basis of response data obtained 
from thermocouples placed in the graphite liner one-half inch from the 
inner wall. The rates of temperature increase indicated by these thermo- 
couples were used with Heisler transient conduction heat transfer charts 
to estimate the heat flux into the wall, which was found to be relatively 
constant at approximately 1.6 Btu/in . 2 /sec . This value derives essen- 
tially from convection since radiation io negligible in this region and 


m 


contributions from condendation are snail (observed LiP deposits on the 
wall were 0.01 to 0.03 inch thick). 

Heat loss in the nozzle region was estimated from the Bartz equation 
(Ref.A-6); average heat flux in the converging section was approximately 
2 Btu/ in . /sec . 


m. 


For each firing, heat losses to the chamber walls were computed for the 
various regions of the thrust chamber from the heat flux estimates given 
above. Total heat loss affecting the observed measurements of chamber 
pressure was then obtained by summation over the appropriate areas: 

ft = Z (q/A)A (A-4) 

where 

7+ = total heat loss 

q/A = estimated heat flux over given region 

A = area applicable to each q/A value 

The enthalpy of the injected propellants was reduced by an amount equivalent 
to the total heat loss to the chamber walls, and performance calculations 
were made to determine theoretical values of c* under the modified condi- 
tions. While this procedure lumps all heat lost and confines it to the 
injection station instead of gradual loss along the chamber length, the 
difference is not significant as far as this application is concerned. 
Comprrisr r » of the modified values of c* with the corresponding original 
values provided a measure of the applicable correction factor, f HL . Since 
heat f i ux in the Pg/Li combustion chamber represents by far the largest 
heat loss, the magnitude of the correction factor is a direct function 
of the length of that chamber. In the present series of data-producing 
firings, this length was either 20 inches (runs No. 6 and 7) or 5 inches 
(runs No. 8 and 9). With the 20-inch chamber, f^ was ~ 1.07 ; with 

the 5-inch chamber, f^ was ~ 1.02. 
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CORRECTED C* EFFICIENCY BASED ON THRUST 


Specific impulse efficiency is defined by the following equation: 


rj 


I 

9 


< P vac> 

(*.,) TO 


a 7 theo 


(A-5) 


where 


P 

vac 


P 


measured thrust corrected to vacuum conditions by the 

equation: F 53 F + P A 

^ vac a e 

measured thrust, pounds 




theo 


ambient pressure, psia 

area of nozzle exit, in.^ 

total propellant flowrate, lbm/sec 

theoretical vacuum shifting-equilibrium specific impulse 
at applicable expansion ratio, lbf/lbm/sec 


By application of suitable corrections to measurements of thrust made at 
sea level, corrected values of vacuum thrust may be obtained. With these 
values, which include allowances for all important departures from ideal 
flow, theoretical thrust coefficients may be used for calculation of c* u 
That is, Cj, efficiency is 100 percent if there is no combustion in the 
nozzle, if chemical equilibrium is maintained in the nozzle expansion pro- 
cess, and if energy losses from the combustion products are accounted for. 


Applicable corrections to measured thrust are specified in the iollowing 
equation: 


H 

e* 


< F * P » A ,> %V «k 


(w T ) (I ) 


s ' t heo 


(A 6) 


where 

p = measured thrust, pounds 

P * ambient pressure, psia 

a 
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A = area of nozzle exit, in.*" 

e 

w = total propellant flowrate, lbm/sec 

( I g ) th eQ “ theoretical vacuum shifting-equilibrium specific impulse 
L i€ ° at applicable expansion ratio, lbf/lbm/sec 

ip = correction for frictional losses 

FR 

<0 = correction for nozzle divergence 

DTV 

<0 « correction for heat losses to chamber and nozzle walls 

HL 

The correction factors in Eq. A-6 were applied to vacuum thrust (F + P a A 0 ) 
instead of to measured site thrust (f) because , for convenience, the cor- 
rection factors were calculated as changes in efficiency based on theo- 
retical vacuum parameters so that the total correction was of the form 


Correction for Frictional Drag (^p^) 


This factor corrects for the energy losses caused oy arug forces resulting 
from the viscous action of the combustion gases on the thrust chamber 
walls. Its magnitude, which is the integral of the local friction forces 
over the chamber inside wall, was estimated by a boundary layer analysis 
utilizing the integral momentum equation for turbulent flow. This analy- 
sis accounts for boundary layer effects from the injector to the nozzle 
exit by suitable description of the boundary layer profile and local skin 
friction coefficient. A computer program was used to carry out a numeri- 
cal integration of the equation including effects of pressure gradient, 
heat transfer, and surface roughness. The program required a potential 
core solution of th rt nozzle flow which was obtained from the variable- 
property, axi symme t r ic method of characteristics calculation of the flow 
field outside the boundary layer; corresponding properties for the sub- 
sonic combustion chamber flow field were also calculated. 


Computed values of <0^ ranged from 1.007 to 1.009 for the conditions of 
the experimental firings. 
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Correction for Nozzle Divergence (V^jy) 

The one-dimensional theoretical performance calculations assume that flow 
at the nozzle exit is uniform and parallel to the nozzle axis. The correc 
tion factor, allow 8 for nozzle divergence (i.e., for nonaxial flow) 

and for nonuniformity across the nozzle exit plane. It was calculated 
from the standard equation (Ref. A-7) : 

«W ’ [ I (l + c ° 9 a) ] (A ' 7> 

where a is the expansion angle (l5 degrees in the present case). This 
gives = 1.017. 


Correction for He at Losses (<^^) 

The effect of heat loss to the chamber walls on measured I g was determined 
by including this loss in the calculation of the theoretical values, as 
described above for the factor f^. Magnitude of ^ is directly dependent 
upon the length of the PgA. combustion chamber; for the 5-inch chamber, 

<p was about 1.02, and for the 20- inch chamber, was about 1.08. 


Sunmary 


The correction factors applied are sunmarized in the following listing. 


Applicable to 


M 


Applicable to 


M 


f 

- 

1.006 



p 





f 

TR 


0.9«9 

to 

0.993 

f „ 



0 . <M 



DIS 





f„. 

__ 

1.000 



FR 





‘hl 


1.017 

to 

1.076 

g FR 

= 

1.007 

to 

1.000 

°div 

= 

1.017 



^1L 

Tt 

1.019 

to 

1.078 


A-10 




REFERENCES 

A-l. Timoshenko, S. and J. N. Goodier, Theory of Elasticity , McGrav-IIi 1 1 
Book Company, New York, 1951* 

A-2. TIU65-107, Performance Characteriatics of Compound A/Hydrazine Pro- 

pellant Combination , Prepared for Air Force Rocket Propulsion Lab- 
oratory, Edwards Air Force Base, California, by Rocketdyne, a Division 
of North American Aviation, Inc., Canoga Park, California, May 1965. 
CONFIDENTIAL. 

A-3. Hatch, J. E. , et al., Graphical Presentation of Difference Soultions 
for Transient Radial Heat Conduction in Hollow Cylinders with Heat 


Transfer at the Inner Radius and Finite Slabs With Heat Transfer at 


One Boundary , NASA Technical Report R-56, I960. 

A-4. Singer, H. A., Discharge Coefficients of Nozzles Having Sharp Throat 
Curvature , LAP 66-518, Rocketdyne, a Division of North American 
Aviation, Inc., Canoga Park, California, 18 October 1966. 

A-5. Arbit, H. A. and S. D. Clapp, Fluorine-Hydrogen Performance Evaluation , 

Phase I, Part I: Analysis, Design, and Demonstration of High-Performance 


Injectors for the Liquid Fluor ine-Gase6us 


NASA CR-5Ji978, Prepared for NASA by Rocketdyne, a Division of North 
American Aviation, Inc., Canoga Park, California, July 1966. 

A-6. Bartz, D. R. , "A Simple Equation for Rapid Estimation of Rocket Nozzle 
Convective Heat Transfer Coefficients," Jet Propulsion , 27 , ^9-51, 
January 1957. 

A-7. Sutton, G. P. , Rocket Propulsion Elements , Third Edition, John Wiley & 
Sons, Inc., New York, 1963, p. 65. 


A-ll/A-12 




APPENDIX D 


MEASUREMENT ANALYSIS PROGRAM 


INTRODUCTION 

The primary purpose of a measurement analysis program is to provide a 
function which relates observed sensor outputs to estimates of corre- 
sponding system inputs, together with quantitative indications of the 
precision of this conversion. The function and the precision estimates 
are established on the basis of sensor calibration history, that is, upon 
a sequence of periodic calibrations of the sensor and its associated meas- 
uring and recording system against known inputs. 


Because calibrations must of necessity be made at a time differing from 
the actual firing time by several hours to several days, the changes in 
random sensor error with time must be established. In the Random Walk 
measurement analysis program (Ref. B-l ) this is accomplished by assuming 
that the input-output ratio at a particular input level performs a random 
walk in time which has normal distribution and variance. It assumes also 
that there is a random measurement error in the observed datum which is 
independent of the random walk and which is also normally distributed. 
Mathematical foundations and development of the program are given in 
Ref. B-2 and B-3 . 

On the basis of the sequence of periodic calibrations, the Random Walk 
program provides the following: 

1. A function, either linear or cubic, which converts observed 
system outputs into estimates of true system inputs 


B-l 


2. Coefficients of short-term and random walk variations as well 
as a combined value valid at specified times 

3. A decision based upon the calculated coefficient of variation 
and a prespecified imprecision limit as to whether the sensor 
should be used as is, recalibrated immediately, or discarded, 
and the maximum allowable interval to next calibration 


MEASUREMENT PROGRAM OUTPUT 

A typical Random Walk computer program output is shown in Table B-l . The 
first line of output gives the test stand name and number (Yoke, 0018), 
recording system (Beckman), transducer serial number (1641002), range 
(200 psi), ID number for data cards (018049), and the physical parameter 
being calibrated (P - 1A) . 

The next set of numbers ("Latest Output") is the most recent raw calibra- 
tion data. On the left are the readings (in Beckman counts) for the 
listed calibration input steps ("Input"); on the right are the precali- 
brate throw zero (Zl), the calibrate throw reading (CT), the postthrow 
zero (Z2), the precalibration zero (Z3), the postcalibration zero (Z4 ) , 
and the date of calibration ("Time"). 


The first two zeros (Zl and Z2) are averaged and subtracted from the throw 
to get a reduced throw. For each calibration step, a linear interpolation 
is made between the last two zeros (Z3 and Z4) and the interpolated result 
is subtracted from the reading to get a reduced reading. Each reduced 
reading is then divided by the reduced throw to get a scaled output. All 
scaled output values from all calibrations in the system history are then 
listed ("Scaled Output") under the appropriate input pressures, with one 
calibration per line and its date ("Time") at the right of each line. 


(I 
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The next line gives the result of the second test, which checks whether 
or not the input-output model is consistent with the estimate of O m 
(the root-mean-square estimate for the calibration curve fit andff^ 
should be approximately equal). If it is, then the model is labeled 
"SATISFACTORY": if not, the model is labeled "UNSATISFACTORY," indicating 
a significant intercept or an error in the input data. 


The following item indicates the ability of the system to meet the speci- 
fied imprecision requirement. On the basis of the calibration data, three 
situations are recognized: 

1. The system can never meet required precision, and should be 
replaced ; 

2. The system will fail the requirement within the next 2 days 
and should be recalibrated immediately; or 

3. The system will meet the requirement up to a certain date (30 
days maximum), on or before which it should be recalibrated. 

In this case, the estimated data reduction imprecision is given 
for test data taken 2 days after the most recent calibration 
and on the specified recalibration date. 


The final item is a 2 by 2 matrix, denoted by R, which is used to esti 
mate data reduction imprecision at any other time of interest and for 
any scaled output by the following expression. 


P 


V + S 2 (h a 2 + <r m 2 ) 


l/2 


(B-l) 


where 


P = estimated standard deviation for a reduced datum 
9 = scaled output 
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h 


V 



= number of days after most recent calibration 
V = matrix product: (s, s^) R ( S 3 ) 

Application of the results of the measurement analysis program tc esti- 
mation of random experimental errors and to measurement reliability is 
given in Appendix C. 
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TABLE B-l 



TYPICAL COMPUTER OUTPUT FOR TRANSDUCER CALIBRATIONS U 

MEASUREMENT ANALYSIS PROGRAM 


* YOKE 0018 BKM 1641002 200 018049 

* LATEST OUTPUT 


* 


* 

364 

789 

1103 

1633 

2162 

* 

INPUT 





* 

30 

70 

100 

150 

200 


SCALED 

OUTPUT 





0.1577 

C. 3706 

0.5278 

0.7932 

1.0581 


0.1577 

0.3689 

0.5265 

0. 7929 

1.0572 


0. 1567 

0.3704 

0.5280 

0. 7923 

1.0566 


0.1556 

0.3682 

0.5283 

0.7919 

1.0565 


0.1551 

0.3692 

0.5283 

0.7924 

1.0580 


0.1570 

C. 370 1 

0.5276 

0. 7927 

1.0578 


* MEASUREMENT VARIANCE IN INPUT-TO-SCALED OUTPUT RATIO 

* RANDOM WALK VARIANCE IN INPUT-TO-SC ALEO OUTPUT RATIO 

* RATIO OF SHORT-TERM VARIANCE TO RANDOM WALK VARIANCE 

* COEFFICIENT OF SHORT-TERM VARIATION 

* COEFFICIENT OF RANDOM WALK VARIATION 

* REQUIREMENT FOR COEFFICIENT OF VARIATION OF REUUCEO DATA 

* SYSTEM NOW PASSES TEST FOR LINEARITY I TYPE I ERROR*. 05). 

* DATA REDUCTION FORMULA IS 

* (INPUT) * (1.8920E 02)*(SCALED OUTPUT) 


* ABCVE OUTPUT-INPUT MODEL IS SATISFACTORY (TYPE I ERROR*. 05 


♦ 

SYSTEM SHOULD BE CALIBRATED 

CN OR BEFORE 9-19-65 


♦ 


COEFFICIENT 

OF 

VAR I AT 

♦ 


COEFFICIENT 

OF 

VARI AT 

♦ 

OATA REDUCTION MATRIX 

( 2.40921E-C2 

0 . 


* 


( 



♦ 


( 0 . 

0 . 



£'C~/ 


CALIBRATIONS USING RANDOM WALK 
PROGRAM 

PC 1 A 

12 l 3 Z4 TIME 

49 49 49 8-20-65 


018049 

Z1 CT 

49 2046 


TIME 


8-20-65 

8-13-65 

7-27-65 

7-27-65 

7-27-65 

7-27-65 


00 
01 

0? (DAYS) 

(PERCENT) 
(PERCENT/0AY**.5I 
(PERCENT) 

OR*. 05). 


OUTPUT) 
ERROR*. 05). 


RATIO 

RATIO 

ANCE 


CEO OATA 


0.28773F- 

0.23258E 

0.12318E 

0.283 

0.0806 

1.500 


5 

Of VARIATION OF REDUCED DATA 
OF VARIATION OF REDUCED DATA 


ON 9-19-65 * 0.53 PERCENT 
ON 8-22-65 * 0.32 PERCENT 


0 . 

0 . 


) 

> 

I 
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APPENDIX C 


MEASUREMENT ERROR ANALYSIS 


INTRODUCTION 

Because it is not possible to measure the true value of any physical 
property or parameter, the error limits, or uncertainty interval, asso- 
ciated with any experimental measurement must be specified. It is the 
purpose of this appendix to indicate the reliability of the experimental 
determinations of performance efficiency by estimation of the errors 
inherent in the data acquisition processes and in the calculation pro- 
cedures. This permits estimations to be made of the range within which, 
at a given confidence level, the true values of the measured or calculated 
parameters may be expected to fall. 

If error be identified with departure of an experimental measurement f*om 
the "true" value, its magnitude can never be completely known; if it were 
known, it would become a correction which could be systematically applied. 
Hence, error limits can only be stated within probability limits. The 
estimation is made by an error analysis procedure which, in the present 
application, consists of the following steps: 


Estimation of the uncertainty intervals in the individual trans- 
ducers including the measuring systems in which they are used 

Combination of the uncertainty intervals of duplicate or redundant 
sensors into an uncertainty interval for the measurement 

Combination of the uncertainty intervals of several measurements 
pressure, temperature, and flowmeter frequency) into an 
uncertainty interval for the parameter they determine (e.g., 
flowrate) 


4k 


Combination of the uncertainty intervals of the measurements 
entering into calculation of the value of the desired variable 
(e.g., characteristic velocity) to estimate the uncertainty 
xnterval of the calculated result 

Two types of errors are possible for any measurement: 

1. Systematic errors . These are associated with the particular 
system, with the experimental techniques employed, or w'ith the 
calibration procedures. They cannot be estimated by statistical 
methods, and are minimized primarily by careful calibration with 
the best available standards, by requirements for consistency and 
traceability of the experimental and calibration techniques, and 
by critical examination of experimental data. 

2. Random errors . These arise from unpredictable and unknown 
variations in the experimental situation and are generally assumed 

to follow a normal distribution to permit simple statistical analysis. 
Error analysis is concerned only with random errors and implicitly 
assumes that systematic errors can be eliminated in a carefully con- 
ducted experimental program. 

From the properties of the normal, or Gaussian, distribution 
function, the probability of a system error exceeding ±1 times 
the standard deviation (cr) is about yi percent, the probabi lity 
of exceeding ±2 times the standard deviation is about 4.6 percent, 
and the probability of exceeding ± 3 times the standard deviation 
is about 0.3 percent.. The value of y in a result expressed as 
(x t y) is generally taken as 2a , thus setting the confidence 
level at 93 percent so that the estimated odds of the true value 
of the result falling within the range (x - y) to (x + y) are 
20 : 1 . 



SENSOR PRECISION 


The precision of a measurement obtained as the output of a physical instru- 
ment or sensor is a quantitative estimate of the uncertainty associated with 
that measurement. This estimate is made by statistical analysis of the 
outputs of the sensor when repeatedly acted upon by known inputs. By sensor 
is meant not only the transducer itself but the complete system which con- 
verts the transducer signal to a numerical value of its physical parameter 
analog. The known inputs, of course, have uncertainty limits of their own, 
but for practical purposes it is assumed that they are accurate (i.e., 
identical to true values) within the limits required by the experimental 
situation. Ultimately, these inputs must be directly traceable to estab- 
lished standards, such as those of the National Bureau of Standards. 

When a sensor is calibrated against known inputs, precision may be con- 
sidered as the certification of an error band within the calibrated 
interval and within a given confidence level. Thus it provides a measure 
of "closeness to truth'* of the reduced data. Precision may be numerically 
expressed as the standard deviation of a measurement, which has the same 
units as the measurement itself, or as the coefficient of variation, which 
permits valid comparisons between measurements in different units. Coef- 
ficient of variation (C y ) is the standard deviation expressed as a percentage 
of the mean, thus making it dimensionless: 



where 

C7 = standard deviation 
(i = sample mean value 
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Pressure 


The coefficients of variation of tne pressure transducers were obtained 
by application of the Random Walk measurements analysis program to the 
calibration data, as discussed in Appendix B. The values obtained ranged 
from 0.1 to OA percent, for static calibration made on a pressure manifold 
mounted on the thrust stand. 

Other errors in pressure measurements may arise, in addition to the ran- 
dom, statistical uncertainty limits. Thus in the measurement of chamber 
pressures through a drilled wall tap, erroneous values of stream pressure 
may be indicated because of the effect of the hole itself upon the flow. 

The following estimated magnitudes (by percentage) of this tap error, 
which is a function of stream velocity, are based on experimental data 
obtained with water and gas (Ref. C— i) : 

(Mach No. as 0): 0.00 

(Mach No. * 0.3): 0.05 

(Mach No. s* 1.5): 1-30 

Coupling errors, arising from effects of the tubing joining the pressure 
taps to the transducers, are not significant in the present series of 
firings since precise dynamic or transient response pressure measurements 
were not required (ttef.C-2) and tubing lengths were not greater than 3 feet. 


Thrust 


Values of coefficient of variation obtained by application of the Random 
Walk measurement analysis program to thrust calibrations were in the range 
0.3 to 0.8 percent. A possible source of error in the thrust measurements 
arose from the necessity of taking system prerun and poatrun zeros with the 
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same degree of inlet line chill and LLi inlet line heating that 

existed during the firings. On the basis of thrust calibrations made 
with inlet lines at varying temperatures (Hef. C-2), the estimated C y 
value due to variation in zero readings is 0.1 percent. 


Throat Area 

Geometric throat diameter was measured with an expansion micrometer by two 
observers before and after every firing. Maximum coefficient of variation 
of the calculated areas was 0.3 percent. 


Volumetric Flowrate 


The coefficients of variation of the turbine flowmeters used to measure 
LFg flowrate were determined from calibration data. Observed C y values, 
which refer only to flow bench water calibrations, were 0.1 percent for 
both meters. Corrections for thermal and viscosity effects in converting 
these calibrations to cryogenic LFg factors are discussed in another 
section of this report. In addition, however, there are unpredictable 
water-to-cryogenic calibration shifts (Ref. C-3) which introduce addi- 
tional sources of error. The coefficient of variation arising from this 
source is appioximately 0.3 percent (Ref. C-4). 


Estimated C y value for the magnetic flowmeter used for liquid lithium is 
0.5, based on several calibrations of the e lectromotive force (emf) output. 


Temperature 


Resistance Temperature Sensors . The platinum resistance thermometers were 
precision calibrated by the manufacturer. These calibrations were checked 
by taking several emf readings with the sensors immersed in LNg end in U)g 
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at atmospheric pressure; these were correct within the limits of reada- 
bility with <0.02 percent coefficient of variation in a series of four 
measurements. Root-sum-square (RSS) error limits of these sensors based 
on specifications for repeatability, insulation, time lag, friction heating, 
and interchangeability were 0.1 percent. Voltage readout of the transducers 
was adjusted to calibration values by means of a standard decade resistance 
box with error limits of 0.2 percent. 


Thermocouple s « Iron-constantan thermocouples were used to measure tempera- 
tures of GH^ in the venturi plena, and chrome 1-alumel thermocouples were 
used to measure lithium system temperatures. Because the latter were not 
involved in performance measurements, they will not be considered in this 
section. Estimated error to be expected with new iron-constantan thermo- 
couple wire at ambient temperatures is 0.7 percent (Kef. C-l). Thermocouple 
calibrations were electrical only; i.e., the emf readout was adjusted on the 
assumption that thermocouple-generated electromotive forces correspond to 

standard values. Total estimated C is 1.0 percent. 

v r 


COMBINED ERROR ESTIMATION 


Redundant Measurements 


Two independent transducers were used to measure most of the important 
parameters (except thrust) in order to increase measurement reliability. 
The most probable value of a redundant measurement is a weighted average 
in which the weight (VT) assigned to an individual determination is given 
by 


1 

a 2 


( 02 ) 


where <j ^ is the standard deviation associated with the i tn measurement and 
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n is the variance. 
x 


The variance of the weighted mean, O m » is given by 



Clearly the variance of the weighted mean is less than any of the indi- 
vidual variances; for example, in the particular case of two measurements 
with equal variances, the variance of the mean is half the individual 
variances. 


Combined Measurements 


The standard deviation of a parameter which is a function of two or more 
independent measurements is taken as the RSS of the standard deviations 
of the independent measurements. Thus LF^ flowrate is a function of flow- 
meter frequency and fluorine temperature (assuming no significant error 
in conversion of fluorine temperature to equivalent density): 


where 

f = flowmeter frequency 

T = oxidizer temperature 

o 

The standard deviation of the oxidizer flowrate is then: 


w 


o 



(C-S) 


In the same way the standard deviation of hydrogen flowrate, which is a 
function of measured pressure and temperature in the venturi plenum, is 
given by 



(C-6 ) 
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where 


<7p = standard deviation of GH^, pressure measurement 

“2 

<7 = standard deviation of GH 2 temperature measurement 

H 2 

Standard deviation is converted to coefficient of variation by Eq. C-l 


When several measured variables are combined algebraically to yield an 
experimental result, the standard deviation of the result, which takes 
into account the propagation of the individual errors, is given by 


a R l 3x x a l 


\ 5x 2 2 


a 2 j + # * ‘ + \*x °n 


where 


= standard deviation of calculated result 


x 1 » x 2 > • 


,x = measured variables 
n 


= <P (* 1 »x J2 . • • > X n ) 

# a = standard deviations of x u x 0 , 
n l*: 


. x q , respectively 


The following example illustrates the application of this type of error 
analysis: 

Run No. 8 (Step l) 


Chamber Pressure : CJp -0.46 psia 

c 

Calculated by Eq. C-l and C-3 from calibration values (0.28 and 0.13 
percent) of the redundant sensors. 




Throat Area: nr A = 0.03 in. w 

A t 

Cy leu lilted from eight measurements of throat diameter. 

Thrust : a p = 1 Ibf 

Calculated from HSS of calibration C y (O.bl percent) and zero point 
variation (0.10 percent). 

Propellant Flowrate : a ■ = 0.04 lb/sec 

W T 

Calculated from the coefficient of variation of the individual i lowrates: 

C v (*LF„) “ °’ 35 percent ’ C v (* GHy ) = 0,31 P ercent * C v (*LLi) = 05 ° P ercent * 
Correction Factors : F = 0.003 

Obtained from HSS of estimated variances of the individual corrections. 


c*, Based on P c : ft/sec 


Obtained by application of Eq. C-7 to the expression 

P c A t , c (C.F.) 

C * = — 

*T 

where (C.F.) is the net correction factor. 


(C-8) 


The resulting expression is: 


‘) 

A t g c (C.F.) 

2 

-c «c (t ' F ) , 


a P 

. W T C J 


A 

W T 1 J 


>■, \ *c « F ‘> 


nr. 


w 


P A g 
c t c 


J C ,F. 


Substitution of numerical values gives (? c# . 


(C-9) 


c-9 




► 


At the 95-percent (2a) confidence level, this value of a c * corresponds 
to an error band of approximately ±1.0 percent. 

I Based on F: a T =2.2 seconds 

s I 

s 

Obtained by application of Eq .C-7 to the expression 


I 

s 


F (C.F.) 


w 


T 


where again (C.F.) is the net correction factor. 


(c-io) 


The resulting expression is: 


2 

F 

2 

C.F. 

2 

F (C.F.) 

°I = 
s 

* T °C.F. 

+ 

w ^F 

T ‘ 


2 

w T W T 


(C-ll) 


Substitution of numerical values gives . At the 95-percent (2 <t) confi- 

dence level, the error band is observed f S and therefore in c* based on 

s 

thrust, is approximately ±1.3 percent. 


DYNAMIC PRECISION 

The estimates of expected standard deviations calculated above are based 
on static calibrations of pressure and thrust sensors, and hence may not 
be strictly applicable to the dynamic system represented by a firing 
thrust motor. It is generally assumed, however, that such calibration 
data may be extended without significant change to dynamic systems 
oscillating at very low frequencies and amplitudes, and that steady- 
state stable combustion is such a system. 


C-10 


Ul illitii li 
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